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Radiant Interchange Within an Enclosure 


R. V. DUNKLE? | Part 1 Absorption and Emission Behavior of Gases 


mi IS THE First of three papers relating to a general 
analysis of the radiant interchange within an enclosure filled with 
an absorbing and emitting medium, i.e., combustion products. 
The problem is approached initially by reviewing the basic 
physics of the absorption (or emission) of thermal radiation by 
gases (Part 1). The band approximations of Elsasser [6]* are ex- 
amined and compared with the experimental results of Howard, 
et al. [12]. The inapplicability of the exponential absorption law 
on a band energy or total energy basis is briefly discussed and the 
deviations from experiment pointed out. As an engineering ap- 
proximation to the use of a monochromatic energy analysis, the 
absorption of gases on the band energy basis is incorporated in the 
geometrical interchange relationships between two surfaces. 

In Part 2, the network method of Oppenheim [16] is used to 
derive the basic network relationships between one node and the 
other nodes of an enclosure. This is first done on a monochro- 
matic basis and then reduced to the band energy approximation 
using the results of Part 1. The usual method of analysis, the 
gray radiation approximation, is shown to follow from the band 
energy and to be generally an approximation of lesser accuracy. 
A reiterative method of solving the network node equations for as 
many nodes as deemed practicable is presented in Part 3 and 
specific solutions are given as examples. These examples are 
used to illustrate the relative accuracies of the band energy and 
the gray radiation approximations. 

{n order to obtain a solution of the enclosure problem, a num- 
ber of assumptions have been required. These assumptions were 
necessitated by the inadequacy of present knowledge. Also, suf- 
ficient information does not exist regarding the radiative charac- 
teristics of materials to permit a numerical analysis of most en- 
closures. These gaps in analytical and physical knowledge are 
indicated and are promising areas for both theoretical and ex- 
perimental study. 


Basic Transfer Equations 


The general expression for the net attenuation of a beam of 
monochromatic energy in an absorbing and emitting medium is 
given by [6]: 


dl, = k{Iw — I,)dm (1) 


! Shell Development Company, Emeryville, Calif. Mem. ASME. 

? Associate Professor, Mechanical Engineering Department, Uni- 
versity of California, Berkeley, Calif. 

* Numbers in brackets designate References at end of Part 3. 

Contributed by the Heat Transfer Division of Tae American 
or Mecuanicat EnGrineers and presented at the ASME- 
AIChE Heat Transfer Conference, Storrs, Conn., August 9-12, 1959. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, April 29, 
1959. Paper No. 59—HT4. 
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This relation implies the use of Kirchhoff’s law equating the 
emission and absorption of the medium. A more familiar form 
of equation (1) may be obtained by assuming the emission of the 
medium to be zero: 


dl, = —k,I,dm (1a) 


The mass of the medium used in either of these expressions refers 
to the molecules which are active in the absorption or emission 
process at the frequency of the incident energy. This mass is re- 
lated to the concentration of the active molecules in the path by: 


dm = pdr (2) 


Equation (1) may be integrated for an isothermal homogeneous 
medium for the usual boundary conditions: 


m=0; 1, (3a) 
to give: 


I, = Iw[l — exp(—km)] + Lelexp (—km)] (3b) 


The restriction to isothermal homogeneous media cannot be re- 
moved until a specification is made of the temperature depen- 
dency of k,, m, and I,,. The latter is known since it is given by 
Planck’s law of ideal radiation; the first is a property of the 
medium; and the second is a variable of the particular problem. 
Specific conditions can be established which would permit the 
integration of this equation but a general form is not readily ob- 
tainable. Any approximate analysis of nonisothermal media 
must begin by satisfying the requirements of equation (1) or by 
ascertaining the magnitude and direction of the error by some 
method. The exponential factors appearing in equation (3b) are 
conveniently used to define the transmission and absorption (or 
emission) of the medium: 


T, = exp (—k,m) (4) 
a, = 1 — [exp(—k,m)] = 1-17, (5) 


As a first step in developing the radiation exchange equations, 
two surfaces separated by an absorbing and emitting medium will 
be considered (Fig. 1). The gas volume involved is that volume 
determined by the two areas and the outermost envelope formed 
by straight lines connecting the boundaries of the two areas. That 
is, for any point on surface one, the gas volume considered lies 
within the solid angle defined by surface two. The converse is 
true for any point on surface two. The intensity of radiation 
leaving an opaque surface consists of radiation emitted and re- 
flected by the surface. The irradiation of surface one from the 
direction of surface two is given by: 


Ai J As 


(6a) 
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ABSORBING 
MEDIUM 


Fig. 1 Geometry of radiation exchange between two surfaces with an 
absorbing medium 


Similarly for surface 2, the incident energy from the direction of 
surface 1 is: 


Ai J Ar JO 


(6b) 


The irradiation of a surface from the direction of a second surface 
in an enclosure filled with an emitting medium is given by equa- 
tion (6a) or (6b). This is composed of the radiation emitted by 
the medium and the transmitted radiation from the bounding 
surfaces. The emission from the gas is determined by the same 
geometry as the transmitted radiation. The radiation originating 
(emitted plus reflected) at either surface 1 or 2 and absorbed by 
the intervening medium is given by: 


cos 8, cos 8 
Tun = f f f + I) dA d Addy 
Ai JA: 


(6c) 


The net heat transfer to the medium by radiation is obtained by 
subtracting from equation (6c), the energy emitted by the medium 
and directly striking surface 1 or 2. The net energy exchange be- 
tween the two surfaces is the difference between equations (6a) 
and (6b): 


= f TALn = dA,dA;dv (7a) 
Ai An JO 


This equation should be clearly understood to apply only to the 
net radiation flux between two surfaces and only for black sur- 
faces does it represent the actual heat gained or lost by the surface. 
The radiosity of a surface is equal to 7 times the intensity for a 
diffuse surface, and for a black surface (emissivity of unity) the 
radiosity becomes identical with the emissive power. The re- 
mainder of the discussion will be restricted to diffuse surfaces, 
although many of the relations given will be equally valid for 
nondiffuse conditions. Equation (7a) may be rewritten with the 
use of the relation between monochromatic radiosity and mono- 
chromatic intensity: 


Ai J Ar JO 


2 / FEBRUARY 1960 


For a transmissivity of unity and black surfaces, i.e., the radiosity 
equals the monochromatic emissive power (£,), this last expres- 
sion becomes the conventional radiant interchange equation: 


= — T:*) 


(8a) 


= AFn = f f 
Ada ™ 


= f, ” Edy (8e) 


Theoretical Gas Absorptivities or Transmissivities 


Equation (3b) may be integrated over all frequencies to deter- 
mine an absorption or transmission in terms of the concentration 
of the absorbing medium lying along a straight line path. The 
resulting expression may be substituted in equation (7a) since 
the integration with respect to the frequency is independent of 
geometry. This integration may be written as a transmissivity 
or an absorptivity, thus: 


T= 


— exp( —k,m)dv] 
a= 


J dv 


The transmissivity and absorptivity of a medium are quantities 
relating the effect of the medium upon the radiation from an ex- 
ternal source. The expression for the emissivity of the medium is 
of the same form as that for the absorptivity but with the mono- 
chromatic emissive power corresponding to a Plankian radiator 
at the temperature of the medium substituted for the monocliro- 
matic radiosity of the external source. Therefore, in general, the 
total absorptivity and emissivity of a medium are only equal if 
the monochromatic radiosity of the source is equal to the mono- 
chromatic emissive power of a black body at the temperature of 
the medium at all frequencies where the medium absorbs radia- 
tion. The media of principal engineering interest are gases and 
the remaining discussion will be further restricted to these media. 

The emission or absorption of radiation by a gas does not ex- 
tend continuously over the whole spectrum as does that of a solid 
body. Instead, the emission or absorption of radiation is re- 
stricted to specific regions of the spectrum. A complete descrip- 
tion of molecular spectra is beyond the requirements of this paper 
and only an approximate description is necessary to demonstrate 
the structure of this type of spectra. 

The emission or absorption of energy by a gas is caused by 
changes in the energy levels of the molecules composing that gas, 
in accordance with Planck’s law: 


AW = hv 


(10) 


(11) 


The energy levels of a molecule may be written as the sum of the 
energies of the possible modes of oscillation: 

W=W,+W,+W, (12) 
The magnitudes of the electronic, vibrational, and rotational 
energies are such that as a generality, electronic transitions give 
rise to lines in the visible and shorter wave lengths, vibrational 
changes are in the infrared and rotational transitions are in the 
far infrared frequencies. However, the vibrational energy levels 
of a molecule are intimately coupled to the rotational levels and 
if the frequency of the energy has the correct value, a rotational 
transition may occur simultaneously with a vibrational transition. 
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Since the vibrational energies are much larger than the rotational 
energies, the frequency determined by the vibrational energy dif- 
ference is dominant and the rotational differences cause perturba- 
tions of this basic frequency. As a result, the spectrum of a gas 
in the region of a vibrational frequency will be composed of a 
number of absorption lines in the vicinity of this vibrational fre- 
quency. This is commonly referred to as a vibration-rotation 
band and the individual absorption or emission lines within the 
band are the rotation lines. If more detailed information is de- 
sired, reference should be made to texts on molecular spectra [9]. 

The emission or absorption of energy by molecules as a result of 
electronic energy changes in combination with vibrational and 
rotational energy changes (the rotation energy changes are not 
necessarily required) give rise to lines in the ultraviolet, visible, 
and near infrared. This is the band spectra observed from a gas 
in a highly excited state, e.g., in shock tubes or flames. These 
bands consist of many lines closely spaced but occur in regions for 
which instrumentation of high resolving power is available. 
Possibly as a consequence, the analyses of these bands are gen- 
erally restricted to the examination of the lines within the band 
relative to each other rather than to the band itelf. To the best 
of the authors’ knowledge, these bands have not been studied in a 
manner similar to the studies of the vibration-rotation bands 
which will be described. At the present time, treatment of these 
electronics bands in a manner analogous to the procedure to be 
described for vibration-rotation bands is not possible. This is 
due to the lack of experimental and analytical information re- 
garding the electronic bands. 

These prefacing remarks were necessary to partially justify 
definitions of total line and band absorptions. A total line ab- 
sorption would be for the total absorption of a single line, either 
rotation or electronic, and a total band absorption would apply 
to the total absorption of a vibrational-rotation band: 


av, Yell — exp(—kem) 
a, = (13) 
J dv 


J,(l — exp( —k,m) 
a, = —— (14) 
J dv 


Arp 


J,{exp( —k,m) 
i= (15) 
J 


ay, lexp( —kym) 
tp = (16) 
Arp Jv 


The total transmissivity and absorptivity of the gas [equations 
(9) and (10)] are given by: 


a J dv + J dv 


(17) 


1 
f, J dv =I ) 


(18) 
Inspection of the necessary equations will show the usual identi- 


ties associated with the absorptivity and transmissivity have been 
retained: 
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a,+t,=1 (19a) 
+ tz =1 (19) 
a+r=1 (19¢) 


Inclusion of both line and band absorptions in equations (17) 
and (18) permits the consideration of line and band effects which 
may occur in combustion and shock wave processes. For many 
engineering purposes, the importance of the line absorptions is 
negligible and only the influence of the bands must be considered. 

Absorption of a Single Isolated Line. The determination of the ab- 
sorption of a single line outside vibration-rotation band requires a 
specification of the variation of the absorption coefficient with 
frequency. This is generally accepted to be given by the Lorentz 
equation: 


6 
E — + (20) 


where 
S= f (21) 


Substitution of equation (20) into equation (13) gives an equation 
which cannot be integrated without a specification of the fre- 
quency. However, Elsasser [6] observed that the expression 
could be integrated by assuming: (a) The absorption of the line 
is essentially complete in the interval Ay, i.e., Av, >6, and (6) 
the intensity of the incident energy is constant over the interval 
Av,. The second assumption permits the monochromatic emis- 
sive power to be removed from the integral and yields an integral 
in closed analytic form. The result for a homogeneous medium 
is: 


a, = — —8)) (22) 
Sm 
B= 


where Jp and J; = Bessel functions of the zeroth and first order. 
This expression can be expanded to obtain values of the line ab- 
sorption for small and large values of [6]: 


Sm 
B<1; qu = (24a) 
2 
B>1; a, (Smé) (24b) 


The corresponding expressions for the transmission of a single line 
are: 


Av, 

wy B>1 (25e) 
Av, 


The latter expressions should be compared with usual assumption 
of exponential transmission of a gas. Only for the case of small 8 
{equation (25b)] would the two agree in mathematical form. 
The assumption of constant incident energy required by Elsas- 
ser to make this integration does not introduce serious error. The 
width of either an isolated rotational or a single electronic line is 
quite small, e.g., for a rotational line in the infrared, the width is 
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less than 0.01 micron. Consequently, the incident energy will 
not vary rapidly within such a small interval. 

The Absorption of a Band.‘ A vibration rotation band is com- 
posed of a number of rotation lines spaced about a central fre- 
quency and the intensity of an individual line decreases with the 
distance of the line from this central frequency. The intensity of 
any specific line depends upon: (a) The probability of the energy 
transition and (b) the fraction of the total number of molecules in 
the required energy state for this transition. The first is de- 
pendent upon the molecular configuration of the molecule and the 
second is a function of the thermodynamic state of the gas. 

The three gases of primary interest in engineering are carbon 
monoxide, carbon dioxide, and water vapor. The observed spec- 
tra of the first two are regular in appearance [2], whereas, the 
spectrum of the water molecule is quite irregular [19]. This ir- 
regularity arises from the overlapping of adjacent vibrational 
bands corresponding to two oscillation modes with central fre- 
quencies differing by less than the vibrational band widths. Asa 
consequence, the theoretical analyses of carbon monoxide and 
carbon dioxide have been approached in a different manner than 
the analysis of water vapor. 

The exact integration of the total absorption of carbon monox- 
ide has been performed by Penner [17] using numerical computa- 
tional techniques. Elsasser obtained an analytic band absorp- 
tion for carbon dioxide which is also applicable to carbon monox- 
ide, using an idealization of the spectrum [6]. The model con- 
sidered was a band composed of periodic equal intensity lines uni- 
formly spaced with frequency. This assumption permitted the 
absorption coefficient to be expressed as: 


S 


(v — nd)? + & 


(26) 


In addition to the assumption given, the assumption made by 
Elsasser for a single line regarding the incident energy intensity- 
frequency distribution was also used, i.e., constant incident energy 
over the band width. Substitution of equation (26) into equation 
(14) gives for a homogeneous medium: 


ag = sin hy [exp( sin hy coth (27) 


S 
(38) 


This gives for two particular cases: 


‘Sm)'/* 
d > O(lines far apart); a, = erf (29) 


d< &lines close together); a, = [: — exp (=) (30) 


The spectrum of the water molecule appears as a number of 
irregularly spaced absorption lines having an almost random in- 
tensity distribution [2]. Goody [8] observed from the experi- 
mental work of Randall, et al. [19] that the apparent random 
nature of the absorption was suitable for statistical analysis. 
Using this postulate and a Lorentz line shape [equation (20)], 
Goody obtained for the transmission of a water vapor band: 


Sméb 
d 


‘ The lines referred to in this section are the rotational lines which 
combine to make up the vibration-rotation band and should not be 
confused with a single rotation or electronic line discussed in the pre- 
ceding section. 


(31) 


4 / FEBRUARY 1960 


The integrated intensity (S) and the line spacing (d) are statistical 
quantities which must be determined from experiment and are 
not related to the physical geometry of the molecule by known 
functions. This analysis also assumed a constant incident energy 
within the frequency limits of the band. 

The analyses by Goody and Elsasser have been the basis for 
extensive elaborations and studies in meteorology. The recent 
work of Plass, et al. [18] exemplifies such an extension of the El- 
sasser model. Unfortunately for the application to engineering 
problems, these analyses do not offer additional refinement of the 
basic relations obtained. The work of Goody and Elsasser do 
furnish an opportunity to speculate qualitatively upon the effect 
of pressure and temperature on gaseous absorption. In both the 
Elsasser and Goody analysis, the half-width (6) of the lines ap- 
pears and the influence of temperature and pressure on this factor 
from simple collision theory for small concentrations may be 


written as [15]: 
T,\'* 
6 = ( (32) 


The reference half-width is dependent upon the nonactive mole- 
cules in the gas and upon the concentration of the active 
molecules, i.e., foreign broadening and self-broadening influences. 

The integrated intensity of the lines (S) appearing in equa- 
tions (28) and (31) are based upon an idealization of the spectra 
and are not related to the physical geometry of the molecules. 
However, if this factor may be considered to be representative of 
the true intensities in the manner of an average value, the effect 
of pressure and temperature may be considered qualitatively. 
As noted previously, the intensity of a rotational line in a vibra- 
tion-rotation band is dependent upon the friction of the total 
number of active molecules in a specific energy state. This 
fraction is determined by the partition function and is tempera- 
ture dependent [9,17]. Quantitatively, a description of the ef- 
fect of temperature upon the partition function is difficult. 
Qualitatively, the fraction of molecules having an-energy which 
gives rise to the lines in the wings of the band is increased at 
the expense of the lines in the center of the band as the tempera- 
ture is increased. This results in a broadening and flattening 
of the vibration-rotavion band, i.e., temperature broadening. 

The exact definition of the mass of the medium traversed by 
the radiation is given by equation (2) and is based upon the num- 
ber of molecules per unit volume (concentration) multiplied by 
the path length. The quantity often used in the literature to 
represent m neglects the effect of temperature upon the concen- 
tration and is the product of the partial pressure and the path 
length: 


w = pr (33) 


The preferred definition is given by equation (2) where p is ob-. 


tained from the ideal gas law: 


p 


= — 34 
p RT (34) 


Since all reported experimental measurements have been made 
under isothermal conditions, the temperature effect has been in- 
herently removed as a factor. The definition given by equation 
(33) will be used in the remainder of the paper, but the more exact 
definition should not be ignored if nonisothermal conditions are to 
be considered. 

Temperature affects the absorptivity of a band by changing 


5 The usual definition of the integrated intensity includes the total 
number of active molecules multiplied by this fraction. The total num- 
ber of molecules has been included asa basic factor in the definition 
of m. 
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the half-width of the rotational lines composing the vibrational- 
rotational band; by altering the fraction of absorbers at the 
various possible energy levels which changes the relative intensi- 
ties of the lines and the width of the band; and by affecting the 
concentration of the absorbers for a fixed partial pressure. The 
relative strength of two vibration bands is similarly affected by 
the temperature. The total pressure also affects the half-width 
of the rotational lines and for increased pressures may cause 
these lines to overlap. Another effect of total pressure is upon 
the number of absorbers in the radiation path for a fixed partial 
pressure. The magnitude and form of the corrections for band 
absorptivities necessary to account for these influences will re- 
quire experimental investigation. Hottel [10] gives a procedure 
for correcting the emissivity (absorptivity) of certain gases for 
temperatures which has been corroborated theoretically by Pen- 
ner [17]. However, the total emissivity includes the tempera- 
ture dependency of the energy distribution of a Planckian radia- 
tor and the various absorption bands acting in toto. Conse- 
quently, the correction by Hottel cannot be used, ad hoc, for band 
absorptivities. 


Experimental Gaseous Absorptivities 


No experimental determinations of the individual band ab- 
sorptions of carbon monoxide, carbon dioxide, and water vapor 
have been reported which are adequate for engineering use, al- 
though total emissivities have been available for a number of 
years [10]. A comprehensive study of the individual bands would 
have to consider the variables: (a) Concentration of active 
molecules, (b) path length, (c) temperature, (d) total pressure, (¢) 
effect of mixing of two or more active gases, and (f) effect of 
different nonactive, diluent gases. A study was recently re- 
ported by Howard, et al. [12] which considered to a greater or 
lesser degree, all of the variables given except temperature. This 
study was for meteorological purposes and is not adequate for 
most engineering applications. This work does provide an op- 
portunity to ascertain the validity of the theoretical band ab- 
sorptions. A summary of the range of variables considered by 


Table 1 Range of variables for the study by Howard, et al. | 12) 
Al] date obtained at temperature of 531°F with - 0, 
mixture (correspunding to sir) as « diluerit. 


Partial Pressure 
Range of Absorber 
(atmos. ) 


Total Pressure 
( atmos.) 


Concentration 
(atmos. ft.) 


0.033 = 26.3 0.0003 = 0.026 0.026 = 0.96 

- 51.3 0.001 = 0.013 

51.3 0.001 = 0.066 
0.001 = 0.045 


0.02 0.099 


0.0033 - 0.033 
0.0026 = 0.033 
0.0026 - 0.033 
0.004 = 0.051 
0.00 = 0.037 
0.00h = 0.037 
0.007 


0.0033- 0.98 
0.0026- 1.0 
0.0026 - 1.0 
0.004 0.97 
0.007 - 0.97 
0.013 0.97 


=~ 0.037 0.013 - 0.97 


ee’ 
“net Given 
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Howard and his co-workers is given in Table 1. The correlations 
obtained for both water vapor and carbon dioxide were of the two 
alternative forms: 


ag = cw/XP + py (35) 
az = C + Diogw + G log (P + p) (36) 


The first expression is valid for band absorptions less than an ex- 
perimentally determined value of absorption and the second is 
valid for absorptions greater than this value. If the value of w 
were increased indefinitely, the absorption given by these corre- 
lations would not reach a limiting value as is physically necessary. 
Consequently, a third expression should exist for larger values of 
the variable w. 

The form of equation (35) for the variable w can be shown to 
follow from the Elsasser model by expanding equation (29) and 
retaining the first term of the expansion. A qualitative explana- 
tion for the form of equation (36) has been given by Elder and 
Strong [5]. The dependency of both equations upon total and 
partial pressure is attributed to the overlapping of adjacent lines, 
self-broadening, and total pressure broadening. Total pressure 
broadening follows from equation (32) but, experimentally, an 
exponent was found which was smaller than predicted theoreti- 
cally. The partial pressure broadening enters the two expressions 
in a double dependency since theoretically and experimentally, 
self-broadening is stronger than total pressure broadening [12, 
15}. The correlation of both carbon dioxide and water vapor with 
equations of the same form differs from the results indicated by 
the theories of Elsasser and Goody. The experimental results 
were matched satisfactorily to the theories of Elsasser and Goody 
for specific bands but the forms of the correlations chosen by 
Howard were found to be more accurate for both gases over a 
wider range of variables and the same form was retained for both 
gases. The constants of the correlations for the various bands 
are given in Table 2 to indicate the magnitudes and relative im- 
portance of the various bands. The band limits given were ar- 
bitrarily chosen by Howard and are not necessarily indicative of 
true physical limits. Howard, et al. [12] have also considered the 


Table 2 Summary of results obtained by Howard, et al. [12] 
Carbon Dioxide 


Nominal Band 
Center Gavelength 


(Microns) 


“Weak* Band Fit "Strong* Band Fit 


Band Limite c 
(Microns ) 
12.5 = 16.19 
5-05 = 5035 
4-63 - 5.05 
4.0 - 
2.63 - 2.87 
1.92 - 2.11 
1.53 - 1.67 
1.38 - 1.50 


B 


4.88 - 8.70 
2.99 - 3-57 
2.27 = 2.99 
1.87 1.69 - 2.08 
1.38 1625 
1.1 1.06 - 1.20 


0.9% 0.87 = 0.99 0.0042 0.27 


_ The weak" band fit is of the form ag = ow/2(p + p) and is 
limiting value of absorptivity given. For \wities 
limit, the “strong band fit must be used: a, = C + D log w + G leg (P + p). 
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| 
j 
(Microns) 
1s 26 0.272 0.220 0.188 
4. Carbon Dioxide 
5.2 22 0.40 — — 
1s 
% 
5.2 n 
bed 267 0.43 0.213 
; 
2.0 0.39 “1.191 0.307 0.253 
2.7 \ 
1.6 oO. ao 
8. Water 
17.7 = 263 0.066 0.099 0.099 = 1.0 — 
thes 3B. Water Center Wavelength Band Limits 
0.060 0.30 18 0.062 O.2h2 0.17h 
3.2 0.09 = 101 ‘ 
2.7 = 167 
0.045 0.32 #9 0.063 0.232 
1.87 1.32 = 167 
ea 0.021 0.30 2 232 0.1m 
a 1.38 2.2 = 79 
0.016 0.30 23 0.369 0.307 0.232 
AS 0.00k7 0.26 2 
0.9% 6.6 -101 


problem of mixtures of carbon dioxide and water vapor which 
have overlapping spectra in certain frequency regions. The re- 
sults for mixtures were correlated for only one band (2.7 mi- 
crons): 


a,(H,0 + CO,) = ag(H,O) + (37) 


The factor € is empirically determined from experiment and is a 
function of the band absorption of the water vapor. This is 
similar to the correction given by Hottel [10] for the total absorp- 
tion of a mixture of these two gases. 

The experimental technique used by Howard, et al. [12] was 
based upon the theory of Nielsen, Thornton, and Dale [15] which 
predicted the band absorption to be independent of the resolving 
power of the spectrometer under certain limitations. The basic 
restriction of this theory is the requirement of a constant incident 
energy over the whole band. Experimentally, this was not found 
to be an important restriction if a source was chosen which did 
not vary significantly over the band, e.g., an ideal radiator at a 
high temperature. This condition corresponds to the theoretical 
restrictions of Elsasser and Goody. 

The effect of temperature has not yet been measured experi- 
mentally but, as indicated previously. the expected behavior 
may possibly be of the form: 


a, = + p)e ( (38a) 
or 


T. 
ay = C+ Dog w + Glog (P + p) + H log ( (38b) 


Engineering Applications 


The most important result obtained from the theoretical 
analyses of Elsasser and Goody and the limited experimental 
work of Howard, et al., is the nonexponential behavior of band 
absorption. Also, the form of the expression for absorption de- 
pends upon the value of the absorption and not directly upon the 
path length. For a given value of the variable w, the total ab- 
sorption would be given by: 


B=% 
f, Jv =1 Av, =i Av, 
L=¥ 
3 
a (39) 


L=1 


An analogous expression for the transmission may be written 
using the definition presented previously. The absorption given 
by equation (39) does not include the bands which have reached 
maximum absorption but this can be adjusted by adding a third 
series and defining another total band absorption similar to 
equations (38a) and (38b). Equation (39) corresponds to the ex- 
perimental total absorptions measured by Hottel [10] and the 
theoretical values computed by Penner [17] with V = 0. 

The foregoing definitions of band absorptivities and transmis- 
sivities may be applied to the problem of the net exchange be- 
tween two surfaces separated by an absorbing medium. This 
results in the definition of geometrical band or line absorptivities 
and_transmissivities: 


(dza)e, = f 2 (40a) 
Ai Ae mr 


6 
(dz)e, = f (es)e dA,dA, (40b) 
Ai 
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The expressions for the transmissivities can be written by sub- 
stituting band or line transmissivities for the absorptivities. 
Then, by the use of equations (19a) and (19b), the relationship 
between transmission and absorption is found to be continued: 


(dg = 1 (41b) 
a, + =1 (Alc) 


Similarly, the total geometrical absorption and transmission 
may be written: 


B=% 
J dv =1 
0 
B=® L=¥ 
+ J dv + a Jair} (42a) 
=1 
B= 
= 


f. Jay 

=1 Avg 
B=, L=¥ 

J dv + f (42b) 


&a+7=1 (42c) 


Substitution of the foregoing expressions into equation (6b) gives: 
= (Jn — | 


— AF: [ f, f, | (43a) 


or 


The integrated radiosity of a totally absorbing surface, i.e., an 
ideal radiator, is ¢7'* and equations (43a) and (43b) reduce to the 
conventional exchange equations for this condition. If there is no 
absorbing gas, 


72 = = 1. 


Discussion 


The problem of the effect of gaseous absorption upon the net 
exchange between two surfaces has been examined from present 
theory and from the limited available experimental information. 
The assumption of exponential absorption (Beer’s law) for the at- 
tenuation of a band has been shown to be incorrect for a single 
line or for an assembly of such lines as occurs in a vibration band. 
The best available experimental data indicates the absorption of 
a band follows a square-root function of the optical “thickness’’ 
(w) for a weak absorption and a logarithmic function for a strong 
absorption. The logarithmic dependency for large values of w 
may be shown with existing data for gaseous absorption [10]. 
For example, a plot of the logarithm of total absorption versus w 
for CO, at a constant temperature yields a straight line for w 
greater than 0.5 atmosphere-foot. This corresponds to equation 
(39) with ®, and L equal to zero. 

The various bands of a specific gas for a given temperature and 
pressure will not all behave according to only one of these two 
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types of functions. Consequently, the behavior of the gas over 
the whole energy spectrum will be a mixture of these two func- 
tions, superimposed upon unity transmission in those regions of 
the spectrum where the gas is not active. The total effect of the 
gas must be considered therefore as the sum of the effects of the 
individual bands. The complexity of this approach is readily ap- 
parent since the number of computations necessary are multiplied 
by the number of vibration bands occurring. 

An advantage is gained, however, by considering the band or 
line absorptivities (ag or a,) or the geometrical band or line ab- 
sorptivities (4, or d,). Theory which has been substantiated by 
experiment shows that these quantities are not strongly dependent 
upon the source but are primarily functions of the temperature 
of the absorbing medium. The source dependency is introduced 
when the comparable total factors are considered, i.e., a or &. 
This is important in the coraputation techniques available for 
analysis. 

The calculations of gaseous absorption are complicated when 
complex reflecting geometrical configurations are conside ed. The 
network technique by Oppenheim [16] is well suite. to this. 
In the last part, a procedure of solution will be described which is 
amenable to either band or digital computer calculation. 

The use of band absorptions for gases and computations using 
this information also requires the use of spectral emittances and 
reflectances for the surfaces involved. Information of this type 
is not available in quantity but a recent paper [1] has described 
a method for performing such measurements. The importance 
of considering spectral emittance rather than total for calculating 
the net exchange between surfaces with no intervening absorbing 
medium has been demonstrated by Goodman [7]. The use of 
band absorptions for calculations, rather than total absorptions, 
is only an additional step beyond the work of Goodman. The 
data for band absorptions for carbon monoxide, carbon dioxide, 
and water as a function of temperature and total pressure are not 
now available but an experimental program has been initiated 
to partially eliminate this deficiency. 

The importance of the nonexponential absorption upon the 
present methods of analysis is difficult to assess. An example 
will be given in a subsequent paper of the errors involved for a 
simple geometry. The historical work by Hottel and his co- 
workers [10] remains the best and only engineering approach 
presently available. Penner [17] computed the total emissivity 
of carbon monoxide by considering the individual contributions 
of each line in the several vibration bands and has confirmed 
theoretically the experimental measurements of Hottel. How- 
ever, the assumption of exponential attenuation for such analyses 
as the volume emission of a gas [13] is questionable. 


Nomenclature for Part 1 


(The units of the various quantities should be in a consistent 
system of, units.) 


area 

constant, equation (36) 

constant, equation (36) 

emissive power of an ideal radiator 
geometrical form factor 

constant, equation (36) 

constant, equation (38b) 

intensity of radiation 

radiosity 

total pressure 
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universal gas constant 

integrated intensity, equation (21) 

absolute temperature 

energy (AW = energy change) 

total absorptivity of a single line or band, equations (13) 
and (14) 

constant, equation (35) 

spacing between adjacent lines, equation (26) 

constant, equation (35) 

Planck’s constant, 6.623  10-* erg-sec 

absorption coefficient 

mass per unit area of the radiating components in the 
medium lying in the path of the radiation 

partial pressure 

heat rate 

path length traversed by the radiation beam 

total transmissivity of a single line or band, equations 
(13) and (14) ' 

optical thickness, equation (33) 

constant, equation (38a) 

geometrical absorptivity, equation (40) 

geometrical transmissivity, equation (41) 

radiation passing from one surface to another 

number of vibration-rotation bands having w'" be- 
havior 

number of vibration-rotation bands having log w behavior 

number of lines external to a vibration-rotation band 

total absorptivity 

parameter, equation (23) 

parameter, equation (28) 

line half width 

correction factor, equation (37) 

parameter, equation (28) 

number of lines, equation (26) 

angle 

frequency 

3.14159... 

density of the absorbing (or emitting) medium in the ra- 
diation path 

Stefan-Boltzmann constant = 0.1714 XK 10-* Btu/hr ft? 
(deg R)* 

total transmissivity 

total number of vibration-rotation bands 

frequency width, i.e., the frequency region within which 
absorption (or transmission) is complete 

total geometrical absorption, equation (42a) 

total geometrical transmission, equation (42b) 


e 


2 aa 


F 


vibration-rotation band 
line 

incident 

ideal (black body) radiation 
electronic 

center or reference 
rotational 

vibrational 

designate surface 

medium 

surface 1 plus 2 

from surface 1 to 2 

single frequency (monochromatic) 


ODD 
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Pat General interchange Equations 


Introduction 


HIS IS THE second paper of three dealing with the 
solution of the radiant exchange within an enclosure containing 
an absorbing medium. 

The starting point of the analysis to be given is the network 
analogy for radiative heat transfer given by Oppenheim [16]. 
The network analogy utilizes a resistive network to represent the 
interplay between the geometrical and physical characteristics of 
the radiation exchange system. An illustration of the conduct- 
ances coupling two surfaces of an enclosure containing an ab- 
sorbing medium is given in Fig. 2(a) for monochromatic radia- 
tion. The two surface node potentials (Z,, and E,.) are the 
monochromatic emissive powers of the two surfaces, respectively. 
The radiosities of the surfaces are given J,,, and J,,. and are the 
potentials of the inner nodes of the network. The quantity 
W,..m is the monochromatic emissive power of the medium lying 
between the nodes 1 and 2 in accordance with the description 
given in Part 1 for the radiant exchange between two surfaces. 

The solution of the simple system given in Fig. 2(a@) may be ob- 
tained quickly by basic electric network concepts. The direct 
solution of a system of only four or five nodes can become com- 
plicated and tedious and recourse is generally made to matrix 
algebra. A reiterative (relaxation) procedure for approaching a 
solution to a multinode network will be given in the final part. 
Prior to this, however, a generalized network for an enclosure will 
be used to examine the physical bases for solution and the ap- 
proximations necessary to obtain the solution. 


Monochromatic Energy Analysis 


Consider a radiation system of M surfaces with a part or all of 
the M surfaces forming an enclosure; the balance of the M sur- 
faces lie within this enclosure. An absorbing and/or emitting 
medium is contained within the enclosure. The number of sur- 
faces is not specified, the aggregate surfaces of the enclosure are 
considered to have been subdivided into M surfaces such that 
each of the M surfaces is isothermal; each is assumed to follow 
Lambert’s law of diffuse radiation and to be in equilibrium with 

* The term node will hereafter refer to an inner node such as J»,1; 


surface node will designate a network point with a potential such as 
Ey; and gas node will be used for such points asWy,,12. 


e 


Gyi2AiFi2* Ae 


Aifie® Ter 


Fig. 2(¢) Mutual conductances coupling two surface nodes in an en- 
closure containing an absorbing medium 
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its surrounds (i.e., no heat storage). The radiation characteristics 
of the medium are assumed to be isotropic and the medium to be 
in equilibrium, The nth node of such a system is shown in Fig. 
2(b). The potentials of the gas nodes have been labeled by the 
subscript referring to the surface which determines this node in 
conjunction with surface n (i.e., W,,, = W,,). 

The sum of the currents flowing into node n must be equal to 
zero by Kirchhoff’s law of electrical networks or since the currents 
are equivalent to the heat fluxes, by the conservation of energy. 
If this is done, the following expression is found (n < m): 


M 
Cyn 
Ag (Brn Ina) = Gan ~ 


m=1 
M 
+ — (44) 


Solving for the radiosity of the nth node (J,,,) gives: 


M 


m=1 


M 
= A, + >, 


Tym m=1 


M 
+ (45) 


m=1 


This equation may be simplified by noting: 


M 
= Ae (460) 
m=1 
Equation (45) then reduces to: 
M 
Cyn Cyn 
Jordy + = A, + Pan 


M M 
— Pam +) (47) 


m=1 m=1 


Fig. 2(b) The ath node of a radiation network for an enclosure con- 
taining an absorbing medium 
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An alternate procedure of analysis may be developed from equa- 
tion (47) and is given in the Appendix. The method to be given 
in the main text retains the definitions given in the first part. 

The total radiant energy is given by integration of equation 
(47) over all frequencies: 


M 
m=1 0 


W,, nay, amA (48) 


m=1 


This expression represents an exact formulation of the energy in- 
terchange at node n. The solution depends upon the evaluation 
of the various integrals and the simultaneous solutions of the M 
equations of which equation (48) is the nth. Assuming the radi- 
ant properties to be independent of frequency will give the 

gray” radiation case. The approximation which is between the 
extremes of the exact and the gray cases, is the band energy ap- 
proximation based upon the definitions given in the preceding 
paper Prior to presenting this, the total energy balance within 
the system should be obtained. Consider equation (44) summed 
over all nodes: 


= 


M 
n=1 m=1 
M 
+>: > (Join — (49) 

n=l 
This expression assumes there is no energy storage within the 


system since equilibrium was assumed for the surfaces and the 
medium. 


The definition of the monochromatic heat flow from the surface 
to the node used to derive equation (44) was: 


where the heat flowing inéo the node is taken as positive. Conse- 
quently, the right-hand side of equation (49) is the total mono- 
chromatic heat flowing out of all the nodes: 


M 
= >> (51) 


n=1 


The first term on the right-hand side of equation (49) is equal to 
zero. The double sum will cause a reversal in the sign of the 
quantity within the parentheses and this term reduces to zero 
when combined with the identity valid for monochromatic energy, 
ie., 


A simple expansion of this term for M = 2 will demonstrate this 
readily: 


inp 


n=1 m=1 
+ — + — = 0 
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The second term on the right-hand side of equation (49) re- 
lates the radiosities of the various nodes and the emissive power 
of the enclosed medium. Considering first the radiosity portion 
of this term, the double sum is the sum over all nodes of the total 
geometrical absorption of the medium, i.e., the total monochro- 
matic energy absorbed by the medium. Similarly, the total 
monochromatic emitted energy from the medium is represented 
by the remaining portion of this sum. The difference between 
these two factors is the total monochromatic energy absorbed by 
the medium: 


M MOM 
n=1 n=1 m=1 


No other sources or sinks were postulated and equation (53) 
merely states this in terms of the radiosities, heat fluxes, and 
radiant characteristics of the system. The right-hand side of 
“eine (53) may be integrated over all frequencies: 


n=1 0 


n=1 m=1 
(54a) 


or 


ff,” = 0 (54b) 


Band Energy Approximation 


As noted previously, equation (48) is an exact formulation of 
the radiation interchange problem. The solution of the M simul- 
taneous equations is a formidable one and requires a specification 
of the frequency characteristics of each factor in the equations 
plus the simultaneous solution of the equations. As an inter- 
mediate approximation between the exact solution and the as- 
sumption of gray radiation, the procedures of the first part may be 
used. This method has been termed the band energy approxima- 
tion. 

In order to apply these procedures in a systematic manner, 
equation (48) is rewritten in its complete form: 


0 
n m 
m=1 An 
af cos 6, cos 
m=1 


The integration of equation (55) over all frequencies is written 
as the sum of integrals over frequency intervals (bands). These 
intervals are the natural subdivisions determined by the absorp- 
tion (or emission) characteristics of the intervening medium. 
The media of usual interest in engineering are gases and the 
procedures outlined in the first paper are applicable. The entire 
spectrum is divided, therefore, into Q intervals of which «, inter- 
vals have no absorption and w: intervals have absorption.’ 
Thus the @ intervals are the line absorptions and vibration-rota- 
tion frequency bands. The energy within each division must be 
conserved and equation (55) may be written in general form for 
one interval as 


7In application, certain « intervals may be further subdivided if 
the spectral characteristics of the surfaces vary rapidly within one or 
more of the w intervals. 
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(56) 


Q=w, +a: = Zw 
a=0 for w =o (57) 


The last three terms of equation (57) are simplified by the use of 
the following expressions from the first part: 


Av. ay 


where the factor a is a function of the path length R and 


= £00 8, dA.dA, (59) 
Aa 


This gives for these three terms: 


M 
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m=1 
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M 
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M 


m=1 
These approximations are based upon the following assumptions: 


1 The radiosity or emissive power varies slowly within the 
interval. 

2 The radiosities of the various nodes follow Lambert’s 
“eosine’’ law. The second assumption is inherent in the network 
method of analysis but is no more restrictive than other pres- 
ent methods of analysis. 

The remaining two terms require the definition of two new 
quantities involving the emittance and reflectance properties of 
the surfaces. Using the same assumptions as for the gaseous ab- 
sorptivities, these quantities are: 


( : f Jna( dy 
Av. Judy 
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The remaining terms of equation (56) become, with the aid of 
these definitions: 


Tun (P)w.n ar, 
A, E, dv Yo.nA, E, (66) 
av, Tym ary 


Substitution of the appropriate foregoing relations in equation 
(56) yields: 


A 
Jy = E, ,dv 
(P)w,n ar, Av, 


M M 
w 


m=1 


M 
ar, 


There will be w, equations having the form of equation (67) and 
@, equations without the last two terms, which must be solved. 
The Q equations are connected by the temperatures and/or heat 
fluxes of the various nodes, e.g., the monochromatic emissive 
power of node n in one frequency interval is related by the tem- 
perature of that node to the monochromatic emissive power of 
the other frequency intervals. The heat flow of a node for a spe- 
cific band is given by: 


Jen = fs. Qv.ndv = Yond, | any E, ,dv - ary Jr ait | 
(68) 


The total band energy may be obtained by starting with equa- 
tion (49) and proceeding with the same approximations as used 
to obtain equation (55). This gives: 


f Jy 
1 an, 
- Wa (69) 
ary 


Certain comments are required regarding the approximations 
given by equations (63) and (64). When equation (64) is applied 
to determine the quantity y.,., the monochromatic emissive 
power of the node may be used to determine this quantity more 
accurately than the definition given. The monochromatic emis- 
sive power is not always a known quantity for every surface. 
The definition given by equation (64) will permit general values 
of Y.. to be determined from spectral emittance or reflectance data 
and these will be temperature dependent only as the temperature 
affects the spectral emittance and reflectance, i.e., independent 
of the temperature dependency of Planck’s equation. This same 
argument holds for equation (63) and is more important to the 
analysis since the radiosity is generally unknown. The approxi- 
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mation given by equation (58) has already been discussed in the 
first paper. In each of these approximations [equations (63) and 
(64)], the accuracy is greater the smaller the frequency interval 
(Ay,,). Fortunately, the variation of spectral reflectance and 
emittance is usually small within the natural frequency intervals 
of CO, and H,O vapor (Part 1 and [1]). This physical behavior 
increases the accuracy of the approximation. Note should be 
made of the fact that neither equation (63) nor (64) correspond 
to physically observable quantities but arise from the method of 
analysis. 


Gray Radiation Approximation 


The gray radiation approximation is the one most frequently 
used in analyses of radiant systems. The approximation is based 
upon the assumption of frequency independent radiation charac- 
teristics of the surfaces and the gases involved. This assumption 
is a result of the inadequate knowledge of the spectral behavior 
of these characteristics and is a necessity born from ignorance. 

Equation (48) may be quickly reduced with the gray radiation 
assumption by simply dropping the subscript denoting frequency. 
Alternatively, an evaluation of this assumption may be obtained 
by simplifying equation (47) in a systematic manner to provide 
expressions which can be used to desermine the errors involved. 
Multiplying equation (47) by the monochromatic reflectance and 
integrating over all frequencies gives: 


a M 
0 0 m=1 0 


M 
m=1/0 


W,, ml yp, (70) 


The following expressions for the individual terms in equation 
(70) are now necessary to make the equation useful for computa- 
tion: 


Ady = 


A,e,oT,! A, Ey ndv 
- 


Equation (71) is actually a definition of the integrated radiosity 
and does not introduce an error in itself. Equation (72) is the 
definition of the total emittance of a surface: 


Ey, (76) 


Equation (73), however, does introduce an approximation which 
cannot be evaluated. The spectral energy distribution of the 
radiosities of the various nodes is not known nor can it be de- 
termined by this analysis since the analysis is predicated on a total 
integrated energy. Instead, a further approximation is introduced 
by writing: 


= = 1 — 
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(77) 


Equation (77) is the first term requiring the use of the gray radia- 
tion assumption. 

A similar series of approximations are necessary in equations 
(74) and (75) with the added complication of the gaseous absorp- 
tivities. Each factor in the two equations is considered to be in- 
dependent of the others. Thus the geometrical absorption factor 
(Gam’) is based upon the concept of mean beam lengths which in 
turn is derived from a specified transmission law. The geometrical 
shape factor (A,F,,,) is obtained from calculations made without 
absorption of the intervening medium and corrected for absorp- 
tion by the factor 4,,,’.. The reflectance of the surface is deter- 
mined separately in a manner analogous to equation (77). 


These factors may be written as: 
0 


Ey 
E, »dy 
f f (78) 
a 
An identical equation may be obtained for the gas by substituting 
W,.» for Jy, in the last bracketed term. The quantity 4,,,’A,F am 
is generally evaluated by the procedures given by Hottel [10]. 
In brief, the procedure uses a “mean beam length’’ to determine 
the effective path of the radiation and this in turn is used to de- 
termine the emission and/or absorption of the gas from data col- 
lected by “parallel beam’’ experiments. The geometrical shape 
factor (A,F,,,,) for the two surfaces with no intervening medium 
is then multiplied by this absorption factor to obtain the product 
Ginm'AnF'nm- Recently, Hottel and Cohen [11] have made an 
important step toward improving the values of the geometrical 
factor. The absorption factor was introduced within the shape 
factor integral as a part of the integrand and values of the integral 
determined numerically. 
The gray radiation approximation may be written using the 
terminology given in equations (71) through (75) but with the 
additional approximations outlined in the discussion: 


M M 


m=1 m=1 
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The total energy balance for this approximation is: 
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n=i m=1 
Discussion 


Three methods of analysis of the radiant interchange within an 
enclosure have been presented in order of increasing approxima- 
tion and decreasing complexity. The first is an exact formulation 
of the problem [equation (47) or (48)] and is only restricted by 
the assumption of Lambert’s law (diffuse radiation). Unfor- 
tunately, application of this method to a physical situation is a 
practical impossibility. Solution of the M expressions given by 
equation (47) would require a solution at each frequency with 
the number of solutions approaching infinity. Solution by means 
of equation (48) is not feasible since the integrations required in- 
clude factors which are not sufficiently described physically as a 
function of frequency and temperature (e.g., the gaseous absorp- 
tion factor a,) and/or include factors which are to be determined 
(e.g., the radiosity J,,m). These expressions do provide the initial 
step in making the approximations necessary for a practical solu- 
tion and, consequently, a measure of the errors involved. 


(80) 
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The second procedure is the band energy approximation as given 
by equation (67). The various approximations used in this 
analysis are based upon the fundamental assumption used in the 
theoretical analysis of gaseous absorption and verified experi- 
mentally, i.e., the absorption may be adequately described by 
assuming the incident energy to be constant over the band. 
Equations (63) and (64) are extensions of this assumption to the 
radiant characteristics of the surfaces. Since the frequency de- 
pendency of these factors is generally less than that of a gas, this 
extension should be at least as valid. The effect of this assump- 
tion is to treat the radiant system as gray over small frequency 
intervals and to approach the actual system [equation (48)] in a 
series of steps rather than in one large one [equation (79)]. 

The band energy approximation may be improved by using the 
monochromatic emissive power of the surface in combination with 
the spectral characteristics of the surface as indicated in the first 
expression of equation (64). The monochromatic emissive power 
of the surface may also be used as an approximation in equation 
(63). If this is done, however, the values of y,, and p,, become 
doubly temperature dependent as noted previously. 

The application of the band energy method is not possible 
within present knowledge. Neither the magnitude nor the func- 
tional behavior of the absorption within the various frequency 
intervals has been studied experimentally over the range of varia- 
bles required for an engineering analysis. Methods of measuring 
the spectral characteristics of the surfaces are feasible and a 
limited amount of information is available [1]. 

The final technique of analysis described is the one most com- 
monly used in calculation of radiant exchange, the gray radiation 
approximation. This procedure assumes that all of the surfaces 
and the gas have absorption characteristics independent of fre- 
quency. The number of equations which must be solved simul- 
taneously is reduced to M and only one set is involved. 

The assumption of frequency independent characteristics is 
obviously wrong whenever a gaseous medium is included. The 
errors involved are significant for the case without the gas [3] 
and the supposition must be made that a gas would add to these 
errors. The inaccuracies of the methods presently available for 
evaluating the geometrical absorption factors are not known. 
Until an analysis is made based upon the band energy approxima- 
tion, the errors of the geometrical absorption factors cannot be 
estimated. 

The particular case of refractory surfaces (surface with g, = 0) 
commonly used in the gray radiation analysis must be imposed as 
a special condition for the band energy approximation. For a 
specific band, the refractory surface does not exist except as a 
very special case, since by definition a refractory surface is: 


Q 
Jen = 0 


Equation (79) for the case of a refractory surface becomes: 


(81) 


M 


M 
pe 
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m=1 


(82) 


. The total energy balance provides a means for separately 
checking a solution of either the gray or band energy approxima- 
tion. In addition, a primary condition given for many physical 
problems is the total energy provided by combustion gases. 
For the band energy approximation, the total energy emitted by 
the gases is obtained from equation (69): 
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The analogous expression for the gray radiation approximation is 
found directly from equation (81): 


M M 
H = -Q= — Su) 


n=1 m=1 


(84) 


Solution of the M simultaneous equations may be secured by 
established methods of matrix algebra. These methods become 
difficult and time consuming for even digital computers as the 
number of nodes are increased. An alternative method of solu- 
tion will be described which is analogous to a relaxation of the 
system of equations. This method approaches the solution of the 
equations by a series of approximations and permits the selection 
of a desired accuracy of solution. This technique and numerical 
examples of its application will be described in the final part. 


APPENDIX 

The band energy approximation may be obtained in a form 
slightly different from that given in the text. This alternative 
derivation requires approximations of slightly different form but 
of no greater degree of approximation than those given in the 
text. The choice of the two forms is subject to the particular 
conditions of a problem since conditions are conceivable whereby 
the method given here would be preferable to that of the text. 

The initial starting point is equation (47) of the text which has 


been altered by multiplying by the monochromatic reflectivity of 
node n: 


m=1 


> Ty, gl am + Ty. mW (85) 


Proceeding in the same manner as before, the band energy equa- 
tion for node n and frequency interval Ay, may be obtained: 


f 
ar, 
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cos 6, cos 0 
m=1 


The approximations given by equations (58) through (62) are, 
for this system: 


(87) 
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The approximations involving the spectral characteristics of the 
surface [equations (63) through (66)}] are not required as before. 
These have been replaced by a fundamental approximation in- 
volving the product of the spectral reflectance of node n and the 
radiosity of surface m or gas node m. This also replaces the as- 
sumption involving the constancy of the incident energy over the 
gaseous absorption band. That is: 


Av. 


Js. 


Analogous relations hold for the other approximations. 

The inherent difficulty of this derivation is the lack of know!l- 
edge regarding the spectral reflectance of the m surfaces in a par- 
ticular problem. The frequency variability of the absorption 
characteristics of the gas are sufficiently established to approxi- 
mate the absorption within a band under the assumption of inci- 
dent energy varying as a black body. In the text this has been 
extended to the radiosity of the various nodes. The validity of 
this extrapolation is subject to serious question since the absorp- 
tion of a gas is not uniform within a band and a distortion of the 
frequency-energy distribution is automatically introduced in the 
radiosities. The addition of the spectral reflectance as a factor to 
this compounds the unknowns within the approximation. The 
spectral variation of the reflectance could increase the distortion 
of the radiosity frequency variation just as well as reduce it. 
As noted in the texi, the spectral reflectance for many surfaces 
is uniform within a limited frequency band so as not to affect the 
approximation. This cannot be taken as a generality for all sur- 
faces or for all frequency intervals to provide a firm basis for in- 
cluding the reflectance as a factor in the approximation. When 
an absorbing medium is not involved, this procedure would not 
incur this difficulty. 
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The reduction of equation (86) to a form analogous to equation 
(67) gives: 
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This expression must be solved for each node and for each Ay,, 
as previously. The gray radiation approximation follows from 
this relation by assuming the spectral characteristics of the surface 
to be independent of frequency. 


Nomenclature for Part 2 


A = area of surface; A, = area of surface n 

E = emissive power; E,, = monochromatic emissive power 
of surface n 

F = geometrical shape factor; F,,,, is defined by the following 
expression : 


H = total energy emitted by the enclosed medium; see equa- 
“tion (83) or (84) 
J = radiosity of the surface; J,,, = monochromatic radiosity 
of surface n 
M = total number of surfaces composing a specific enclosure 
Q = net energy flowing from all surface nodes, equation (540); 
Q, = monochromatic energy, equation (51); Q, = 
energy within the band, equation (69) 
R = length of the line joining the two differential areas of sur- 
faces A, and A,, 

T = absolute temperature; 7, = absolute temperature of sur- 
face n, T,,, = absolute temperature of the medium 
(gas) involved in the exchange of energy between sur- 
faces n and m 

W = emissive power cf the medium; W,,,, = monochromatic 
emissive power of the medium involved in the ex- 
change of energy between surfaces n and m 

a = band absorption of the absorbing medium, equation (58); 
Gn. = geometrical band absorption, equation (59); 
4,,,’ = gray radiation approximation of the absorption 
of intervening medium, equation (74) or (75) 

e = emittance of the surface; r,,, = monochromatic emit- 
tance of surface n; e, = total emittance of surface n, 
equation (26) 

m = number of surfaces composing a general enclosure; for a 
specific enclosure m = 1, 2,3,....M 

n = a specific surface of the enclosure 

q = energy flow; 9,,. = monochromatic energy from the sur- 
face node n to the inner node n of the radiation network, 
equation (50); gu. = band energy from a surface n, 
equation (68) 

r = reflectance; r,., = monochromatic reflectance of surface 
= gray radiation approximation reflectance of 
surface n which is a function of the radiosity of surface 
m, equation (74) and (75); r, = ,,7, form = n, equa- 

tion (27) 
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a = absorptivity of the intervening medium; a, = mono- 
chromatic absorptivity; &,.. = geometrical mono- 
chromatic absorptivity of the medium between surfaces 
n and m, i.e., 


1 cos 6, cos 6,, 
AF am f. awk? fA, 


= band energy approximation to the ratio of emittance and 
reflectance; y. = for a specific band, equation (64) 

difference; Ay, = frequency interval for a specific band 

angle between the normal and the line R; 6, = value at 
the surface n 

frequency; Av, = frequency interval for w 


Ayam = 
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| * FIRST PART in this series reviewed the principles 
of gaseous absorption and emission of thermal radiation. The 
second developed the analysis of the radiant interchange within 
an enclosure using the concepts of the network analogy. This 
final part will be used to demonstrate the application of the pro- 
cedures discussed in the first two and to illustrate the errors 
involved in the band energy and gray radiation approximations. 
The solution of a simple cubical enclosure will be used to il- 
lustrate the relaxation or reiterative method of solving the multi- 
node network. The inadequacy of present knowledge of the radia- 
tion characteristics of surfaces and gaseous media plus the desire 
to demonstrate the procedure in a direct manner dictated the use 
of a simple system for the example. The last part of the paper 
presents a comparison of the band energy and gray radiation ap- 
proximations for a gaseous medium between two parallel infinite 
surfaces. These examples are an extension of the analysis given 
by Goodman [7] and provide magnitudes of the errors involved 
in the basic methods of analysis of radiation exchange for specific 
systems. These errors are in addition to those introduced by the 
uncertainties in the radiation characteristics of the surfaces and 
media within the radiant system. 


ae (Reiterative) Solution of a Radiation Network 
(14) 


The example selected to demonstrate the method of solution 
is a cubical enclosure (box) with a heated top, a cooled floor, and 
refractory walls. The assumption of gray radiation and no ab- 
sorbing medium will be made in order to reduce the number of 
computations. The top is assumed to be at a temperature of 
2000 F and to have an emittance of 0.90; the bottom is main- 
tained at 1000 F and has an emittance of 0.60. The temperature 
distribution of the side walls will be approximated by dividing 
these walls into four horizontal bands of equal area. The geo- 
metrical system is illustrated in Fig. 3 and the equivalent radia- 
tion network is given in Fig. 4. The temperature distribution 
over each of these six areas (top, bottom, and four bands) is as- 
sumed to be uniform as required by the network method. 

The six equations which must be solved are obtained for this 
problem from the expressions given in the Appendix of Part 2, 
rearranged to give the radiosity of the nth surface in terms of the 
radiosities of all other surfaces. The radiosity of a gray, heat- 
transfer surface of uniform temperature and irradiation, is given 
in terms of its emissive power and the radiosities of all other sur- 
faces within the enclosure by: 
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Part 3 A Method for Solving Multinode Networks 
and a Comparison of the Band Energy 
and Gray Radiation Approximations 


3.14159... 

p = the quantity analogous to reflectance required by the 
band energy approximation; (p),., defined by equa- 
tion (63) for the specific band n and surface m 

7 = transmissivity of the intervening medium; t = mono- 
chromatic transmissivity; rT, + a, = 1 

Q = total number of frequency intervals (Av) used in sub- 
dividing the energy spectrum for the band energy ap- 
proximation; 22 = w, + w: 

w = specific frequency intervals (Av) used to subdivide the 
energy spectrum; w, = intervals with no absorption 

by the intervening medium; w, = intervals with ab- 

sorption by the intervening medium 
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Fig.3 Cube with refractory walls used as an example for the reiterative _ 
method 
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Fig. 4 Radiation network for system shown in Fig. 3 
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The heat flux from this surface is expressed as: 


M 
en 
(q/A)n = (‘:) (E, =(1 — 7 J nm 
(n¥m) 


(95) 


The comparable equations for the four refractory surfaces 
must be modified by the definition of a refractory surface, i.c., 
(q/A), = Oor E, = J,. Using these conditions gives: 


M 


m=1 
(n¥m) 


In an enclosure of M surfaces, there will be a total of M equa- 
tions of the form of either equation (94) or (96) The particular 
form will be determined by the conditions of the surface, i.e., 
whether a heat-transfer surface or a refractory one. The proce- 
dure of solution is initiated by assuming radiosities for all nodes 
and then, starting with node 1, solving for the radiosity of node | 
using the assumed radiosities for nodes 2 through M and the per- 
tinent equation for this node. This new value of the radiosity of 
node 1 is used along with the assumed radiosities for nodes 3 
through M to find a new value for the radiosity of node 2. The 
procedure is continued until new values for the radiosities of al! 
nodes have been found. A second trial is made by repeating the 
computation with these new radiosities and the trials continued 
until the desired accuracy is achieved. This accuracy may be 
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= 3,650 
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Table 3 Solution of radiation network shown in Fig. 3 


established in a number of ways such as requiring for each node 
the difference between the radiosities of successive trials to be 
less than a fixed percentage of the last trial radiosity. Accuracy 
may also be determined using the heat flux computed from the 
radiosity and emissive power of the node. The convergence of 
the procedure described is rapid and generally only a few trials 
are necessary for a moderate accuracy. Judicious guessing from 
experience of the initial radiosities will reduce the number of 
trials required. 

This procedure is rapidly applied to the sample problem. The 
potentials (radiosities) of all six internal nodes are unknown as 
are the currents in the internal conductances. The currents to the 
external nodes of surfaces 1 and 6 are equal and opposite since 
the remaining surfaces are assumed to be refractory and no other 
heat sources or sinks exist in the system. The external potentials 
of nodes 1 and 6 are fixed while the potentials of the remaining 
nodes are “floating.’’ The expressions for the radiosity of nodes 
1 and 6 for the particular conditions of the example are obtained 
from equation (94): 


6 
J, = 56,664 + 0.10 (97) 
m=2 
5 
Je = 3,100 + 0.60 JFem (98) 
m=1 


Similar expressions for the remaining nodes are obtained from 
equation (96). The expression for node 2 is, for example: 
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1-2 0. 360 50, 00 19,000 46,000 17,500 1-2 51,66 :19,625 51,865 19, 709 
0.205 43.000 8,200 37,300 7,650 1-3 47,684 9,775 47,926 9,826 
0.135 30,000 4,050 30,800 4,160 1-4 ub, 218 5,969 ub 6,001 
2 1-5 40, 353 3,226 0,560 3,265 
t 
28, 390 36,418 3, 909 
Jy 37,300 ul,100 7,926 48 
0.135 6,150 8,120 60, 300 8,150 62,252 8,256 061,172 258 
u-2 0.070 46,000 9220 7,400 4-2 51,865 3,631 51,978 3,638 1 
33,77 33,6 ‘| 
‘ 


The solution is shown in Table 3 using the following initial as- 
sumptions for the radiosities: J, = 60,000; J. = 50,000; J; = 
40,000; J, = 30,000; J; = 20,000; J, = 10,000. Only trials 1, 
2, 7, and 8 are shown in Table 3 for brevity. The difference be- 
tween the values of trial 7 and 8 are seen to be less than 0.25 per 
cent indicating the values to be very close to the correct ones. 
An additional check shows the heat fluxes (currents) of surfaces 1 
and 6 are equal within 1.3 per cent and for this particular ex- 
ample, the true value would be very near the arithmetic average 
of these two values. 

The particular example chosen is not difficult to solve by 
established methods [10] and provides an opportunity for com- 
parison. The equivalent circuit between nodes 1 and 6 may be 
drawn as shown in Fig. 5 where A,F,-+ is the value of the conduc- 
tance between the two nodes. The value of this conductance 
given by Hottel [10] is: 


| = 0.536 
Using the results obtained on the eighth trial (Table 3), we have: 


(q/Ah 15,980 
= 7 7 61,172 — 31,508 
The comparison between these two values is satisfactory. 

The radiosities of the four areas composing the walls of the en- 
closure may be used to determine the temperatures of these areas. 
Since these are refractory surfaces, the radiosities are equal to the 
emissive powers of a black body: 


J = 
Calculation of these temperatures gives the following values: 
= 1884 F = 1797 F 
= 1838 F ts = 1744 F 


These temperatures may be assumed to be the average values of 
the particular areas, or that temperature at the mid-height of the 
area. The equivalent radiant temperature of a surface is equal 
to the temperature of a black body which would have a total emis- 
sive power equal to the radiosity of the surface. This radiant 
temperature is the temperature which would be registered by a 
radiation pyrometer sighting at the surface. Errors using an 
optical pyrometer are liable to be even larger than when using a 
total radiation pyrometer, and the effect of nongrayness may 
either increase or decrease the error, depending upon the spectral 
characteristics of the surfaces and media in the system. For the 
problem just worked, the equivalent radiant temperatures found 
for the top and bottom surfaces are 1980 F and 1609 F, respec- 
tively. When compared with the true temperatures of 2000 F 
and 1000 F, the influence of the surface characteristics and of the 
surrounds upon radiation pyrometry is evident. 

In the preceding example, either the emissive power (tempera- 
ture) of a surface or the heat flow into the surface was known 
(zero for the refractory surfaces). This type of problem is a 
special case of a more general problem which may be considered. 
The unknowns in any problem must equal the number of nodes in 
order to insure a single solution. The number of known or 


Fig. 5 Equivalent circuit for radiation network between surfaces 1 and 6 
of the cube 
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specified quantities may be assigned to the nodes in any manner 
compatible with this restriction. Thus a node may have both 
the heat flux and the emissive power specified, either the heat 
flux or the emissive power may be given, or neither of these quan- 
tities may be given. 

The application of the reiterative method to high speed digital 
computers is readily apparent and avoids the difficulty of matrix 
inversion required by other methods of solution of the radiation 
network. The limitation of the use of reiteration for hand cal- 
culation is the patience of the computer. 

The use of the band energy approximation for nongray systems 
results in a large increase in the complexity of the problem. If 
the temperstures of all surfaces and the gas are known, the pro- 
cedure is straightforward. Each band is solved separately, the 
heat-transfer rates are obtained for each surface by summing the 
bands for that specific surface. If only a few surfaces are in- 
volved, the radiation network can be solved analytically; but if 
many surfaces are involved, it is better to handle each band by 
the relaxation method as outlined for the gray system. If many 
bands and nodes are involved, the use of a computer is advisable. 

If any of the surfaces in the system have specified heat rates 
(e.g., refractory surfaces with zero heat flux) and unknown tem- 
perature, or both unknown temperature and heat rate, the band 
energy approximation is much more difficult to apply. Particular 
attention should be given the refractory surfaces since the sum of 
the heat fluxes for all bands may be zero, but the value for any 
one band is not necessarily zero. The solution of a problem may 
be facilitated by first solving on the basis of gray conditions and 
using the results obtained to start the reiterative solution. This 
procedure should provide better initial estimates than experi- 
enced guesses for the unknown radiosities to be used to start the 
solution. 

Solving an enclosure using the band energy approximation is 
presently restricted by the lack of adequate knowledge of the 
radiant characteristics of the surfaces and gaseous media. In 
addition, the complexity of such a solution would have obscured 
the procedure used in solving the network. Consequently, the 
reiterative method of solution has been separated from a discus- 
sion of the band energy and gray radiation approximations. This 
latter comparison has also been restricted to simple conditions 
in order to minimize the complexity. 


Comparison of the Band Energy and Gray Radiation 
Approximations 


The first case to be considered is two infinite parallel plates 
separated by a nonabsorbing and nonemitting medium. The ra- 
diant exchange will be computed using the exact solution, the band 
energy approximation and the gray radiation approximations. 
This is similar to the analysis made by Goodman [7] for the exact 
solution and the gray radiation assumption. The exact solution 
and the band energy approximation differ only in the number 
and size of the frequency intervals used for a particular problem. 
When an absorbing medium is present, these intervals are deter- 
mined by the characteristics of the medium and for this example 
the band limits of carbon dioxide have been used according to the 
work of Howard, et al. [12]. In addition, the “window” region 
between the wavelengths 5.3 and 12.5 microns (1870 to 800 em~') 


+ 


Fig. 6 Radiation circuit for parallel plates with no absorbing medium 
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has been subdivided into three regions. This procedure has been 
found to increase the accuracy of computations with an absorbing 
gas whenever a wide spectral region of zero absorption exists. 

The radiation network for parallel plates consists of three con- 
ductances in series (Fig. 6) and the solution of this system on a 
monochromatic basis is: 


— 


= — E, 2) (100) 


The net exchange between the two surfaces is obtained by in- 
tegrating equation (100) over all wavelengths using previously 
published procedures [3]: 


a @ v 
f s—2d = A,oT, f d 
~ f ( d (=) (101) 


The evaluation of equation (101) requires two numerical or 
graphical integrations, and it may be seen that a simple reci- 
procity relationship between the conductances from surface | 
to 2 and from surface 2 to 1 does not exist. This is in contrast with 
the gray radiation case where: 


E,, v 


The heat exchange for the gray radiation case is given by: 
= — T:*) 


The band energy approximation can be applied to this problem 
by replacing the integrals in equation (101) by the summation 
of individual band conductances. The band conductance in this 
case is approximated by using the arithmetic mean values of the 
emittances of each surface over the given band interval as indi- 
cated in equation (104), 


(103) 


MONOCHROMATIC REFLECTANCE 


60} — 
50}- 


30 
20}—- 


2 


TITANIUM ALLOY 


ANODIZED 


(104) 


This approximation is slightly easier to use in the analytical ex- 
pression for the conductance of this simple radiation circuit than 
the general expressions given in Part 2 of this paper for the multi- 
node network. This approximation is also used later on for the 
parallel surfaces with absorbing gas between them. 

Two examples of the evaluation of the net exchange are given 
in Table 5. In both cases, one surface was titanium alloy, C- 
110M which had been oxidized for 303 hours at 1000 deg F [1] 
and was assumed to have a temperature of 1260 F. The first 
example used a surface of the same material at 660 deg R. For 
these conditions, the following results were obtained: 


Substitution of these values in equation (101) gives a net heat ex- 
change of 2004 Btu/hr ft*. The computation with the gray radia- 
tion approximation used the computed values given in the litera- 
ture [1]. The results of these two methods of computation and 
the band energy approximation result is summarized in Table 5. 
Similar calculations using anodized aluminum for the colder sur- 
face are also presented in Table 5. The spectral data for both 
surfaces are shown in Fig. 7. 

The results tabulated in Table 5 demonstrate the magnitude of 
the errors involved in the band energy approximation and the 
gray radiation approximation relative to the exact calculation. 
The errors given are indicative only and apply to the particular 
cases considered. No significance should be attached to the par- 
ticular magnitudes shown as may be seen by the shift in the sign 
of the error in the two cases. However, in both instances the 
band energy approximation is significantly less in error. 
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Table 4 Band intervals selected for calculations with carbon dioxide 


absorption (12! 


Path Length = 1 foot 
Pressure (partial and total) = 1 atmosphere 
Temperature = 532°R 


Bominal Band Limite Band Limits Band Absorption 
Wavelength Wave Wavelength 
~wicrons- 
1s 5500 800 12.5 0.652 
600 «1000 12.5 10.0 
1000 1433 10.0 6.97 
1433-1870 6-97 5-35 
5.2 1670 «1960 5.35 5.05 0.022k 
4.8 1980 «2160 5-05 0.055 
260 62500 4-62 & 00 0.052% 
2500 4.00 62.87 
2.7 2.87 2-63 0.0570 
300s 2.63 2.10 — 
2.0 4750 §200 2.10 1.92 0.109 
$200 000 1.92 1.67 
1.6 6000 «(6550 1.67 1.52 
6550 6650 1.52 1.50 —_ 
lb 6650 7250 1.50 1.3% 0.0095 


Table 5 Comparison of band energy and gray radiation approximations 
with “exact” solution of radiation exchange between two infinite par- 
allel surfaces 


Surface 1 
600°F 


Surface 2 


Temperature 200°F 


4. Titanium Alloy (oxidised 303 hours at 1000°F) 


Method Heat Transfer Rate Per Cent Errer 
“Exact Solution® 2004 Btu/he te? 
Band Energy Approxiaat ion 1977 
Orey Radiation Approximation 1400 


B. Titanium Alloy C-L1OM (at 600°F) te Ancdised Alumioun (at 200°7) 
Method Heat Transfer Rate Per Cont 


“Exact® Solution 1646 
Band Energy Approximation 1626 
Grey Radiation Appraximation 22 +29 


The final example will consider the case of parallel surfaces of 
the same titanium alloy as in the previous example separated by 
carbon dioxide gas. The hot surface is taken as 1500 deg R, the 
cold surface at 530 R, and the gas at 530 R. This temperature of 
530 R corresponds to the temperature for the results of Howard, 
et al. [12] and is imposed upon the example by the lack of any 
other suitable data. The carbon dioxide is assumed to be at one 
atmosphere total pressure and the effective partial pressure path 
length product is equal to one atmosphere foot. 

The radiation network for these conditions using the gray 
radiation approximation is shown in Fig. 8(a). This circuit is 
readily solved to yield the desired heat-transfer rates. The net 
radiation loss from surface 1 is found to be: 

(+) €¢27-20(T,* — T;) 


1 — 


1 — 


+ 


(105) 


The net heat loss from surface 2 is of the same form as equation 

(105), but with the subscripts 1 and 2 interchanged throughout. 

The net heat gained by the gas, g,, is equal to the sum of the net 

heat losses from the two surfaces or, since A; = Az: 

( + — T,*) 
A/, 1 — rr2(71-2)* 


+ 


(106) 


1 — 
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Fig. 8(c) Radiation circuit for parallel plates with an absorbing medium: 
Gray radiation case 


WV 


dv 


Fig. 8(b) Radiation circuit for parallel plates with an absorbing medium: 
Band approximation 


For the conditions of this problem, T; = 7, so that the heat 
transfer to the gas becomes: 


( ) _ + T,*) 
A 


1 — 


The values of the radiant properties of the two surfaces and the 
effective absorption of the medium were taken from the litera- 
ture [l, 7] as: e, = 0.655, r; = 0.345, eg = 0.560, r, = 0.440, 
&-2 = 0.110, and 7 = ,-,0.890. The net exchange for these con- 
ditions was 1090 Btu/hr ft? between the hot surface and the gas. 

Equation (107) may be written on the basis of monochromatic 


energy exchange: 
(+) + rr - dy (108) 

A 0 T 1-2)? 

The difficulty of performing this integration has already been dis- 
cussed in Part 1. However, the band energy approximation may 
be applied to the original network and the various potentials and 
radiant properties evaluated over discrete frequency intervals 
[Fig. 8(b)]. The band energy approximation was computed using 
the band absorption data for carbon dioxide shown in Table 4. 
The band limits shown in this table are a function of the total 
pressure and temperature of the gas. For conditions other than 
those shown, these band limits would have to be adjusted. The 
availability of these data determined the conditions selected for 
the present example. Experimental data for carbon dioxide at 
higher temperatures and pressures are presently being measured 
and should be published shortly [4]. 

Utilizing the band energy approximation, the calculated heat 
transferred from the hot surface to the gas was 670 Btu/hr ft’. 
The result obtained using the gray radiation approximation was 
63 per cent greater than this. The result is in agreement with the 
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qualitative statement given by Hottel [10] to the effect the gray 
gas assumption will yield heat-transfer rates that are generally 
larger than the actual ones. An exception to this is the case of an 
enclosure with completely absorbing (black) walls where the gray 
assumption does not affect the results. Examination of equation 
(108) with e,., = 1, r,. = 0, and r,.. = 0 will show that under the 
conditions of “black” surfaces, only the characteristics of the gas 
are important and these characteristics may be shown to cor- 
respond to the measured values reported in the literature {10}. 


Discussion and Summary 


The procedure for a reiterative or relaxation method for solving 
a multinode radiation network has been described and a simple 
example presented. The exact solution, band energy approxima- 
tion, and gray radiation approximation have been compared for 
parallel surfaces with no absorbing intervening medium using two 
different surface combinations. In both cases, the band energy 
approximation was shown to be more accurate than the gray 
radiation method. A similar system with carbon dioxide gas as 
an absorbing medium between the parallel surfaces showed a 
difference of 63 per cent between the gray and the band energy 
approximations under the particular conditions selected. 

The examples comparing the two methods for approximating 
the exchange have shown that the band energy method is the 
more accurate. Physically, this should be the general result, 
since the actual spectral characteristics of the surfaces are more 
correctly approximated. The band energy method is more 
laborious and time consuming and, consequently, will probably be 
restricted to problems requiring more accurate results than ob- 
tainable with the gray radiation assumption. The application 
will also be limited by the availability of the data required for 
such computations. However, the significant fact which should 
be noted is that the inaccuracy of the gray radiation approxima- 
tion is in the method of analysis. Errors in the data used can be 
obscured by the errors of the method. 

The relaxation method for reducing the solution of the network 
to a systematic procedure is more generally applied. This pro- 
cedure is readily adaptable to either digital computers or hand 
calculation methods. Consequently, the analysis of complicated, 
multisurface enclosures can be performed with less restrictive 
limitations on the complexity of the system considered. The 
problem of matrix algebra, which can be difficult and intricate, is 
avoided by the relaxation or reiterative method. 


Nomenclature for Part 3 


A = area; A, = area of surface n 
E = emissive power; EZ,,, = monochromatic emissive power 
of surface n 
F geometrical shape factor (see Part 2); F,,, = geometrical 
shape factor of surface n with respect to surface m 
F factor relating geometry and radiant characteristics of 
two surfaces; $,,.-2 = monochromatic value for surface 
1 with respect to surface 2; F4,:-2 = band energy ap- 
proximation value of $,.-2 over the frequency interval 
Av, equation (104); :-. = value of $,.:-2 integrated 
over all frequencies, equation (102) 
radiosity; J, = radiosity of surface n 
total number of surfaces composing a specific enclosure 
absolute temperature; 7’, = absolute temperature of sur- 
face n; T, = absolute temperature of the medium en- 
closed within the enclosure 
emissive power of the medium; W,,,, = monochromatic 
emissive power of the medium involved in the ex- 
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change of energy between surfaces n and m 

emittance of the surface; ¢,,, = monochromatic emit- 
tance of surface n; ¢, = total emittance of surface n, 
equation (76), Part 2 

number of surfaces composing a general enclosure; for a 
specific enclosure m = 1, 2,3,...M 

a specific surface of the enclosure 

energy flow; q,.m = monochromatic energy flow between 
surface n and m; q,,, = total energy flow between sur- 
face n and m 

temperature 

absorptivity of the intervening medium; @,.nm = geo- 
metrical monochromatic absorptivity of the medium 
between surface n and m (see Part 2) 

frequency; Av = frequency interval for the band energy 
approximation 

Stefan-Boltzmann constant 

transmissivity of intervening medium; 7,,.m = geometri- 
cal monochromatic absorptivity of the medium lying be- 
tween surfaces n and Ty,nm + = am 
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Numerical Solutions to an Inverse 
Problem of Heat Conduction for 


Simple Shapes 


Numerical methods are presented for solving an inverse problem of heat conduction: 
Given an interior temperature versus time, find the surface heat flux versus time. 


The analysis is developed specifically for spheres; it applies to other simple shapes. 


The system is treated as linear, permitting use of the superposition principle. 
essence of the method is the numerical inversion of a suitable direct problem: Given a 


surface heat flux versus time, find an interior temperature versus time. 


Care is required in selecting a time spacing for, if it is chosen too small in relation to 


the conditions, undesirable oscillation results. 


examples. 


Introduction 


N INVERSE PROBLEM of transient heat conduction 
arises whenever a boundary condition such as surface heat flux 
versus time is sought from a knowledge of an interior condition 
such as temperature versus time. Such an inverse problem oc- 
curred as part of the co-operative Quenching Program [1, 2]! at 
Columbia University, where obtaining heat flux and surface tem- 
perature data was one of the purposes of a number of experiments 
in which temperatures were measured inside a body being cooled. 

In contrast with the inverse problem, the direct problem of 
transient heat conduction is relatively straightforward; in a 
direct problem the boundary conditions are known and the in- 
terior temperature is sought. Many cases appear in the literature 
{3}. To the author’s knowledge no general treatment of the in- 
verse problem has been given. For solving of the inverse problem 
a number of approximate procedures exist; some are: 


(a) Trial and error, including electrical analog. 
(b) Simplifying assumptions, when tenable, whereby the in- 
verse problem becomes a direct one. 


(c) Numerical inversion of a suitable direct problem. 


! Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division and presented at the 
Semi-Annual Meeting, St. Louis, Mo., June 14-18, 1959, of Tue 
American Society or MecHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, March 
19, 1959. Paper No. 59—SA-20. 


——Nomenclature 


Simplifying suggestions are presented, and the use of the methods are illustrated by 


All three procedures have been used by the author. An elec- 
trical analog computer method (a) was abandoned because in- 
strumentation problems delayed the work. A graphical proce- 
dure, based on a simplifying assumption (b) is in current use; it 
appears in the Appendix. The numerical inversion method (c) 
described in this paper has been in use since June, 1957 [4]. 

Several publications dealing with special simple methods (6) 
have appeared, many without mention of the restrictions implied. 
The first published use of a numerical inversion (c) seems to be in 
the recent paper of Mirsepassi [5], who presented charts for use 
with semi-infinite bodies. The inverse problem is dealt with as 
an illustration of graphical procedures, without mention of its 
derivation or possible limitations. These important questions are 
explored here in the development of a numerical procedure ap- 
plicable to simple finite bodies. 

Whatever method is used, the solution will suffer from an in- 
herent uncertainty. In a heat-conduction system the effect of 
boundary conditions is always damped at interior points, and the 
inverse problem involves basically the extrapolation of the 
damped datum to the surface, and sometimes differentiation at 
the surface. The damped information cannot be completely re- 
covered, but, as will be shown, this need not constitute a handicap. 

Statement of the Problem. For simplicity the heat-conduction 
system will be treated as linear with the consequent advantages 
of the superposition principle. Constant physical properties and 
no internal heat generation will be assumed. While the de- 
velopment presented specifically applies to spheres, it can be 
easily extended to slabs and cylinders. 


A = area, sq ft P = a function defined by equation Z = aconstant = peV/A, Btu/sq ft 
c = specific heat, Btu/(Ib)(deg F) (13) a = thermal diffusivity, sq ft/hr 
F = response ae defined by Q = total heat flow, Btu 8 = difference operator 
equation (1: 
g=a defined by equation 6 = deg F 
(16) z : \ = time spacing, hr or sec 
h = boundary conductance, Btu/hr, " ™ ~ particular values of r . y = index of series 
sq ft T = a function defined by equation Gini Teen te 
H = heat content, Btu (13) 4 : + 
j = index of series u = unit step function 7 = time, hr or sec 
k = thermal conductivity, Btu/(hr) V = volume, cu ft @ = heat flux out of solid, Btu/(hr) 
(ft)(deg F) x, = root of transcendental equation (sq ft) 


Nai Biot number = AR/k 


(12) 
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If the temperature of the body is initially uniform, and the tem- 
perature 6(r,, 7) at some interior point r; as a function of time is 
known, the statement of the inverse problem becomes 


r) (7,2) = 0; O<r<R (1) 


2 = 0 (2) 


Or, 0) = 8;, constant (3) 
Orn,7) =f(r); O<n<R (4) 


We seek the surface heat flux and surface temperature, ¢(T) = 
—kd0(R, r)/dr and &(R, 7). Note that, if either of these is found, 
the other is readily established, and the two results may be com- 
bined as @[A(R)). 

To facilitate the solution of the problem set forth in equations 
(1) to (4), we consider first a related direct problem set up in such 
a way that it may be easily inverted. 


A Direct Problem 

Suppose that equation (4) is replaced by equation (4a) which 
states: Instead of the interior temperature 6(r,, 7) being known, 
the heat flux ¢(7) is known. Now @(r, 7) is sought. 


d0(R, 7) 


= (4a) 


We can subdivide $(7) into elements as follows: 


Hr) = >> (5) 


2=0 


An element of temperature 6 (r, corresponding to each 
may be found by any convenient method—analytical, finite dif- 
ference, electrical analog, or direct experiment. The desired 
temperature is composed from its elements 


n 
Hr, 7) = 7) (6) 
27=0 
More generally, if F(r, 7)/z is the solution of equations (1), (2), 
(3), and (4a) with ¢(r) = —1 and 8; = 0, then from Duhamel’s 
theorem or the convolution theorem the solution of (1), (2), (3), 
and (4a) as written, i. e., with @(7) arbitrary, is given by 


(6, — Hr, = = — 7')F(r, 


dt Jo 


re) 
= f Hu) = F(r,7 — w)dp (7) 


For convenience the elements of heat flux @,(7) in equation 
(5) were chosen as steps 5¢(j) at equal time intervals \ and per- 
sisting thereafier; the argument denotes the time at which the 
increment in heat flux starts. 

Thus 


=nX (8) 


and 


Hr) = (nd) = 


7=0 


(9)? 


2A more conventional but also more cumbersome notation would 
n 


be o(mk) = — where u(r) is the unit step function. 
j=0 
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Other heat flux elements were tried—straight lines, exponen- 
tials—but showed less desirable properties. 

In performing the subdivision of equation (9) a good approxi- 
mation results if a heat balance is realized at each step: 


» a=} 


j=0 
The temperature response to the step 5@(0) is for a sphere 


10, — Or, ae + E 
10 


sin(zr/R) exp 
(r/R) sin z, 


2 
(11) 


where z, is the vth positive root of the equation 


z, cot z, —1=0 (12) 


The general form of equation (11) can be abbreviated 


F(r, 7) = t + P(r) + T(r, r) (13) 
where the “steady state’”’ part is r + P(r); the “‘transient’’ part, 
T(r, 7) becomes negligible after a short time. 

The behavior of other simple shapes is similar. The response 
to a step in heat flux for a cylinder [6] and a slab [7] can be put 
in the form of equation (13). 

Note that 6; — @(r, 7) in equation (11) is always proportional 
to 5¢(0) whether the latter be positive or negative. Further, if 
the step in heat flux is shifted in time, the corresponding tempera- 
ture is shifted in time by the same amount. Therefore, cor- 
responding to any single increment of heat flux 5@(jA) there is a 
temperature change at time nA, given by 


The total temperature becomes, corresponding to equation (6) 


3=0 


6, — Wr, nd) = (15) 


This could be regarded as a finite difference form of equation 
(7). 

A graphical illustration of the use of equation (15) appears in 
Fig. 1. The given heat flux, Fig. 1(a), is subdivided into ele- 
ments, Fig. 1(b), and the corresponding temperature elements, 
Fig. 1(c), are summed to yield the total temperature, Fig. 1(d). 

The evaluation of equation (15) is easily programmed either 
for a desk calculator or a digital computer. The response function 
F need be obtained only once for each value of r. Since the argu- 
ments of the terms 5¢ and F increase in reverse order, for desk 
calculation the functions are tabulated in reverse order and the 
sum obtained by accumulative multiplication. For a different 
value of nX the tables are shifted accordingly. 

An interesting graphical procedure for the evaluation of equa- 
tions of this type has been proposed by Mirsepassi [5]. 

Two additional shortcuts may speed the processing: 

The transient part of the response function in equation (13) 
becomes negligible after a certain time gX, called a long time in 
this context. For practical calculations gX would be given by 

ag\/R? = 0.2 (16) 
If the times of the calculation exceed gA, we split the summation 
in equation (15) into two parts, from j = 0 toj = n — g and from 
j=n—g+i1toj =n. In the first term Flr, (n — j)A) will 
contain only the steady-state part. Substitution of equations 
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Fig. 1 Graphical i 


1 
6, — Or, nd) = z Qli(n — g)A) + Flr, — anit 


+ >> ir, (n — (17) 


j=n-g+l 


A large number of terms have been replaced by two. 

In some portions of the problem a time spacing coarser than 
may save effort. This is done without altering the procedure by 
choosing integral multiples of \ and setting the heat flux incre- 
ments equal to zero at those times not of interest; then slight 
modification of equation (15) is needed. 

No limitation need be put on the time spacing since a direct 
problem is involved. Smaller time spacing increases accuracy 
but may cost more. 

In the next section, where the direct problem will be inverted, 
numerical considerations will set a lower limit on A, in part de- 
pendent on the location of point r;. 


Numerical Inversion of the Direct Problem 


Returning to the main problem stated in equations (1) to (4), 
examination of equation (15) suggests a solution. Where in (15) 
O(r, nX) was the only unknown, now all of the values of 5@(jA) are 
needed. After noting that F(r, 0) = 0, rearrangement yields 


+ — Hr, A) (18) 


As any value of d@(nA) depends on all the previous heat flux 
increments, the solution has to be generated step-by-step from 
time zero. 

Fig. 1 may be consulted for graphical interpretation. The in- 
formation at hand is 1(d), a measured temperature-time curve 
known at integral multiples of A. We would like to subdivide 
it into components as shown in 1(c), each of which will be a re- 
sponse to the heat flux increments 1(b), which in turn can be 
summed to give the total heat flux l(a). 
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of superposition principle in equations (15) and (18); (a) total 
heat flux, (b) heat flux elements, (c) temperature temperature 


(9) and (10) and rearrangement yields in place of equation (15) 


(d) total 


Now at the end of the first time interval A a temperature 
change 6; — @(r, X) has been observed. If it is agreed that a 
step change in heat flux occurring at time zero, 5@(0), caused this 
temperature change, then 6¢(0) has the value: 


59(0) = — Or, A)]/F(r, A) (19) 


At time 2A, 6¢(0) would produce a temperature change given 
by 6@(0)F(r, 2A)/Z; it differs from the observed value, 6; — 
&r, 2X). To account for the difference a new heat flux increment 
is started at time A, which in the time from A to 2A can produce 
the temperature change — @(r, 2A)] — 45@(0)F(r, 2A)/Z. 
5@(A) becomes 


5G(A) = — Or, 2A)] — 5G(O)F(r, 2A)}/F(r, A) (20) 


In this way the entire heat flux time curve can be generated. 
Equation (18) might have been obtained by induction by con- 
tinuing as in (19) and (20). 

Limitations of the Inverse Method. As pointed out in the intro- 
duction, the calculation is based on truncated data which have 
been damped. Questions of the effect of experimental errors are 
considered separately and will be dealt with in the following. 
Improper selection of A—too small for a given set of condi~ 
tions—will lead to an undesirable oscillation. The oscillation 
shows up as an error which alternates in sign and may keep 
increasing in magnitude at each step. Even if oscillation is 
avoided, the solution will be approximate (determined within a 
certain range) as will be illustrated. 

The oscillation stems from the character of the system. It 
will be recalled that heat conduction systems exhibit damping, 
i.e., at interior points the effect of any surface disturbance is de- 
layed and diminished. Thus the effect of a heat flux change oc- 
curring at time (n — 1) may be undetected in the data at time 
nX, if X is too small, and the calculation yields zero change in heat 
flux. In calculating the next heat flux element, the error is com- 
pensated, overcompensated if \ is too small. Thus in succeeding 
times the error oscillates. 

Whether the oscillation can become unbounded (unstable) and 
under what general conditions is not known at present. Re- 
gardless, it may become much larger in amplitude than the cor- 
rect result; in such a case the answer is useless. 

A general criterion for the optimum value of \, or more gen- 
erally a\/R?, as a function of r/R is under investigation. Usually 
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Fig. 2 Results of performance test using equation (18) compared with 
exact result. Two different time spacings. 
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a few trials will establish the optimum. As might be expected 
from the foregoing, the smaller the value of mn, the larger \ needs 
to be. Empirical results seem to show that Aop: is at or near the 
inflection point in the response function F(r, 7). 

Oscillation may sometimes be reduced when the choice of X is 
not open, by proceeding stepwise in space. I.e., consider rx-—R > 
T: > r:—to be a surface, the “boundary condition’’ found there 
becomes the datum for the second step, etc. But such bootstrap 
procedures do not improve the accuracy of the result sufficiently 
to warrant the effort. 

The arbitrary choice of heat flux element as a step and time 
spacing \ leads to a particular answer. A somewhat different 
answer would be obtained from a different choice. Fig. 2 shows 
the result of a performance test. The particular conditions 
represent a 1.0-in. radius silver sphere (k = 207 Btu/(hr)(ft)(deg 
F), a = 5.33 sq ft/hr) subjected to the rather severe prescribed 
heat flux denoted by the dotted line. The temperatures at r = 
0.90 in. were obtained analytically and became the statement of 
the problem. Two different time spacings, \ = 0.025 sec, \ = 


a 


0.05 sec, were used to calculate the heat flux, which is plotted. 
The smoothed end results agree fairly well with the prescribed 
condition; the worst discrepancies occur at zero time and at the 
peak heat flux. It is shown that the results do differ somewhat. 
Further, for the shorter time spacing an oscillation is seen. 

The effect of a single experimental error in a particular set of 
circumstances is readily determined by use of equation (18). 
Since the system is linear, the effect may be found directly after 
specifying a particular error by carrying out the operations indi- 
cated. Again the magnitude of the error depends on the choice 
of conditions. It diminishes as the location r; nears the surface. 

The foregoing all points to one conclusion: When the choice of 
the location of the temperature probe is free, all other factors 
being equal, the best place to probe is close to the surface. 

Practical Details. Some suggestions as to procedure have been 
presented in the section Direct Problem. The use of variable 
time spacing and the shortcut for long time—equation (17)— 
apply directly to equation (18). 

The final heat flux, in the form of steps, may be smoothed by 
invoking a heat balance, equation (10). The smoothed values 
become the mid-points of the vertical segments; for more detail 
the mid-points of the horizontal segments may also be used. 

An entire problem need not be handled by the inverse method if 
a simpler method can be applied. As the foregoing treatment 
implies an initial condition of uniform temperature (and zero heat 
flux initially), a minor revision is required: 

According to equation (13) the response to a stationary or 
slowly changing heat flux becomes linear in time after a time gX 
given in (16). When the region of interest is preceded by a suf- 
ficiently long period of approximately constant cooling rate, the 
constant heat flux value that would produce this behavior is in- 
troduced to start the calculation. In this way the correct tem- 
perature profile is established. Since the inverse method deals 
with changes in temperature, such a procedure is possible and 
permits the calculation to be started at a suitable time other than 
zero. 

illustrative Example. For a specific application of the methods 
presented a quenching experiment is selected. The experiments 
themselves appear elsewhere [8]. A silver sphere, 1.0-in. radius, 
surface roughness 140 microinches, was heated to 1600 F and 
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curve in Fig. 4. 


Journal of Heat Transfer 


FEBRUARY 1960 / 23 


to 

aes 

= 
a 

Pay 

} 2 

800 1000 1200 i400 


‘ quenched in a “fast oil” at 110 F. The oil was initially quiescent 
in a deep 38 X 31-in. tank, and the sphere was submerged by 2.5 
e— in. Direct measurement of surface temperature would interfere 


seriously with the phenomena being observed. For this illustra- 
tion temperature was recorded at r; = 0.9 in. 

In Fig. 3, showing the final result plotted as boundary con- 
ductance versus surface temperature, the points shown as circles 
were calculated by a simple graphical procedure (appendix) since 
they correspond to low and slowly changing heat fluxes (Ngi < 
0.04). 

The triangles were calculated by the inverse method. The 
measured temperature 0(0.9, r) appears in Fig. 4, from which the 
heat flux $(7) in Fig. 5 was calculated. Application of the direct 
problem using $(7) in equation (15) gave the surface temperature 
&(R, 7) in Fig. 4; the results were then correlated as Fig. 3. 

In the example given the effect of an error in temperature was 

calculated with equation (18). With an isolated temperature 

(0.3, T) error of 1 F, using time spacings of 0.05, 0.1, and 0.2 sec, the cor- 

responding maximum errors in heat flux, in millions of Btu/hr sq 
ft are 0.050, 0.024, and 0.013. 

A complete calculation as illustrated requires approximately 


two man-days using a desk calculator. It is estimated that the 
~~ use of a digital computer such as IBM 650 would replace one day 
by a few minutes. 


TIME, SECONDS This work was performed at the Heat and Mass Flow Analyzer 
Laboratory at Columbia University; the following named in- 

dividuals are present or past members of the laboratory. Dr 
time. H. G. Landau generously gave advice at many stages of the in- 
vestigation. He was also responsible for some of the mathemati- 

cal developments. Mr. J. L. Sullivan performed much of the 

25 numerical work, aided by Mr. F. Ruf and Mrs. A. Heller. The 
author is indebted to Prof. V. Paschkis and Prof. C. F. Bonilla 
for their encouragement and patience, particularly during the 
early phases of the work. Grateful acknowledgment is made for 
the support of the following organizations: Caterpillar Tractor 

20 Co., Chrysler Corp., Curtiss-Wright Corp., Engineering Founda- 
tion, Frankford Arsenal, General Electric Co., Gulf Research & 


TEMPERATURE DIFFERENCE °F 


amining the relative magnitude of the thermal resistances of the 
fluid and solid. If the fluid controls, the Biot number Nj; will 
be low and, as a rule, will change gradually. The body will be 
practically isothermal. For example, when Ng: = 0.04, the dif- 
os ference between surface and center temperature (referred to 
ambient) will be less than 2 per cent. The temperature, now (7), 
can be taken proportional to heat content H(r) and 


¢ Development Co., Ladish Co., Libbey-Owens-Ford Glass Co., 
2 SKF Industries, Texas Co., Timken Roller Bearing Co., Westing- 
«. house Electric Corporation. 

is 

5 APPENDIX 

4 

" Simple Method for Low Heat Flux 

+ | Lo a Considerable computation can sometimes be spared by ex- 


= —(1/A)dH/dr = —(peV/A)dO/dt = —ZdO/dr (21) 
or, if desired, 


h(t) = = —Zd O/dr (22) 


A graphical method for finding the time derivatives is described 
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Fig. 5 Ilustrative example. Calculated surface heat flux ¢(7) versus 
time. 


be 


by Hartree [9]. It consists of plotting @ or In @ versus time and 
obtaining the normals to the curve using a straight edge with a re- 
flecting surface. 98 per cent of the data recorded by the author can 

0.4k 
be processed in this way. It can be used whenever =| {= 
with accuracy of usually 1 per cent, and at somewhat more rapid 
cooling with decreasing accuracy. 
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DISCUSSION 
T. J. Mirsepassi* 


The author is to be congratulated for his interesting presenta- 
tion and useful study of an important problem of heat conduction. 

On numerous occasions during recent years we have used this 
technique both numerically** and graphically** and have found 
it, quite in agreement with the author, a fast and accurate method 
especially adaptable to digital computers. Some cases were also 
calculated by the well-known arithmetic method of Schmidt. It 
was observed that the finite difference form of Duhamel’s integral,’ 
having the following advantages over the Schmidt method, is far 
better-suited for engineering applications: 


1 Simplicity of calculation scheme. 

2 Fewer arithmetical operations. 

3 Independence of operations in direct problems, i.e., possi- 
bility of calculating temperature at a point inside the solid at any 
time without necessity to calculate for previous time intervals. 

4 Lower probability of errors in calculation. 

5 Errors, being local, are detected readily by comparison with 
the adjacent results. 

6 Adaptability to desk calculators; i.e., the operator can carry 
out accumulative multiplications and register the final result 
only, while in Schmidt’s method a large number of intermediate 
calculation must be registered. 


The technique is by no means limited to simple shapes.‘ ** 
The temperature at the point of interest for an arbitrary boundary 

* Director of Research Computer, Space Technology Division, 
Aerojet-General Corporation, Azusa, Calif. Mem. ASME. 

4T. J. Mirsepassi, “Special Report No. 1375,"" Aerojet-General 
Corporation, Azusa, Calif. 

*T. J. Mirsepassi, ‘‘Heat-Transfer Charts for Time-Variable 
Boundary Conditions,” British Chemical Engineering, March, 1959. 

*T. J. Mirsepassi, “Graphical Evaluation of a Convolution Inte- 
gral,”” Mathematical Tables and Other Aids to Computation, July, 
1959. 
Equation (2),* equations (3) and and equation (3).* 
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function can be readily calculated in an irregular configuration, 
provided that the temperature history at that point is known 
when the surface (or the ambient) temperature undergoes a step 
change. Since such problems are customarily solved by an analog 
computer, this technique results in (a) considerable saving of 
analog computer time, and (b) the possibility of solving response 
problems when the analog computer is no longer available. 


E. M. Sparrow® 


Evidence attesting to the interest in the inverse problem of heat 
conduction may be found in an earlier publication on this subject 
entitled, ‘‘A Method for the Experimental Study of the Process of 
Heating a Solid Body,” by N. V. Shumakov.* The problem which 
is posed there is the same one as is considered in the present 
paper, namely, to find the heat-transfer rate at the surface from 
unsteady temperature measurements in the interior. The tech- 
nique used in the earlier paper is called The Method of Successive 
Intervals, and it is applied to the heat flow in a slab. Experi- 
mental temperature measurements are made, and these are used 
in conjunction with the analytical method to predict the surface 
heat transfer. 

It would be interesting if the author would compare his ap- 
proach with that of the earlier paper, pointing out the ad- 
vantages or disadvantages of each. Also, it might be useful for 
him to apply his method to the slab to see if his surface heat-trans- 
fer predictions are consistent with those computed by Shumakov. 


Author’s Closure 


The interesting and valuable comments by the discussers are 
sincerely appreciated. 

The Shumakov paper cited by Dr. Sparrow escaped attention 
until shortly before the oral presentation. The author is pleased 
to elaborate his verbal remarks. 

Shumakov limits his treatment to a minimum dimensionless 
time spacing or Fourier number, N»,*, of 0.5. For such long times 
the transient part of the response function F (r, 7) in 
equation (13) becomes negligible, and the remaining steady 
state part becomes linear in time. The author feels that an 
Ny.* of 0.5 is too conservative; 0.3 is a safe limit, and in many 
cases 0.2 is satisfactory. Such a restriction does not apply to 
the present treatment. The author has been able to use dimen- 
sionless time spacing as small as 0.002, a hundred times shorter. 

Equation (17) of this paper is readily specialized to yield the 
Shumakov treatment. We choose the time spacing long enough, 
ie., 


ar/R* > (23) 
which, when compared with equation (16), permits us to put 
g=1 (24) 


Using this in (17) and noticing that F(r, 0) = 0, the summation 
vanishes. With the definition of Q(n A) in (10) one obtains: 


n-2 
Z(0,, — Or, = — LAIF(r, A) +A 2: @(jA) (25) 
7=0 


Except for differences of notation this is identical with Shu- 
makov’s equation (23). 

Since the treatment is the same for “long times,’’ Dr. Sparrow’s 
suggestion to redo the Shumakov calculation was not followed, as 
identical results should be obtained. However, a much simpler 


®* NASA, Lewis Research Center, Cleveland, Ohio; at present 
Professor of Mechanical Engineering, University of Minnesota, 
Minneapolis. Assoc. Mem. ASME. 

* Soviet Physics—Technical Physics (translated by American In- 
stitute of Physics), vol. 2, 1957, p. 771. 
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calculation is possible. If Shumakov’s temperature-time data are 
replotted on semilog co-ordinates, one obtains very nearly a 
straight line, indicating a constant Biot number, Ng; = 0.42. 
The corresponding heat flux versus time relation is within 10 per 
cent of the Shumakov prediction. 

A more instructive comparison would perhaps be obtained by 
attempting to recalculate the results of Figs. 2 and 4 of this 
paper by the Shumakov technique. In the examples used, a 
Fourier number of 0.5 corresponds to a time spacing of 2.4 
seconds. Thus the calculation leading to Fig. 2, which spans only 
0.3 seconds, could not be done. Turning to Fig. 5, based on Fig. 
4, we choose for convenience \ = 2.5 seconds and extrapolate the 
curves back by one second, permitting two time steps from 12 to 
17 seconds. We obtain for the first time step a heat flux of 0.098 
< 10° Btu/hr, ft*, quite close to the correct value. From 14.5 to 
17 seconds the result is 0.57, a correct average value, but the 
essential feature of the physical process has been lost. 

The pertinent comments of Dr. Mirsepassi, based on his ex- 
perience with similar problems, are a welcome addition to the 
paper. Especially interesting is his suggestion that the techniques 
described can be applied not only to simple shapes, but to irregu- 
lar configurations as well. One qualification appears warranted. 
The pertinent boundary condition—heat flux or surface tempera- 
ture—must be expressed as a product of two functions, one of 
position alone, the other of time alone, e.g., 


6 or = f(z, y, z, 7) = filz, y, z)fe(T) on S (26) 


where 0/0n indicates differentiation in the direction of the nor- 
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mal to the surface. In such a case Duhamel’s theorem” leads to 
a convolution integral, in the usual sense, as in equation (7). 
There Dr. Mirsepassi’s suggestion would apply. If, however, 
the boundary condition is not separable as in (26), the resulting 
form of Duhamel’s theorem would not readily lend itself to this 
type of approximation. 

It is of interest that in attacking the inverse problem, where 
heat flux was sought, neither Dr. Mirsepassi, using his graphical 
technique on a semi-infinite geometry, nor Shumakov encountered 
the numerical oscillation reported by the present author. For 
the long times used by Sumakov, one would not expect oscilla- 
tions. It is possible that a semi-infinite body, where the response 
function has a different form, would not behave the same as the 
finite geometries considered in this paper. 

In reply to some questions about alternate formulations for the 
inverse problem: More elaborate approximations to either the 
unknown function or the entire integrand in equation (7) are 
possible. The Laplace transform of the same equation could be 
inverted numerically. When ¢(7T) is the unknown, equation (7) 
is classed as a Volterra integral equation of the first kind. Fora 
similar integra] equation Hilderbrand'! presents a least squares 
technique which ought to be of value in reducing numerical 
oscillation. The present approach is recommended by its relative 
simplicity, as Dr. Mirsepassi has also noted. Even simpler 
methods such as an extension of the analysis recently presented 
by Elrod"? are worth while whenever they apply. 

” Reference [3], pp. 21-22. 

11 F. B. Hilderbrand, “Methods of Applied Mathematics,” Pren- 
tice-Hall, New York, N. Y., 1952, pp. 452-459. 


12H. G. Elrod, “Improved Lumped Parameter Method for Tran- 
sient Heat-Conduction Calculations,"” ASME paper No, 59—HT-28. 
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Effect of Thermocouple Cavity on 


JAMES V. BECK 


Associate Scientist 


Heat Sink Temperature 


An analysis is made of problems associated with prediction of and correction for tem- 


H. HURWICZ 


Senior Staff Scientist 

Research and Advanced Development 
Division, AVCO Corporation, 
Wilmington, Mass. Assoc. Mem. ASME 


perature disturbance created by thermocouples placed beneath a surface of a heat sink 
(or calorimeter) exposed to heat flux during re-entry. Two of the important factors 
affecting temperature measurement are discussed: 


(a) The disturbance created by the thermocouple itself, “‘hot spot.” 
(b) The fact that the thermocouple has to be placed at some distance from the heat flux 


surface and thus not measuring the surface temperature. 


The magnitude of the surface hot spot caused by the presence of the thermocouple is 
determined, and optimum location of the thermocouple is found where the undisturbed 
surface temperature may be read with least over-all error. 


HE FLIGHT TEST EVALUATION of a re-entry vehicle 
performance requires the determination of the aerodynamic heat- 
ing at the heat sink surface. Thermocouples may be placed near 
the surface of the heat sink to obtain a temperature history which 
then may be used to find the surface heat fluxes caused by aerody- 
namic heating. This analysis is concerned with the case in which 
the thermocouples possess a much lower conductivity than the 
surrounding heat sink; a valid approximation (which gives the 
maximum temperature disturbance) is to replace the thermo- 
couples with an insulated cylindrical void. The thermocouple 
void would disturb the normally one-dimensional heat flow and 
in so doing would increase the surrounding temperatures. The 
surface heat flux determination would be affected by the increased 
thermocouple temperatures and also by the fact that the thermo- 
couple has to be placed at some distance from the heat flux sur- 
face. Several analyses of these effects are presented to cover the 
geometries of interest. The results may be applied to produce 
more accurately the experimental data which are to be used to 
calculate heating rates. 

The thermal effect of a void or an imperfection near the heat 
sink surface would be obtained from the same boundary condi- 
tions. The presence of the void would result in a hot spot at the 
surface of the heat sink. Though the results are for a cylindrical 
void which extends from just under the surface to the rear of the 
heat sink, they may be used as the (close) upper bound of the ef- 
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fect of a void which does not possess such a large depth. The 
determination of the magnitude of the surface hot spots due to 
voids may be used to formulate heat sink material fabrication 
specifications and to indicate permissible distribution and size of 
imperfections, voids, and occlusions near a heat sink surface. 


Discussion 


A thermocouple embedded in a heat sink, which has a much 
higher conductivity than the thermocouple, disturbs the sur- 
rounding temperatures. The maximum temperature disturbance 
would be caused by replacing the thermocouple with an insulated 
cylindrical void placed just inside the heat sink and perpendicular 
to the heat sink surface, Fig. 1(a). 

The effect of the assumption in which the thermocouple is re- 
placed by a void is assessed at the end of this paper. Until the 
discussion of this assumption the thermocouple effect and void 
effect may be considered identical and a reference to an “upper 
bound”’ or “lower bound’”’ will be considered in relation to whether 
the effect is respectively greater or smaller than the void effect. 
The upper or lower bounds would be caused by boundary condi- 
tions which approximate the disturbance caused by the void. 

The boundary conditions [see Fig. 1(a)] for the insulated cylin- 
drical void problem are 
oT (0, r, t) 

dz 


1) k — qt) 


oT(L, r, t) 
oz 


2) 0 (r > R) 


o7(z, R, t) 
dT (£, r, t) 
— 
5) T(z, r,0) = 0 


3) 0 (E<z<L) 


4) 0 (0<r<R) 


Nomenclature” 


E = distance of thermocouple tip 
below sink surface 


g=R/E dimensionless ratio of posi- 


tion of thermocouple tame 


k thermal conductivity of heat 
sink material 

L thickness of heat sink 

de heat flux at thermocouple 
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surface heat flux 
radius 
thermocouple radius 


response time 

temperature 

thermocouple tip tempera- 
ture 

surface temperature above 

tip thermocouple 


temperature at EZ if thermo- 
couple were not present 

surface temperature if ther- 
mocouple were not present 

2-co-ordinate 

dimensionless position of 
thermocouple tip 

= thermal diffusivity 


—— = dimensionless time 
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Fig. 1(a) Boundary conditions and dimensions for determining 
perature disturbances 


tem- 
caused by cylindrical void in a finite heat sink 


Fig. 1(c) Boundary conditions for the thermocouple error in semi- 
infinite slab 


for the heat conduction equation in cylindrical co-ordinates 
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A literature Review. An upper bound of the temperature above 
the thermocouple junction is given by the fine paper by Masters 
and Stein [1]** who effectively used the superposition principle. 
The heat sink with a surface heat flux was first considered with an 
imaginary void [Fig. 1(6)]. The heat flux ¢,(t) which passes 
through the z = E surface of the imaginary void was determined. 
Consider now the heat flux ¢,(t) introduced in the negative z-direc- 
tion at Z = E [Fig. 1(c)] in an insulated semi-infinite body. If 
the heat fluxes at the boundaries in Figs. 1(b and c) are added, it 
may be observed that the boundary conditions for the thermo- 
couple problem for a semi-infinite body are obtained. The tem- 
perature disturbances due to the void may then be obtained using 
the boundary conditions shown by Fig. l(c). This change in 
boundary conditions still did not permit the analytical solution of 
the problem immediately. By allowing, however, the heat re- 
jected by top of void to diffuse only in the infinite plate of thick- 
ness E, M and 8 were able to obtain an analytical solution. 

In the M and § paper, the effect of several smooth analytical 
time-variable heat fluxes applied at the sink surface was ob- 
tained by determining the variation of these fluxes at the E-level 
in a semi-infinite body and then applying these latter heat fluxes 
at the thermocouple tip. The results at the sink surface hot spot 
(point O in Fig. 1) for a uniform heat flux for g( = R/E) of 0.25, 
0.5, and 1 are shown by Figs. 2 and 3, and are quite similar to 
those obtained with variable heat inputs. It should be noted that 

2? Numbers in brackets designate References at end of paper. 


* Subsequent reference to Masters and Stein will use the abbrevia- 
tions M and 8. 
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Fig. 1(b) Boundary conditions of one-dimensonal heat flow in a semi- 
infinite slab with imaginary void 


in Fig. 2 the per cent change in hot-spot temperature has been 
divided by g? and in Fig. 3 the hot-spot temperature has 
been divided by g. Thus it is observed that per cent change in 
temperature due to the hot spot varies approximately as g* while 
the hot spot itself varies as g. Though the maximum error for g = 
0.25 given by M and S is 1.5 per cent (which is quite small) it 
will be shown that usually the most advantageous position for the 
thermocouple is where g = 1.1 when the maximum per cent error 
given by M and 8 would be about 23 per cent which is rather 
large. The authors obtained the M and § curve for g = 1 using 
the given equation and found that the solution could not conven- 
iently be extended beyond Tyg = 250 since a large number of sensi- 
tive subtractions were involved. 

In Fig. 2, it may be noticed that the maximum per cent error 
occurs approximately at at/E* = 8. For a typical thermocouple 
in a copper heat sink the maximum per cent difference would 
occur at about '/2 second when the temperature and heat flux 
level would be low. The determination of the thermocouple hot- 
spot effect is desired, however, at the time of maximum surface 
temperature which might occur several magnitudes of time later 
after the beginning of heating. From Fig. 3 the absolute hot- 
spot temperature error as given by M and § is seen to grow ap- 
preciably in time and to have no maximum. The use of the M and 
8 results (Fig. 3) for time variable heat inputs which are discon- 
tinuous—such as might be caused by a change from laminar heat- 
ing to much higher turbulent heating rates—is not apparent and 
may result in somewhat larger upper bounds. 

It will also be shown that usually the most important state for 
the thermocouple disturbance is not a transient but a “steady’’ 
state which is said to exist when the disturbance is approximately 
proportional to the surface heat flux. This steady state is only 
possible in a semi-infinite heat sink but will be shown to be closely 
approximated in a finite heat sink; that a steady state exists 
should be emphasized since it is a fundamental point of the 
present paper and distinguishes it from the M and 8 work which 
emphasized the transient aspects. As a result of the steady-state 
disturbance concept, a “null’’ position is located where the 
thermocouple reads essentially the undisturbed surface tempera- 
ture; the location of this position was not given by M and 8. 
Bussell [2] also investigated and obtained a value for the hot spot 
with R/E = */, using a passive element analog; the results are 
shown by Figs. 2 and 3. 

B Calculations and Results. To cover the complete range of in- 
terest, two analytical and two machine solutions were obtained. 
One analytical solution, which is particularly applicable for small 
voids (g < 0.5), employs a point heat source; transient and 
steady-state solutions are given. The other analytical solu- 
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Fig. 2 Comparison of transient thermocouple hot-spot temperature results relative to the surface temperature 
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Fig. 3 Transient thermocouple hot-spot temperature results for a semi-infinite solid 


tion (steady state only) gives reliable results for large g’s(g > 2); 
due to space limitations the solution will be outlined. The machine 
results have been obtained for the intermediate void or thermo- 
couple sizes (0.5 < g < 3). Both steady and transient solution 
were obtained on a digital computer. All these solutions supple- 
ment one another and have their particular range of greatest ac- 
curacy. The results exhibit substantial agreement among them- 
selves and all possess a steady-state solution in distinction to the 
data available in the literature. The point that each of these 
solutions has yielded a steady-state solution is important since 
the steady-state thermocouple errors are shown to be applicable 
to transient thermocouple measurements. 

1 Point Heat Source Analytical Solution. As noted in the previous 
discussion of the M and S paper, the thermal effect of a void in a 
heat sink may be simulated by a heat input at the circular void 
surface equal to the negative of the undisturbed heat flux at that 
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depth [Fig. 1(c)].. As the void becomes small in radius, the top 
of void approaches a point. In this case the thermal effect of the 
void would be approximated by a heat point source of strength 
equal to the heat that would be reflected at the tip of the void. 
This approach is very similar to that used by M and S who em- 
ployed a disk heat source; the solution may not, however, be ob- 
tained directly from their treatment since a disk heat source in the 
limit (as R approaches zero) does not approach the point source 
solution. This may be proved by noting that at the center of the 
disk source the temperature is always finite while at the point 
source the temperature is not bounded. 

The void effect using the point heat source would be deter- 
mined mathematically by solving Equation (1) for a point source 
of strength 1*R*g, between two parallel insulated planes. 4g, is 
the undisturbed heat flux that would be caused at EF by a surface 
heat flux, go. To simplify initially the determination of the void 
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effect, the source strength at E could be considered to be constant. 
For a uniform surface heat flux (qo), ¢, would take a small but 
finite time to reach ite maximum value. For a heat sink of finite 
thickness ¢, will approach 


L-E 


Ge = Go 


In reality, this heat flux is variable at a point in a heat sink 
which has a uniform heat flux at z = 0 and is insulated at z = L 
[4]. Equation (17), however, may be used with the principle of 
convolution to determine the transient hot-spot temperature for 
various values of Z. Fig. 4 gives a comparison for Z = 0.1, 0.01, 
and zero. It may be noted that initially all the Z-values give the 
same dimensionless temperatures; at later times, however, the 
ternperatures for Z = 0.01 and Z = 0.1 fall below the values for a 
semi-infinite body since during this time the heat input at the 
point source is reduced for larger Z’s [Equation (2)]. The 
situation is reversed later due to the presence of the rear sink sur- 
face. It is very important to note that for a semi-infinite body 
(Z = 0) a steady-state temperature disturbance is reached and 
is closely approximated for a void very near the surface for ordi- 
nary times (rz < 10° for Z = 0.01). 

A further conclusion may be drawn from Fig. 4. The hot-spot 
error grows rapidly at the beginning and reaches about 80 per cent 
of the steady-state value in at/E*? = 30 which will be called the 
response time. For ordinary thermocouple placement this dimen- 
sionless response time indicates a very rapid temperature re- 
sponse at the void. For example, the response time in a copper 
heat sink for a 0.02-in. radius thermocouple, 0.02 in. below the 
surface, is computed to be about 0.08 second. This suggests that 
the steady-state calculation for a thermocouple located near the 
sink surface is sufficient for most cases. 

For a semi-infinite sink heated with a uniform surface heat flux 
fand thus time-variable at the void, Equation (18)], the tem- 
perature disturbance caused by a point source is given by Equa- 
tion (19). 

For the hot spot (0, 0) the transient temperatures would be 
from Equation (19) 


To Tre 1 1 
= —erle —= 3 
~ 2° @) 


Equation (3) may be used to obtain the variation of the hot-spot 
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Fig. 4 Effect of the rear surface of the finite heat sink upon the hot-spot temperature 


2 


temperature relative to the surface temperature. The latter 
may be found for a semi-infinite slab with a uniform heat input 
from Reference [4] to be 


/k g? Vx (4) 


To Toe 1 (5) 
which is shown by Fig. 2 and may be compared with previous re- 
sults. 
The steady-state temperature increase at any point is given by 
Equation (20). 
The steady-state temperature distribution at the surface is of 


particular interest as it may be used in determining the interac- 
tion of closely spaced voids; it is given by 


(6) 
r 2 
2 


The heat sink surface steady-state temperatures given by Equa- 
tion (6) are shown by Fig. 5. 

The temperature of the hot spot is of particular interest as its 
determination was one of the primary objectives of this investiga- 
tion; from Equation (3) it is (see Fig. 6) 


(7) 


Equation (7) would be expected to apply to small values of g 
(small radii for a given void depth). 

2 Large g Analytical Solution. For large g ratios the geometry may 
be broken into two parts: The finite cylinder of radius R and 
thickness EF immediately above the void and a semi-infinite hol- 
low cylinder of radius R. An exact solution may be obtained if 
the two cylinders have their boundary conditions properly 
matched atr = Randz = O0to£F. This is not convenient to do; 
however, for large g’s it can be deduced that the temperature dis- 
turbance along this interface will be approximately constant. 
Klamkin [5] solved analytically the finite cylinder problem 
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Fig. 5 Steady-state numerical and point heat source solutions for the heat sink surface 
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which has the uniform heat flux g- at the z = E surface and a 
constant temperature at r = R. The heat flux was then found at 
r=Rfromz = £. 

The semi-infinite cylinder also was solved analytically by 
Klamkin [6] for a uniform heat flux at r = R from z = 0 to E; 
the temperature disturbance at the point r = R and z = 0 only 
was obtained using a close upper bound approximation to 
evaluate a certain integral. With this later solution it was then 
possible using the principle of superposition to introduce the space 
variable heat flux at r = R found from the disk. By then adding 
the temperature disturbance at r = R, z = 0 in the semi-infinite 
cylinder to the disturbances found in the finite disk, the steady- 
state temperature disturbances due to the void are obtained. 
The results for the hot spot, the thermocouple tip temperature 
disturbance, and the null point are shown by Figs. 6, 7, and 8 re- 
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spectively; these results constitute an upper steady-state bound 
of the temperature disturbances. 


3 Digital Solutions. The solution of the thermocouple problem 
was correctly thought by M and § to be somewhat difficult using 
a digital computer. Several considerations led to this conclusion, 
The thermocouple is a rather small device and the area above it 
occupies a very small part of the volume of the heat sink; this 
leads to a very small cell in the vicinity of the thermocouple which 
results in a large number of finite difference cells and in short 
computing intervals. 

Another difficulty is adequate simulation of the boundaries. 
A finite section of the heat sink must be considered in a numerical 
solution while from the point source solution, it has been shown 
that the thermal effects are not strictly local. The authors have 
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successfully adapted digital techniques to yield both transient 

and steady-state numerical solutions for this problem. 
(a) Transient Digital Solution. In a transient digital solution, 03 F atve* 

a “direct’’ approach (superposition not used) may be employed. Sy 
The boundary conditions used in this solution were a uniform 02 

heat input at the front surface of a homogeneous slab with an in- 

Since a finite section of the heat sink only can be 

considered in this approach, a sufficiently large cylindrical section 

of the heat sink was considered not to affect the temperatures 

about the void. The difficulties mentioned by M and 8 were or 

overcome by (1) using small cells only near the void and larger 

cells elsewhere and (2) by not attempting to go to large physical -O! 

times which actually were not required since the thermal re- 

sponse of the points about the void are rapid and also since the 
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Fig. 9 Transient thermocouple temperatures for a semi-infinite solid 
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final values (steady state) for a semi-infinite sink are given by a 
separate numerical solution. The transient numerical results for 
g = */s, 1, and */; for the hot spot are given by Fig. 3 until the 
radial boundary temperature disturbance would increase 2 per 
cent according to the point source solution. 

The transient variation of the temperature at the thermocouple 
tip relative to the undisturbed surface temperature is shown by 
Fig. 9. 

In locating a thermocouple it would be desirable \ if the hot-spot 
temperature were not excessive) to place the thermocouple so 
that the thermocouple reading would require a minimum correc- 
tion. Provided the response of the temperatures about the void 
were rapid compared to the variations of the surface heat input, 
the thermocouple could be placed at the null point where 7, — 
Tie ™ 0 (see Fig. 9). At at/E* = 30 the null point is within 10 
per cent of its steady-state g value of 1.0. 

The response time for both the null point and hot spot might 
then be given by the dimensionless time 

ad, 
30 (8) 

As given in a previous example, the response time ¢, is very 
rapid (0.08 sec) for a typical placement of a thermocouple in 
copper. In cases when this occurs, the thermocouple may be 
placed at the null point to read the surface temperature and the 
steady-state solutions may be used to determine the temperature 
errors about the thermocouple. If the surface heat flux varies 
widely in the response time, then the transient results could be 
used (see below). The “transient” thermocouple corrections are 
not usually required, however, since the response time is usually 
quite small. 

(b) Steady-state Digital Solution. As has been observed the 
steady-state solution of this problem is frequently more im- 
portant than the transient solution. One use of the numerical 
steady-state solution is to determine the steady state null point 
which cannot be accurately determined using the known analyti- 
cal solutions. The steady-state numerical results cannot be con- 
veniently obtained using the straightforward numerical ap- 
proach, that have been desc: ‘ibed, because the temperature at any 
point is continually rising though the temperature disturbance 
reaches a maximum in a semi-infinite body. A superposition 
method would be desirable to obtain the steady-state temperature 
errors directly. Since it is known that the void effect extends to 
infinity, some mechanism must be used to simulate infinity. 

The point source solution is, however, valid at large distances 
from the void. The boundary conditions used to simulate infinity 
were the point source temperature results given by Equation (20) 
atr/R = 6 and z/E = 6forg = 2/3, 1, and */,. Forg = 3 and 
also */, the rear boundary was extended to z/E = 12. Asin the 
point source solution, the front surface is insulated. The cylindri- 
cal void is insulated at its radial surfaces and has a uniform heat 
flux over the circular top surface of the void. 

The numerical steady-state hot-spot temperatures are shown 
by Fig. 6. It may be observed from Fig. 6 that the numerical re- 
sults for g = 1.5 with z/E = 6 is about 6 per cent below the value 
given using z/E = 12; the latter result is probably more accurate 
since then the boundary conditions obtained from the point source 
solution would have less effect at the hot spot. As small values of 
g the numerical solution should approch the point source results; 
at larger g values the numerical results approach the large g 
analytical results. 

The steady-state surface temperatures for g = 0.667, 1, 1.5, 
and 3 are shown by solid lines of Fig. 5. The point source solution 
[Equation (6)] is shown to approximate the numerical results 
more closely at [1] larger r/R’s for all g’s and for [2] smaller g’s at 
allr/R’s. The greatest disagreement, as would be expected, is in 
the immediate area above the void. 
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4 Comparison of Results. A comparison of the transient results of 
M and § the authors made in Figs. 2 and 3 indicates differences 
in the magnitude of the effects investigated. The hot-spot tem- 
perature disturbance given by M and § is larger than that given 
here (Fig. 2); the shape of the curves, however, given by Fig. 3 
are similar. 

The apparent lack of close agreement between the M and 8 and 
the present results is attributed to the analysis of the problem. 
M and 8 are concerned only with transient results and do not show 
steady hot-spot temperature. Also, they were interested only in 
the upper bound of the void effect. 

In contrast, solutions presented here give the actual hot-spot 
temperatures and other effects and show the upper and lower 
bounds. In addition, through analysis of the transient and steady- 
state effects, it is shown that a steady-state solution is applicable 
to transient problems since rapid response times are involved. 
Such steady-state solutions have been obtained with a digital 
computer as well as by means of theoretical analytical approaches 
offering means of independent verification of results; the steady- 
state solutions in turn were verified by numerical transient calcu- 
lations. 

Since several solutions presented here demonstrate the validity 
of the authors’ results and since the M and § solution is ad- 
mittedly an upper bound, the results of the present investigation 
are recommended. 

5 Interpretation of Results. Transient results for the hot spot are 
given for g = 0, ?/s, 1, and */; (Fig. 3). The g = 0 values are ob- 
tained from the point source solution while the other were 
obtained by a numerical solution employing a digital computer. 
The authors’ results (Fig. 3) are consistent and recommended to 
obtain the transient hot-spot temperature (7) — 7'..) from y = 0 
to 1.5. The time-variation of the thermocouple tip temperature 
compared to the undisturbed surface temperature (7, — To.) is 
shown for g = */3, 1, and 4/2 by Figs. 9 and 10. 

In comparison with the transient results, steady-state results 
are given which require a decision which is the most accurate. 
A numerical solution would usually be preferred for a given 
problem with irregular geometry since it may approximate the 
boundary conditions most closely. The steady-state numerical 
solution is about 10 per cent less than the large g analytical solu- 
tion near g = 1 and even closer at larger g’s. If the transient re- 
sults for the hot spot (Fig. 3) are extrapolated parallel to the point 
source solution, the steady-state results labeled “Transient 
Numerical Solution’ (Fig. 6) are obtained. These results are 
about 5 per cent below the numerical results. The point source 
solution is on the average about 15 per cent below the numerical 
solutions from g = 0 to 4. The steady-state numerical solution 
appears to be the most accurate; the steady-state numerical 
results are then recommended in Figs. 5, 6, 7, and 8. 

The null point is shown by the numerical solution of Fig. 8 to 
beg = 1. In locating a thermocouple at the null point of g = 1 
several assumptions are made. The first is that the thermal con- 
ductivity of the thermocouple assembly is much less than that of 
the heat sink. If the thermocouple assembly thermal conduc- 
tivity (k,) is of the same magnitude as the heat sink’s conductivity 
(k) then a larger g should be used. A suggested g-correction 
factor, Cre would be 


Cre = bak (9) 

The next consideration concerns the transient variation of the 
null point. The response is usually quite rapid for many thermo- 
couples; there may be cases, however, when the thermocouple 
would be located relatively far from the surface or is in a relatively 
low conductivity material. Both of these effects would lengthen 
the response time and would cause the null point g to be also 
larger (Fig. 9). 
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Fig. 10 Steady-state numerical temperatures at tip of thermocouple 


The last point is explained by the fact that the determination 
of the null point was based only on the point at the tip of the 
thermocouple. Actually the thermocouple junction would inte- 
grate to find the average temperature of the void surface that it 
contacts. The contacting area probably would not be large but 
would again tend to increase the null g value. Fig. 10 shows the 
variatiou of temperature along the top surface of the void. This 
figure indicates that the effect of the thermocouple junction ex- 
tending over the tip area is small until 2/2 when the null g is 1.07. 
In summation, the null point g-value depends upon the actual 
problem though perhaps the best value for most cases would be 
about g = 1.1. 

C Applications. The hot spot and null point have been found 
based on a uniform heat input while the heat sink would, in gen- 
eral, be exposed to a variable-in-time heat input. If the variation 
of the heat input in the response time, Ago, divided by the average 
do at that time is small, then the steady-state results may be 
used. The criteria would be 

At <K (10) 
qo 
where K would be less than some permissible value, probably 
about 0.1. It should be emphasized that the condition expressed 
by Equation (10) is usually satisfied, consequently, the steady- 
state thermocouple results may be used in an otherwise transient 
heat conduction problem. 

To demonstrate that the steady-state solution may be used for 
a thermocouple which is placed in a heat sink exposed to a time- 
variable surface heat flux, consider the following problem. A 1-in. 
thick copper sink is exposed to a sine-squared shape heat input 
[7] for 12.5 see with a maximum flux of 800 Btu/ft*? sec. This 
heat input is similar to that which a re-entry vehicle might ex- 
perience. The rear surface of the sink is insulated. It is desired 
to use a 0.044-in. diameter thermocouple to record the surface 
temperature accurately without causing the surface tempera- 
ture hot-spot temperature to exceed the undisturbed surface 
temperature by 2 per cent. Calculations are needed to determine 
the location of the null point and to insure that the allowable hot- 
spot temperature is not exceeded. Withg = 1.1, E would be EZ = 
R/1.1 = 0.022/1.1 = 0.02 in. From Fig. 6 the steady-state hot 
spot AT’ would be 

1.1(0.6)(0.022/12) | 


= = 6 = 0. 
To. = Tre k (0.6) 220/36 do = 0.0198 go 
(11) 
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for k = 220 Btu/ft-hr-deg F. ¢(Btu/ft? sec) is the surface heat 
flux at the time the hot spot A7' is required. The undisturbed 
surface temperature is shown by Fig. 12 for constant properties 
and zero initial temperature (ac, = 4 ft?/hr). The hot spot A7’s 
computed using Equation (11) are also shown by Fig. 11. 

The hot spot may also be expressed as per cent of the undis- 
turbed surface temperature. This may be misleading as shown 
by Fig. 11 since the maximum per cent increase occurs at early 
times when the temperature level is quite low. Though at very 
early times the per cent temperature increase at the hot spot 
may be as high as 12 per cent (Fig. 2), the per cent is shown by 
Fig. 11 to decrease rapidly and to be about 0.6 per cent at the time 
of maximum surface temperature; the thermocouple hot spot 
then does not exceed the 2 per cent allowable temperature in- 
crease at the time of maximum surface temperature. 

In using Fig. 6, two assumptions concerning transient effects 
were made. Use of the steady-state results assumes a small 
variation of go in the response time. From Equation (8) the re- 
sponse time may be computed for this case to be 0.075 sec. The 
permissible value of K [Equation (10)} would not be exceeded for 
this response time for the given heat flux except at the lowest 
fluxes. This indicates that the steady-state results may be used to 
obtain the hot-spot temperatures for the variable-in-time heat flux 
in this example. The other assumption is that the Z ratio (where Z 
= E/L) is such that the rear surface does not cause the hot spot to 
appreciably exceed its steady-state value (Fig. 4). The value of 
at/E* for the total heating time is 5000 and E/L is 0.02. Since 
from Fig. 4 the Z-value of 0.1 does not exceed the point source 
steady-state value by more than 6 per cent at rg, = 5000 and 
since the Z = 0.01 value does not exceed the steady value until 
Tg = 100,000, it may be safely assumed that the steady-state 
value will be not appreciably exceeded in our example. The hot- 
spot results given by Fig. 11 would then be essentially the same 


using the steady-state solution (Fig. 6) or the transient solutions - 


(Figs. 3 and 4). 

For variable in-time heat inputs the principle of the convolu- 
tion [8] may be used to sum results for uniform heat fluxes to ob- 
tain any desired variable heat flux. The steady-state results 
yields an upper bound for all cases except when the rear surface 
affects the hot-spot temperature (Fig. 4). In this case and also 
when the surface heat fluxes vary rapidly in time, the convolution 
may be applied; it is, however, more cumbersome than using the 
steady-state results. 
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Summary and Conclusions 


The temperature disturbances caused by a thermocouple 
placed normal to the heated surface of a heat sink or calorimeter 
are analyzed. The problem is similar to the determination of the 
temperature errors caused by an insulated cylindrical void near 
the surface since the effective thermal conductivity of the ther- 
mocouple is much lower than a heat sink material. The present 
report gives analytical and machine analyses, which investigate 
a comprehensive range of the important parameters. An im- 
portant contribution of the present investigation is the inclusion 
of a steady-state as well as transient solutions. 

1 Hot-Spot Temp e Disturb Though the temperatures 
are varying with time in the heat sink even for a constant in-time 
surface heat flux, itis demonstrated that the temperature disturb- 
ances caused by insertion of the thermocouple rapidly acquire 
values which are directly proportional to the surface heat flux; 
and when this occurs the disturbances are said to be in “‘steady 
state.’’ Thus the temperature disturbances are usually given for 
transient temperature measurements by a steady-state solution. 
The transient surface hot-spot temperature for constant surface 
heat flux is shown by Fig. 3. For small voids (g < 0.5), the tran- 
sient solution may be given by Equation (3) which indicates, as do 
the numerical solutions which are for larger voids, that the tem- 
perature disturbances near the thermocouple reach 80 per cent of 
the steady-state values when rg = 30 which is called the response 
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TIME (SECONDS) 


Exampie showing the effect of the thermocouple upon the hot-spot temperature 


time. In most transient thermocouple measurements, it will be 
found that the surface heat flux does not vary significantly in the 
response time. 

The steady-state hot-spot temperature disturbance is shown by 
Fig. 6 for various locations of the thermocouple tip relative to the 
heated surface. A convenient equation for the hot-spot tempera- 
ture for small voids (g < 0.5) is Equation (7). 

2 Null Point. Both the transient and steady-state solutions in- 
dicated that a minimum error in the thermocouple measurements 
of the undisturbed temperature would be present if the thermo- 
couple were located at the null point. When the thermocouple 
is located at this poisition, the temperature disturbance at the 
thermocouple due to its presence equals the undisturbed tem- 
perature drop from the sink surface to the thermocouple tip. 
From the consideration of the null point, the optimum place- 
ment of the thermocouple was concluded to be when g = 1.1. 

Effect of Void Distribution. The steady-state surface temperature 
disturbance for one void is given by Fig. 5 as a function of radial 
distance from the void. By properly superimposing the tempera- 
ture disturbances found from this figure the effect of an arbi- 
trary distribution of voids may be determined. For small voids 
the steady-state surface temperature disturbances may be ap- 
proximated by Equation (6). 

The order of magnitude of the hot-spot error depends upon the 
sink material and surface heat flux. For the typical case of a cop- 
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per heat sink which was considered, the maximum temperature in- 
crease at the hot spot was 16 F; this value could have been much 
higher for another material since the temperature increase is in- 
versely proportional to the thermal conductivity. The maximum 
per cent temperature disturbance at the surface hot spot occurred 
at very early time and was about 12 per cent for the geometry of 
most interest; at the time of maximum surface temperature, the 
per cent temperature disturbance was less than 1 per cent. 
The per cent surface disturbance would also increase for a sink 
material of a lower conductivity. 
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APPENDIX 
1 Constant Heat Source [3] 


Consider a constant point source of strength Q(Btu/hr). Inte- 
grating with respect to time, the solution for an instantaneous 
source given by Carslaw and Jaeger [4] is 


T(z, y, 2,0) = f (12) 
where 
b? = (x — + (y — yo)® + (z — 2)? 
The solution [4] of Equation (12) is 


Q 


b 
y, z,t) = - rfc —— 13 
T(2, y, z, t) pare, c (13) 


The insulated surfaces at z = 0 and L may be simulated by 
placing heat sources along the z-axis at +2, +(2L — E), +(2L 
+ EB), =(4L — EB), =(4L + EB), to give 


1 
T(r, z,t) = erfc + erfc 
(14) 
where 
b,? = r? + (2 — E — 2nL)* 
b,? = r? + (z + E — 2nL)? 


For large times the heat flux g, at a point Z in a slab with a 
uniform heat flux go at z = 0 and insulated at z = Lis 
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(15) 


For a void of radius R the point source strength would be— 
L-E 


Q = rR4, = (16) 


Equation (14) then becomes 


T 1-Z |B, 1 |B,| 
- = erfe + 
/k 4 2Vr, |B:l 
(17) 
with 


2 Time-Variable Heat Source 


For a semi-infinite body with a uniform heat flux at its surface, 
the heat flux at any point z = E is 
UP 


— = erfe (18) 
qu 


2 Vr. 

Using the procedure outlined, but with this variable heat flux 
and one image reflection to simulate an insulated plane at z = 0, 
the temperature disturbances for a void in a semi-infinite body 
are: 


opR/k (19) 


with B,* and B,? as given below Equation (17) with n = 0. 
The steady-state temperatures in a semi-infinite body are from 
Equation (19) 


T(r, z) 1 ( 1 1 ) 


1 J. L. Masters and 8. Stein, “Effect of an Axial Cavity on the 
Temperature History of a Surface Heated Slab,” Review of Scientific 
Instruments, vol. 27, 1956, pp. 1065-1069. 

2 B. Bussell, “Analog Study of a Thermocouple Instrumentation 
Error,”’ Dept. of Engineering, University of California, Los Angeles, 
August, 1957. 

3 M.8. Klamkin, “Heat Transfer Due to Source Between Two 
Parallel Insulated Planes,"” AVCO RAD-2-TM-58-128, November, 
1958. 

4 H.S8. Carslaw and J.C. Jaeger, “Conduction of Heat in Sciids,” 
Oxford University Press, London, England, 1947, p. 252. 

5 M.8. Klamkin, “On a Heat Flux Determination in a Cylinder,” 
AVCO RAD-U-58-12, April, 1958, p. 4. 

6 Private Communication, August 9, 1957, M. S. Klamkin, 
AVCO RAD. 

7 Private Communication, R. E. Mascola, AVCO RAD. 

8 H. Hurwicz and M. 8. Klamkin, “Conductive Heat Transfer in 
a Re-entry Body-Analysis of Computational Methods and Their 
Reliability,”’ presented at this symposium. 


Transactions of the ASME 


= 
r \? 2n \? 
Be =(—) + 
4 

; 
j 
References 
a 


G. HORVAY 


General Electric Research Laboratory, 
Schenectady, N. Y. Mem. ASME 


Freezing of a Growing Liquid Column 


The one-dimensional problem is treated in which a liquid metal column, of the original 
temperature Jo, is supplied with additional metal at the top, of temperature 3,, at the 


rate v ft/min. 


At the bottom of the column a heat sink of temperature 0, is main- 


tained. The case in which 3, exceeds the fusion temperature is solved for v = const, 
both rigorously and approximately; the case in which 8, is below the fusion tempera- 
ture so that a freezing front sweeps upward in the column is next treated by an approxi- 
mate method. The foregoing analysis may represent, for instance, the filling stage of 


a casting process. 


The present analysis will be adapted to the more important con- 


tinuous casting process in a succeeding paper. 


| Introduction 


URING THE FILLING STAGES of a casting process 
heat flows, through a growing liquid column and solidified metal, 
to the heat sink at the base. The difficulty of analytically han- 
dling this problem was pointed out by Chao and Weiner [1].?_ In 
their study they therefore restricted themselves to determining 
the temperature J in a semi-infinite solid which is in perfect 
contact with a perfectly mized liquid column (i+., one whose 
temperature is uniform throughout); the column is growing 
linearly with time. (It is assumed that the initial solid tem- 
perature #3, exceeds the freezing temperature of the liquid #,. 
The liquid is poured at a temperature 3, > 3,.) In the Chao and 
Weiner paper the quantity of primary interest is the interface 
temperature of the liquid and solid. 

In the present paper we reconsider the problem of Chao and 
Weiner from a somewhat different point of view. We drop the 
assumption that the fluid is perfectly mixed, but assume, on the 
contrary, that heat is transferred across the fluid by one-dimen- 
sional conduction. On the other hand, we assume in Sections 
II to V that the semi-infinite solid is replaced by a heat sink 
of fixed temperature 3,. In section VI we drop the assumption 
that 3, > d, (thus we allow freezing to take place), and also re- 
place the fixed heat sink by a conducting slab of finite thickness. 
We treat the problem in an approximate manner, by the method 
of lumped parameters. This method, first introduced for this 
class of problems by London and Seban [2], was later generalized 
to problems of continuous casting by Horvay [3a], Giaever and 
Horvay [3b], and Horvay [4]. 


Il Rigorous Determination of the Temperature 
in a Growing Liquid Column 

In Fig. 1 we illustrate the problem on hand. A plate (mold) 
of original temperature 3, and isothermally maintained lower 
face, is brought in contact at the time ¢ = 0, with a liquid column 
of original height Yo and original temperature Jo. Thereafter 
the height of the column is increased at the rate v ft/hr by addi- 
tion of liquid at temperature 3,. We seek to determine the tem- 
perature distribution J in the liquid phase, in the frozen phase, 
and in the plate. The two most important quantities to be 


1 We assume one-dimensional heat flow. The container side walls 


are assumed to be insulated. 

? Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of Tae AMERICAN 
Society or MecnanicaL ENGINEERS and presented at the ASME- 
AIChE Heat Transfer Conference, Storrs, Conn., August 9-12, 1959. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, May 
14, 1959. Paper No. 59—HT-11. 
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determined are: The location X(¢) of the frozen-liquid inter- 
face; and the temperature of the liquid surface. 

If we assume the metal plate as infinitely thin, and the ambient 
temperature as exceeding the fusion temperature, then the mathe- 
matical problem reduces to that of solving the equation 


or dn? 
subject to the initial and boundary conditions* 
r=0: U=U,,H =1 
7 =0: U=0 


_ dH , aH ou 
Us dr 


(1) 


= H: 


ar On 
In the foregoing 


= U(H,r) (8a) 


U(n, r) = @-%), 


are dimensionless temperatures (A = latent heat of fusion,‘ c = 
specific heat, y = specific weight, k = conductivity, x = diffusiv- 
ity = k/ye) 


= (3b) 


is dimensionless time, and 


n=v/Y, H=Y/Y, (3e) 
* For simplicity we treat in this section only the problem in which 
the initial temperature J of the column equals the pouring tempera- 
ture &. In sections III to V we drop this restriction. 
* When no freezing is involved then one may use for \ any refer- 
ence value which has the dirnensions Btu/lb. When there is freezing 
then the use of the heat of fusion for \ is mandatory. 


UNIT AREA 


2 


FROZEN 


: PLATE 4, 


Fig. 1 Freezing of a growing liquid column 
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are dimensionless distance and height of the liquid column, 
respectively. The dimensionless parameter, representing the 


pouring speed 


(H = dH/dr) (4) 


is usually referred to as the “Peclet number.” The boundary 
condition (2c) means: The heat supplied to the liquid column at 
its surface (per unit time) equals the heat gained by the column 
through material growth, plus the heat gained by conduction. 

It is shown in Appendix 1 that the solution of the problem (1), 
(2) is 


U(n, 1)/Us = 1— 2 By(r) sin (5a) 
1 


where 


t,(r) = root of tH cot ({H) + HH = 0 


2 
Han + +H 


A,(r) = exp [- | ¢,(r) exp | dr 
( 


For the case of linear growth of the liquid column, 


B,(r) (5c) 


5d) 
H=1+pr (6) 
(p is the Peclet number) the use of the approximate 4,(r) formula 
10 


Up 


RIGOROUS 
——— APPROXIMATE 


r 10 2 5 10 


Fig. 2(a) 


SIMPLIFIED 


RIGOROUS 


PRIMITIVE 


Fig. 2(b) 


Fig. 2 Determination of (dimensionless) temperature U in the growing 
liquid column versus (di jonless) time r+, for initial value U(n, 0) = 
Uy (pouring temperature) and initial (dimensionless) height H(0) = 1; 
at levels » = H (top of column), 7 = H/2 (mid-height), and 7 = 1 (original 
height); when the bott of the col is maintained isothermally at 


U(O, r) = C, and the Peclet number (dimensionless velocity of pouring) 
is 1 
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(41a) in place of (5d) leads to results which, in a plot like Fig. 
2(a), are indistinguishable from the true solution. In Table 2 of 
the Appendix we tabulate A,(r) for thecasep = 1. In Fig. 2(a) 
we plot, in full lines, the temperature U(n, r)/U, at 9 = H(r) (the 
top of the liquid column), at 7 = H/2 (the mid-height of the 
column), and at 7 = 1 (the original height of the column). In 
dashed lines are shown the results for 7 = H, H/2 by the “ap- 
proximate” method, Equation (11) of Section IIT. These latter 
curves are derived from the lumped circuit representation, 
Fig. 4(a), of the liquid column in Fig. 1; the agreement seems to be 
very satisfactory. In Fig. 3 we give the results by this method 
also for p = 2, 4, 8, */s. 

In Figs. 4(b and c) we use still simpler lumped circuit repre- 
sentations for Fig. 1; we refer to these approximations as the 
“simplified” and “primitive” models. The dotted curve in Fig. 
2(a) gives U(H, r)/U, by Equation (13b), Fig. 4(b). The latter 
approximation method is evidently inferior to that based on 


titiitiiis 


Fig. 3(b) 
Fig. 3 Variation of U(H, r) and U(H/2, r) with r for various Peciet 
numbers. (Approximate method.) 


Fig. 4(c) The “approximate” (2-node) model; (b) the “simplified” 


——— model; (c) the “primitive” (1-node) model for the growing liquid 
mn 


(a) 
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Fig. 4(a), but also simpler to program on a computer, and there- 
fore will be adopted for the more complicated problem dealt with 
in Section VI, the combined growth and freezing of the liquid 
column. The solution by Equation (13b) is compared in Fig. 
2(b) with the still poorer approximation, Equation (14), Fig. 
4(c). 

Figs. 2 and 3 show that the temperature at the top of the liquid 
column begins to sink, from its original value, U,(0) = U,, toa 
minimum near pr = 1 (i.e., near the time at which the original 
column doubles its height; the conduction to the heat sink is most 
strongly felt in this early period, pr < 1); then rises again toward 
U> as the column becomes very high, i.e., the conduction to the 
heat sink becomes vanishingly small. When pouring is done very 
slowly (p< 1) the minimum sinks close to the heat sink tempera- 
ture (U:, min —> 0), when pouring is done very fast (p > 1) then 
the minimum of U, barely drops below U,,. 


Approximate Solution 


We assume that the entire heat absorbing capacity of the 
liquid column ycY (per degree) is concentrated at the mid-height. 
Then we obtain the two-node model illustrated in Fig. 4(a).* 
The lower terminal of the column is maintained at the constant 
temperature #,. Kirchhoff’s law of current balance applied to 
the nodes #; and #, states that 


0= Q Qn 
” = qn + Qu — Ge 
Q = — 8,)dY/de 


(7a) 
(7b) 
where (8a) 


is the heat supplied per unit time, 
d 
[ye — 


is the charging current, 
qn = 0,)/Y, fia => 2k(d, 0,)/Y 
are the gonduction currents, and 


d(ycY) 


(8d) 
is the eurrent due to mass transfer. [A more elaborate explana- 
tion of this term is given after Equation (27).] 
Introducing (8) into (7a) and (7b) and multiplying through 
the equations by cY./Ak yields, on using the notation (3) 
2 dH 


dH 
- 7 U- W- 


9a 
H (9a) 


d(U,H) _ 2 


9b) 
dr H 


Addition of the two equations leads to the single differential 
equation 

dr 


The solution of this equation, for linear growth 


H=1+ pr 


and initial value 


‘From Fig. 3(a): pr < 1 at the minimum for small p, > 1 for 


large p. 
* We follow the method developed by Giaever and Horvay in an 
earlier work [3b]. 
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U,(0) = Us 


pr 


2 
1 + pr 1+ 
(11a) 


From (9a) we then obtain 


(11d) 


The curves U,(r), Ui(r) are plotted in Figs. 3(@ and 6) for the 
initial value 


Up = U, 


and the Peclet numbers 8, 4, 2, 1, */2, '/«, 1/s.7 It may be men- 
tioned that the calculation of Figs. 3(a and b) took but a small 
fraction of the time required for calculation of the rigorous solu- 
tion in Fig. 2 for the single Peclet number, p = 1. 


1¥ Simplified and Primitive Solutions 


If, in contrast to Fig. 4(a), we approximate the liquid column 
by two capacities, one resistance, then we obtain the “‘simplified’’ 
model of Fig. 4(b). Since the lower capacity is maintained at 
the ambient temperature #,, there is now effectively only a single 
node, 3,, involved in the problem. Proceeding as before we find 
the differential equation 

1 d(U,H) dH 1 


dH 


(12) 


(13a) 


The last term represents the heat flow due to mass transfer. 
The solution for linear growth, H = 1 + 77, is 


—2r 
+ 
U(r) —=—_ Es + 1 


(1 + pr)? 


2/p 
1+ pr 


ep | 
pr 2 [ ( 
+ — |B 
(1 + pr)? p* 


1+ pr 


| (136) 


U, is plotted in Figs. 2(a and b) for p = 1, Us = U,. It is seen 
that the one-node approximation is inferior to the approximation 
of Section III, but the results still may be used in a cautiously 
quantitative fashion. 

An inferior 1-node approximation is illustrated in Fig. 4(c).! 
Here all the capacity is concentrated in the condenser at the top. 
The governing differential equation is 


d(U,H) aH 

dr dr 
7 In the figure we write U(H, r), U(A/2, r) in place of Ui(r), Ui(r). 
* This ‘primitive’ approximation offers no advantages in the pres- 


ent problem, but may do so in more complicated problems, see dis- 
cussion of Fig. 11. 


(14a) 
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if 
| (106) 
is . 
—2r | —2r ] 
+ pr + pr 
1+ pr 1+ opr 
—2r 
U(r) = = 
fi 
(8b) 
(8c) 
i 
i 
(6) 


and the corresponding solution is, see Fig. 2(b) 


Multinode Approximation’ 


In Fig. 5 we illustrate the case of the n + 1-node approxima- 
tion, forn = 3. The Kirchhoff’s equations now become 


dH 6 dH 
0= U, or H U)) U; (15a) 
1 d(U,H) 6 4 3 
- H (U, — Ui) H (Ui — U2) 


1 d(U.H) 
2... dH 1 dH 
+ 3 O dr 3 Us dr (15¢) 
1 d(U.H) 6 dH 
—- = —U; om 1 
3 or H (U: U;) H U + 3 U; dr (15d) 


Because of the excellence of the two-node approximation there 
did not seem to be much point to working out a numerical example 
on the basis of Equations (15). Of greater interest is the estab- 
lishment of the limiting form which equations of the type (15) 
assume for the case n 

Equation (15c) may be written as 


H dU; 
3 dr dt 
2 
3 . (3 i) 
= (U, — 2U2+ + U1 (16) 


Denoting furthermore (A is not to be confused with the heat of 
fusion) 


* Similar ideas are developed also by D. R. Otis, Proceedings, 1956, 
Heat Transfer Institute, Stanford, Calif.; and by W. D. Murray and 
F. Landis, Journat or Heat Transrer, Trans. ASME, series C, 
vol. 81, 1959, pp. 106-120, using the Landau equation (18) as starting 
point, and converting it into a system of difference equations. 


Fig. 5 The multinode model for the 
growing liquid column 
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H H = Had, U.-Us=dU 


(17) 
Equation (16) becomes 
Han 2 yay 
i.e., 
au _ wv _ 1 = 
or H ddA dr H* odd? 
subject to the initial and boundary conditions” 
r=0: U = U,(A) (19a) 
A=0: U=0 (196) 
A=1: (19¢) 


where H is a prescribed function of r. Equation (18) was previ- 
ously derived, in a different manner, by H. G. Landau [5]. 

Whereas the finite difference approximations (9) or (15), to 
(18), are very convenient both for numerical and for analog com- 
putations, the partial differential equation (18) is exceedingly 
awkward to use, and therefore no attempt was made to numeri- 
cally verify that solution of the problem (18), (19) coincides 
with that of (1), (2). 


Vi Freezing of a Growing Liquid Column 


We represent the system of Fig. 1 by the 4-node circuit in Fig. 
6. Heat enters at the top of the liquid column (which is of 
length Y and original temperature 3)), then passes through the 
frozen column (of length X) and the mold (of length L), to the 
heat sink at the bottom, which is maintained at the ambient 
temperature 3,. 0, is the temperature at the top of the liquid 
column, 0, is the fusion temperature, 0, the temperature at 
the bottom of the solidified portion, and #,, the temperature at the 
center of the mold. The quantities whose time variation we 
want to determine are, in order of importance, X, 3, 3,, and #,,. 

Evidently, we are representing the liquid and frozen columns 
in Fig. 6 by the less accurate Fig. 4(b), the mold by the more 


Condition (19¢) results from (15a). 


Fig. 6 Simplified circuit model for the 
freezing of a growing liquid column 
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accurate Fig. 4(a). Since our main interest 1s the development of 
the simplest circuit which represents the foregoing four quanti- 
ties to correct orders of magnitude (so that the results may be 
used for semiquantitative predictions in continuous casting 
problems), we refrained from using the more elaborate Fig. 4(a) 
for the liquid and frozen portions."' 

The Kirchhoff circuit laws at the four nodes give 


9: gy — Qn — 


+ Qin + Q,* Q, 


= ~ Yom Qn 


Im = Yam ~ Ima 
In the foregoing 
Q, = — = d(X + Y)/dt) 


(21) 
is the heat supplied per unit time 
= hd, — (22) 


is the radiation loss from the liquid metal surface (for which we 
also allow in our analysis; h, is the heat-transfer coefficient to the 
ambient), 


d yc(b, — 3,)Y d yc(0, — 8,)Y 


dt 2 ’ dt 2 
d  d yeld,- 
dt 2 2 


4 = 


4. = (23) 


d 
In = a — 3,)L) 


are the charging currents (we assume the metal properties to be 
the same in the liquid phase and in the solid phase; conductivity, 
density, and specific heat of the mold material are denoted by 
K, r, C), 


Q; = yAdX /dt 


is the latent heat current produced during solidification, and 


11 That would have added two more nodes to our circuit representa- 
tion, at the centers of the two sections. We state the equations for 
this more elaborate model in Appendix 2. 


(24) 


‘77 


(X+Y) 


<4 


‘ 
ax 


(0) 


Fig. 7 Explanation of the mass 
transfer terms due to freezing and to 
growth of column. See text after 
Equations (27). 
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= k(0,-8)/Y, apn = KO, — 8)/X | 
— d,, — 


dam = » Ima 
1 L/2 
(+ + x) L/2 


(25) 
are the conduction currents from / to f, from f to s, from s to m 


(here we included also a contact resistance term 


(26) 


to represent the airgap resistance between frozen metal and 
mold), and from m toa. Finally, 


1 1 
= 2 — Qn = 2 — 3,)dX /dt | 


(27) 
Q,* = = 10, ~9,)dX/at 


are the mass-transfer terms. 

The latter terms are more clearly illustrated in Fig. 7. The 
column in Fig. 7(a) represents the heat content in the frozen and 
liquid sections X + Y at time ¢; dashed lines separate the heats 
assigned to the upper and lower condensers. In time dt the freez- 
ing front progresses through a distance dX. The corresponding 
heat, Q,*dt shown by cross-hatching in Fig. 7(b), is thus trans- 
ferred from the lower liquid condenser to the upper solid con- 
denser. (The two terms, —Q,* + Q,*,in the 0, equation of (20) are 
self-canceling, because of the coincidence of the two condensers.) 
The lower liquid condenser regains heat Q;,/.dt (it regains height 
dX /2) from the upper liquid condenser which, correspondingly, 
loses this amount of heat. Finally, Q,,dt represents the heat lost 
by the upper half of the frozen column to the lower half, because 
of the rise, by dX /2, of the separation line. Besides the changes 
brought about by the motion of the freezing front, there are addi- 
tional changes caused by the rise of the liquid column surface 
by amount d(X + Y). Clearly the separation line of the two 
halves of the liquid metal must be raised by d(X + Y)/2; thus 
the upper half loses to the lower half also heat Q,,;:dt." 

We now introduce, in addition to (3), also the notations 


==X/Y,, == hY,/k (= Nusselt number) 
+=—)L/2 
K K, 
Pr 


and replace definition (4) of the Peclet number, in accordance 
with (21), by 


ycYo 


p=Yo/=-H+E (28) 


We also replace the radiation loss term (22) by the simpler ex- 
pression 


= AO, (28¢c) 


12 Note that when there is melt-out, the direction of the dX flow 
reverses; heat is transferred from the solid phase, at temperature 
0,, to the liquid phase at the higher temperature 8;; the representative 
circuit element for mass transfer must therefore be a battery; see 
[3b]. 

13 It would be equally appropriate to refer to » as the ‘‘Biot number,” 
since the conductivity of the liquid and solid phases is the same. 
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(Since the superheat is usually quite small, this change represents 
an insignificant error.)'* On multiplying Equations (20) by 
cY./Ak we are thus led to the governing equations for the system 
of Fig. 6 in nondimensional form: 


1 d(U,A) U,-U 
9 = pU, — f 
1 d= 


d(U,Z) 


2 dr = Pr 


1U U,-U U 


dr Pi pz 


/ 


- We shall adopt the system (29) in [4] to problems of con- 
tinuous casting. This will involve replacement of some of the 
ordinary time derivatives in (29) by convective derivatives. 
In the present paper we restrict ourselves to exhibiting the prin- 
ciple features of the solution of Equations (29) for the very 
simple case where the column bottom (node #,) is maintained at a 
temperature 3, < J, (the last two equations of (29) are thereby 
eliminated); the growth rate v is a constant 


H + = = p = const 
(30) 
H=1+pr-z 


‘and the radiation loss (vU,) is negligible. Then Equations (29) 
reduce to 


= p(U,-—U,) (31a) 


(U, — = U, (31b) 


+- 


In Fig. 8 we show solution 2 (r), U,(r)/U, of this equation 
system, obtained on the Reeves Electronic Analog Computer 
by R. Noorda, for the initial conditions 


= (0) =0, U,0) = U, (32a) 

the supply speed 
p=! (32b) 

and the parameter values 
(a) U,=10, U,=0 
(b) U, = 1.70, U, = 1.60 ‘ 
(33) 
(c) U, = 0.20, U, = 0.10 
(d) U, = 1.60, U, = 1.60 


Case (a) is identical with the problem treated earlier by the 


“ If desired, we may incorporate some of the change in a suitably 
redefined h, 
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simplified solution, the growing liquid column, without any 
freezing (8, > W,). The excellent agreement between the 
calculated curve Figs. 2(a and b), and the computer curve Fig. 8(a) 
is noted. Case (b) is representative for aluminum, with 3, — 
0, = 65 deg superheat and 3, = 170 deg ambient (k = 120, 
y = 170, c = 0.26, X = 170, 3, = 1220; lb, ft, hr, Fahrenheit 
units); case (c) corresponds to aluminum with 65 deg superheat, 
and an ambient 65 deg below fusion temperature. In case (d) 
the bath is originally at the fusion temperature. 


From (316) one notes that for 7< 1: 


w is known as the “solidification” parameter [6].“ 

In Fig. 8 curves (=*)'/2U, are also plotted. From their initial 

1 

values + v,) U,, at = 0, the curves are 
seen to rise to a slight maximum around the r value where U;,/U, 
is a minimum (conduction at this 7 is most “‘effective’’, hence the 
temperature gradient across the liquid most conducive to speed- 
ing up the freezing front); thereafter (=*)" drops slowly but con- 
tinuously, because of the continuous replenishment of the super- 
heat of the liquid. When pouring is at the fusion temperature, 
see Fig. 8(d), then the replenishment is eliminated, (=*)' = 
constant, and, so far as the freezing front is concerned, this moves 
as if the liquid column were semi-infinite; the 2w Jf T law holds 
now for all times. * 
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Fig. 8(b») Fig. 8(d) 
Fig. 8 Motion =(r) of the freezing front, and the temperature U;(r) of the liquid top are illustrated for initial conditions =(0) = 0, U,(0) = U 
(= pouring temperature), boundary condition U(O, r) = 0, Peclet number 1; and for the parameter values (a) U; = 1, Uy = 0; (b), Uy = 1.7, Us 
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These cases are representative for aluminum when (a) ambient is above fusion temperature; (b) aluminum is poured at 65 deg superheat, ambient 
is at 170 deg; (c) aluminum is poured at 65 deg superheat, ambient is at 65 deg below fusion temperature; (d) aluminum is poured at fusion tem- 


perature, ambient is at 170 deg. 
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APPENDIX 1 


Rigorous Solution of the Heat Conduction Equation in a 
Growing Liquid Column 


On denoting 


H 9) 
A,r) = f - 2] sin (35) 
0 b 


and taking the sine transform of Equations (1) and (2) there 
results (see e.g., Sneddon [8], pp. 309-312) 

dA, 

+ [,7A, = A,(0) = 0 (36) 


which is solved by 


A,(r) = exp [ - | f, exp | dr 


(37) 
By the inversion theorem there follows 


U = H* + 
1 - — =2 
U, + 52) +H 


i.e., solution (5) of the text. 

Formally (see e.g., Sneddon [8]) Equations (5) are identical 
with the solution of the heat conduction equation subject to the 
“radiative boundary condition; only (5b) is replaced in the 
latter problem by 


acota+C=0 (39) 


where C is the heat-transfer coefficient from the ambient atmos- 
phere to the conducting medium (which is assumed to be of fixed 
height). In the present problem 


C =HH (40a) 


Fig. 9 Wlustration of the error involved in using a mean value 4 in- 
stead of a (r) in formula (37). The deviation of the mean of the two 
curves from zero represents the error. 
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depends on time, hence so do also a, and 
=a,/H (40b) 


For constant ¢, the integrations in Equation (5) lead to an ele- 
mentary formula for A,; the nonconstant case does not seem to 
have been treated heretofore. 

When H varies linearly with time, as in (6) (the case we are 
interested in), then a very good closed-form approximation A,(r) 
may be obtained for A,(r) as follows. 

If it be permitted to visualize a, in the A,(r) expression as a 
time-independent constant &, then 


A,(r) = &, exp | 
0 


&,*/p ] e*dv 
- 
P P+ Prd 


Pp 1 + pr 1+ pr 


+ Ei(—&,*/p) (41a) 


For values of rt where a,*/p(1 + pr) > 10, (41a) may furthermore 
be replaced by the asymptotic formula (see Fliigge [9], p. 116) 


+ pr)? 
A,(r) a, + pr &,? 


p\ pr 
- (1-2) en} 


Is it permitted to treat the a,(7) in the A, integrands as suita- 
ble constants &,? We note from Table 1 (which extends the 
table in Carslaw-Jaeger [10], p. 378, for suitable values of C, to 
k = 12) that a, varies from (k — 1/2)m to km as C = p(l + pr) 
varies from 0 to @. Evidently a is the root of (39) most sensi- 
tive to variation in r. Let us write A,*, A,**, A, for A,(r) in 


. formula (41) when we define &, by its smallest value a,(0), by 


its largest value a,(r), and by its mean value [a,(0) + a,(r)]/2, 
respectively, in the interval 0, 7. In Fig. 9 we plot the curves 
(Ay* — A;)/A; and (A,** — A,)/A), for the case p = 1. The 
error involved in the use of A, instead of the rigorous A;(r) is 
represented by the deviation of the mean of the two curves from 0. 
Evidently, this deviation is less than 0.25 per cent for r < 1, and 


Fig. 10 Approximate circuit model for the freezing of a growing liquid 
column 
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Table 1 
Roots a, of aCota+C=0 


5 6 


7 


0 
-07008 
-07279 
-07698 


1.222 -8726 


15.00} 2. 1.992 1.504 1. 183 


40.00} 4.186 3.937 3.302 2.755 2.333 


100.00} 6.054 6.692 6.241 5.592 4.972 


less than 2.5 per cent for r < 100; for higher & the error of using 
A,(r) instead of A,(r) is completely negligible. Thus it is 
demonstrated that A, may be replaced, for all practical pur- 
poses, by A,. The A,(7) functions (k = 1,..., 12) so deter- 
mined are tabulated for p = 1 in Table 2. Once the A,(r) are 
available, one readily determines B,(r) by (5c). Next the 
series (5a) is calculated, but one finds that the series converges 
very slowly, and the 12-term oupmetion we used is insufficient. 


However, one easily extrapolates pe to 2» by Saltzer’s method 
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.009 


. 429 


.04957 
05068 
05359 
- 05844 
.06813 
.09714 
. 1929 

3804 

-6971 - 4454 


1.735 1. 264 1.156 1.064 


3.972 - 261 2. 984 2.712 2.520 


{11] (a graphical extrapolation is adequate), and one then obtains 
the U(n, r)/U, plots, 7 = H, H/2, 1, shown in Fig. 2(a). 


APPENDIX 2 
Alternate Models for the Freezing Problem 


In Fig. 10 we replace the 4-node representation of Fig. 6 by a 
more elaborate 6-node model. We now label the station at the 


interface of frozen metal and mold by #; rather than by #,. 
The current balance equations are 
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1.00] 2.0288 4.9132 7.9787 11.086 14.208 17.336 20.469 23.604 26.741 29.879 33.017 36.156 2 
1.02 | 2.0354 4.9169 7.9811 11.088 14.209 17.338 20.470 23.605 26.742 29.879 33.018 36.156 a 
1.04| 2.0419 4.9207 7.9835 11.090 14.210 17.339 20.472 23.606 26.743 29.880 33.019 36.157 ae 
1.10 | 2.0610 4.9318 7.9908 11.005 14.214 17.342 20.474 23.609 26.745 29.882 33.021 36.159 
‘ 1.20] 2.0917 4.9502 8.0028 11.104 14.221 17.348 20.479 23.613 26.749 29.886 33.024 36.162 . 
1.40 | 2.1483 4.9862 8.0268 11.122 14.235 17.360 20.489 23.622 26.756 29.892 33.030 36.168 
Ee 2.00 | 2.2889 5.0870 8.0062 11.173 14.276 17.393 20.518 23647 26.779 29.912 33.048 36.183 ve 
4.00} 2.5704 6.3540 8.3029 11.335 14.408 17.503 20.614 23.730 26.853 29.978 33.108 36.239 - 
4 8.00 | 2.8044 5.6669 8.6031 11.599 14.637 17.703 20.805 23.896 27.000 30.111 33.228 36.348 ap 
, 15.00 | 2.9476 5.9080 8.8898 11.896 14.925 17.974 21.068 24.141 27.225 30.319 33.420 36.529 
= 40.00 | 3.0651 6.1311 9.1987 12.269 15.342 18.418 21.454 24.520 27.584 30.649 33.715 36.779 4 
e 100.00 | 3.1105 6.2211 9.3317 12.443 15.554 18.665 21.773 24.884 27.904 31.105 34.216 37.326 iZ 
Table 2 
Ay (t) 
1.00} 0 0 0 0 0 0 0 0 0 
ake. 1.02} .03868  .07731 .09037 .08320 .05848 .02819 
‘i 1.04] .073968 .1256 —.09222 05978 02874 
1.10] .1625 .1918  .09829 .06321 .03040 
1.20} .2723 2280 .1473 .1071  .08302  .06892 .03316 
1.40} .4155 .2665 .1712 .1247 .09777 08032 03867 
2.00] .6550 .1768 .1200 05621 ie 
4.00] 1.122  .6866 .4615  .3454 .2741 .2266 
poe. 
8.00| 1.769 . 4447 .2193 


0=Q-%- q— Qu } 
an — + Qn — Wy 

0=Q,+ qs + — Qn 


Im = Vim Uma 


where 


Q, = — (X + Y)/dt, 
hd; ‘ta = hd, 
Q, = ydX/dt, = dlyc(, — 
qm = — 3,)L)/dt 


Qu = ye(d, — + Y)/dt, 
Qy, = ye(d; — 8,)dX/dt 


Qa = — 0,)dX/dt 


qn = 2k(8, — 3)/Y, = — 
qn = 2k(P, — 


2k (dy 0,)/X, din = 2K (0; d,,)/L, 


dno = 2K(0,, — 9,)/L | 


Nondimensionalizing Equations (42) gives 


9 
0= pu, vU, - H (U, U,;) - pu, 


2 2 
(UH) = - Uy) Ui - Uy) 


+ pU, — BU, 


9 
E+ (U; — Uy) — (U, — Us) 


Ir] | 


+ — BU, 


~ 

| bo 

o 

| 
S 


2 , 
= (Ui U) + EU, 


9 
O=- = (U,-U,) - U, — Un 
= pr 


mU, = — 
pi Pz 


value, then Equations 3,, 32 of (44) reduce to 


+ UV; U;) = 2(U, U2) 


=? (U, — U3)" 


+ 2(U, U2) = 2(U; 


and one finds, similarly to (34) 
for 1: 
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When the interface temperature 3, is maintained at a fized 


(42) 


(43) 


(44) 


(45) 


‘ U; U2 1 
Vr, U; 2+ w? 
= U; “a 
E 1 + U; 


1+ 


For c(8, — 3,)/A = 0.2, 2.0 this gives w = 0.31, 0.86, when 3, (0) 
= J,, and 0.28, 0.80 when c[v,(0) — 3,]/A = 0.2; in good agree- 
ment with the rigorous solution quoted in footnote 14, when 
(8, — 8,)/(8, — is small, and — is not too large.” 

Nevertheless, Equations (44) have some unpleasant short- 
comings, in common with Equations (29). This is revealed most 
clearly by Equation 3; where, at time r = 0, a current in the 
infinitesimal element =/2 is equated to a current in the finite 
element p;. (It is always desirable in circuit representations to 
employ sections of, roughly speaking, equal resistances. This, 
however, is patently impossible when one of the sections has zero 
thickness.) In consequence, our approximate analysis implies 
U;(0) = U;(0), followed by a gradual, though rapid, separation 
of U(r) from U,(r), whereas in an exact solution initial values 
U0) # U,(O) are assumed, abruptly. Our approximate 
method therefore cannot furnish clearcut estimates for these 
initial values, which are frequently of great interest.“ A more 
annoying disadvantage of the method is that expressions of type 
0/0, 0: are very awkward to program on an analog com- 
puter, and what is gained by use of a more refined model may be 
more than lost through the inaccuracy generated by computer 
elements which inadequately reproduce 0/0 or 0: type quanti- 
ties. 

We may remedy the situation (at the expense of replacing the 
infinite rate = at r = 0 by a finite rate) by lumping all the ca- 
pacity of the frozen section X, Fig. 10, into a single node at 0,. 
(Compare Fig. 4(c), where the lumping is done into the node at 
v,.) This is a particularly good approximation during the early 
times of the freezing process, where the thickness = is small.” 
Fig. 6 and Equations (29) are then modified to Fig. 11 and 
Equations (47): 


- U, 


Bs 


(47) 
U,-U. ; 


where we used the abbreviation 


A=ZE+p (A(O) = pr) (48) 
Fig. 10 and Equations (44) may be similarly simplified to Fig. 12 


and Equations (49): 


" Similar approximate expressions for w* were previously ob- 
tained by Goodman [12]. 

18 Whatever estimate is obtained depends quite sensitively on the 
thickness of the finite element adjoining the infinitesimal element 
whereas, in reality, only the boundary conditions should matter. 

% Many casting processes (coating processes) are never carried 
beyond this stage. 
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pU, = pU, — vU; — H 


2 
HU; + pu; = pu, + H (U, — Ui) 


Use of these equations will be made in [4]. 

Nore: Because of the importance of formula (46) we also 
state it for the case in which the liquid metal properties k,, c, 
differ from the solid metal properties k,, c,.” Writing 

one finds, by steps entirely analogous to (42) to (46), that at the 
incipience of freezing 


%” When 7: * 7, then convection currents are set up, and the pres- 
ent analysis is not applicable without major modifications. See 
Chambré [13]. 


Fig. 11 Alternate model for 
the freezing of a growing 
liquid column 


0, — 2+ w? 
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| 


(51) 
It is noted that w is (in our approximation) independent of k,/k, 


w= 
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Fig. 12 Alternate model for the freez- 
ing of a growing liquid column : at 
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ROBERT G. VINES? 


Measurement of the Thermal Conductivities 
of Gases at High Temperatures 


Experimental results are reported for the thermal conductivities of air, argon, nitrogen, 
and carbon dioxide at temperatures up to 900 C, and of steam up to 560 C. These re- 
sults are compared with values predicted from correlation formulas based on low tem- 


 * DETERMINATION of thermal conductivities at 
high temperatures is unusually difficult because the large heat 
losses produce thermal instability in experimental equipment. 
Investigations previously carried out in this laboratory have 
shown that measurements above 500 C are often liable to errors 
of up to 20 percent. In the work reported here a system of high 
thermal capacity was employed, and consistent and reproducible 
results were obtained. It thus seems likely that the conductivity 
values given are substantially correct, and that reasonably ac- 
curate high temperature measurements have been achieved. 


Apparatus 


At elevated temperatures a concentric cylinder apparatus is 
necessary to reduce the effects of “temperature jump.’’ The 
conductivity cell was of novel design and was constructed en- 
tirely of platinum (see Fig. 1). The emitter (length ~5 cm) was 
a thick cylinder of 6 mm outside diameter (internal diameter 
<2mm). A heater of fine platinum-rhodium wire, wound on a 
porcelain mandril (diam 1 mm) was inserted into the bore of the 
emitter and cemented in position with alundum cement. The 
cement was prepared with distilled water so that, at high tem- 
peratures,.it still remained a good electrical insulator. Double 
leads (current and potential) were gold soldered to each end of 
the heater winding. 

Two cylindrical heat stations (length 2.5 cm, other dimensions as 
above) were connected to the emitter by small thin-walled tubes, 
being separated from it by a distance of <1 mm. The heaters 
were “double wound’”’ on hollow porcelain mandrils through 
which the leads to the emitter were taken. The leads from both 
the emitter and heat stations were then passed through porcelain 
insulators, and these were enclosed in short platinum tubes 
rigidly fixed to the ends of each heat station. 

The receiver was a hollow cylinder 12.5 em long and of 10 mm 
internal diameter (wall thickness ~1 mm). It was separated 
from the inner unit comprising the emitter and heat stations by 
two circular spacing pieces of fired “lava.” 

The completed conductivity cell was mounted vertically in a 
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perature measurements. 


quartz tube directly connected to a glass vacuum system and gas 
supply. Eight thermocouples (Pt against Pt/10 per cent Rh)* 
were set in small flat pieces of gold and smoothly welded to the 


* Two sample thermocouples identical with those used were cali- 
brated at the melting point of silver (960.8 C). During solidification 
of the silver, the temperature, as measured by the thermocouples, 
remained constant (to within better than 1/25 deg C) over a period 
of ten minutes. At this temperature the emf readings of the two 
thermocouples were 9115 u volt and 9116 uw volt in good agreement 
with the value given in the U.S. National Bureau Tables (9120 4 volt). 
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cell at the positions shown in Fig. 1. External connections to 
the thermocouple (and heater) leads were made through Teflon 
seals in the glassware. 

The quartz tube was itself completely surrounded by a thin, 
closely fitting nickel sheath (to help even out any temperature 
discontinuities) and placed in a constant temperature furnace. 
This was made from a porcelain cylinder (diam ~1 in. and some- 
what less than one ft in length) wound with three sections of 
platinum ribbon—a long middle section with adjacent short 
windings to supply extra heat at each end. By manual adjust- 
ment of the heating currents (a-c with voltage regulation), it was 
possible to achieve uniformity of temperature along a considera- 
ble portion of the porcelain tube. The central winding was also 
connected to a stabilizing d-c circuit, the current in the circuit 
being regulated by a platinum resistance thermometer adjacent 
to the furnace; this provided thermostatic control when used in 
conjunction with a conventional resistance bridge, photocell 
circuit, and thyratron control unit. 

The entire furnace was finally enclosed in a double-walled 
cylindrical brass case, filled inside with magnesium oxide and 
asbestos insulation. To provide a constant surrounding tem- 
perature of 100 C the outer compartment contained water which 
boiled when the furnace was at 500 C or higher. With careful 
adjustment of the electrical system the furnace could be main- 
tained at any temperature between 200 and 1000 C for periods 
of an hour or more; over shorter intervals of time the maximum 
temperature variation was no more than a few thousandths of a 


degree. 
Calibration and Experimental Procedure 


1 The thermocouples were annealed by heating “‘in situ’’ at 
750 C for 36 hr in an atmosphere of argon. 

2 In making measurements the furnace was first adjusted 
until the readings of all three thermocouples on the receiver agreed 
to within 1/30.deg C. Current (d-c) was then passed through 
the emitter and heat station windings. To minimize end effects 
and obtain reproducible (and accurate) results, it was necessary 
to keep the temperature differences between emitter and heat 
stations smaller than 1/20 deg C. This produced a nonuniform 
temperature distribution along the emitter, the center becoming 
slightly hotter than the ends. 

3 The heat introduced to the emitter was determined from 
measurements of the potential difference across the winding when 
a known current was passed through it; the product of these two 
quantities gave the total heat input. 

4 Conductivities were calculated from the measured dissipa- 
tion of heat as determined with the apparatus under equilibrium 
conditions. Heat was lost by both conduction and radiation, 
but measurement of the corresponding heat dissipation when 
the apparatus was highly evacuated enabled corrections to be 
made for end conduction, radiation, and so forth; during these 
vacuum runs the temperatures of the emitter, receiver, and heat 
stations were maintained at roughly the previous levels. From 
a smoothed curve relating the radiation losses with temperature, 
appropriate corrections could be made to each experimental 
measurement. 

5 The temperature difference (AT) between emitter and re- 
ceiver was usually between 5 and 10 C. The greatest accuracy 
was obtained when AT’ was large. Within experimental! error, 
the results of the measurements at any temperature were inde- 
pendent of AT. The “experimental temperature’’ was assumed 
to be the mean of the emitter and receiver temperatures. 

6 The apparatus could not be operated at pressures higher 
than atmospheric. Convection effects appeared to be very small 
for, within experimental error, all results were independent of 
pressure over a wide pressure range. In general, measurements 
with gases were made at pressures close to atmospheric, though 
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when water vapor was introduced the pressure was very much 
lower; since external portions of the apparatus were not heated, 
the pressure was only that of water vapor at room temperature 
(~2 cm Hg). 

7 The cell constant as calculated from the simple formula 


const = 24L/log r2/r 


does not take into account heat losses from the end surfaces of the 
emitter, edge effects, nonradial flow, andsoforth. (Lis thelength 
of the emitter, r,; ites radius, and rz the radius of the receiver.) 
As the constant so calculated is clearly somewhat low, it was ar- 
bitrarily increased by a few per cent so that the conductivities of 
air, argon, carbon dioxide, and nitrogen at ~250 C were in 
agreement with accepted values at this temperature. This ad- 
justed constant was then used to calculate conductivities up 
to 900 C, small corrections merely being made to compensate for 
the expansion of platinum with temperature. 

8 Calculation showed “temperature jump”’ effects to be unim- 
portant with the present apparatus, except possibly for the 
steam measurements at 560 C; even in this case the error involved 
(~1 per cent) was less than the likely experimental error. Meth- 
ods of estimating such effects are outlined in references [8]* and 
{12}. 

9 In calculating conductivity coefficients the temperature 
difference between emitter and receiver must be accurately 
known. This requires knowledge of the variation of thermo- 
couple emf with temperature. The factor dE/dT was obtained 
by differentiating the equations relating E and 7’, as given by the 
U. 8. National Bureau of Standards; although this is not entirely 
satisfactory, no better method is available. 


Results 


The results are shown in Table 1. The values listed at 900 C 
are probably less accurate than the other figures, not only be- 
cause of the magnitude of the radiation correction but also be- 
cause the apparatus was damaged at this high temperature. It 
is interesting to compare the conductivity values in Table 1 with 
those expected from extrapolation of low temperature data. 
Empirical correlation formulas have been developed previously 
[1] which reproduce experimental results obtained over a wide 
temperature range (—200 to 350 C). It may be seen that 
“high temperature conductivities’ as calculated from such 
formulas are in excellent agreement with the present results for 


« Numbers in brackets designate References at end of paper. 


Table 1 Thermal conductivities of selected gases 
K-105- cal sec™'- deg 


Temp C 240 260 270 # 550 560 760° 900° 
Air 9.76 16.7 18.4 
Nitrogen 9.92 13.66 16.1 17.8 
Carbon 

dioxide 8.08 13.900 ... 17.5 19.4 
CO./N: mixtures: 760C 

% COz 0 33.5 47.1 65.9; 100 

K-10° 16.1, 17.5 17.4 17.5 17.1° 


* At high temperatures the “blank” corrections (radiation, etc.) 
are large, being of the same magnitude as the conductivity values. 
As might have been expected, the radiation corrections appeared to 
vary roughly with the cube of the absolute temperature; with argon 
at 760 C approximately half the heat dissipated was accounted for in 
this way. 

* This value for pure carbon dioxide is not quite in agreement with 
that given in the upper part of the table. Damage to the apparatus 
at 900 C, which took place before the measurements on mixtures were 
made, is the probable cause of the discrepancy. 
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simple gases like argon, nitrogen, or air, but not with those for 
carbon dioxide. 
Argon. The correlation equation‘ is 


0.388 

186.7 

10-”/T 


K-105 = 


where 7’ is in degrees absolute. Comparison with the earlier ex- 
perimental results of Kannuluik and Carman [2] is shown in 
Fig. 2. There is close agreement between the predicted values 
at high temperatures and the present measurements. Reference 
to the literature indicates that there is excellent agreement with 
the results of Rothman and Bromley [9] at 600 C, and with those 
of Gardiner and Schiifer [13] at 200 and 300 C. Values listed 
by the National Bureau of Standards [14] agree well below about 
600 C, but at 850 C are ~2'/; per cent lower than the curve in 
Fig. 2 suggests. 

For argon, reliable viscosity data are available. If these are 
used to calculate the ratio between conductivity and viscosity, 
a systematic decrease with increasing temperature is apparent 
in the value of the Chapman-Enskog transport factor—from 
2.5 at 0 C to ~2.45 at 1000 C. Since the estimated experimental 
error is at least +1 per cent this variation is scarcely significant; 
however, as reported earlier [3] helium and neon exhibit similar 
behavior which suggests the effect may be real. If so, it raises 
an interesting problem in kinetic theory which does not predict 
a variation of this kind. 

Nitrogen. The results are shown in Fig. 3, together with those 
of Vargaftig and Timrot [4]. The correlation equation is 


0.604 ./T 

224 

— 
1+ = 0-"*/T 


K-16 = 


It may be seen that the agreement between observed and pre- 
dicted conductivity values is excellent. 

Measurements reported by Rothman and Bromley [9] are 
uniformly ~4 per cent lower than the present results from 300 C 
to 800 C, and are not consistent with probable values for the con- 
ductivity of air—cf. the measurements in Figs. 3 and 4 which 
are reasonably self-consistent over a wide range. The results 
of Stops [15] are ~3'/; per cent lower at 900 C, though the dis- 
crepancies are smaller at 500 C. 

Air. The results are shown in Fig. 4, together with those 
of Kannuluik and Carman [5]. The correlation equation is 


__ 0.632 VT 


K-105 = 
1+ -10-"/T 


Again, there is good agreement between observed and predicted 
values. 

For air there are many measurements up to 500 C, and most of 
these are in accord with the present results. Values in the NBS 
Tables are in substantial agreement up to 750 C, but then be- 
come slightly lower; the curve estimated from the Lennard- 
Jones 6-12 potential, using accepted constants, also agrees to 
-~700 C, but is ~2 per cent higher at 900 C. Figures derived 


* The constants used in the correlation equations are taken from 
reference [1]. 
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Fig. 5 Conductivity measurements for stecm 


from the formula of Stops [15] are slightly lower at 500 C, but 
2'/s per cent lower at 900 C. 

Steam. No satisfactory correlation equation has been formu- 
lated since measurement of the thermal conductivity of steam is 
very difficult; this is emphasized by the general lack of agreement 
between values obtained by different workers. Nevertheless, as 
may be seen from Fig. 5, the present results compare well with 
those of other recent investigations [4], [6], and [7].6 Over the 
entire temperature range the curve in Fig. 5 is higher than values 
listed in the NBS Tables, being as much as 7 per cent higher at 
450 C. 


Carbon Dioxide. The results for carbon dioxide are interesting, ~ 


for the values predicted from low temperature measurements 
(ef. also Vines [8] ) are much in excess of those actually observed; 
the latter are, however, in excellent agreement with the recent 
results of Vargaftig and Timrot [4] (see Fig. 6). It is quite clear 
that the extrapolated curve in Fig. 6 is incorrect, for other experi- 
mental values of Stops [15], and Rothman and Bromley [9], 
which are reasonably consistent in themselves, are ~2 per cent 
lower than the present measurements at 400 C, and 5 per 
cent lower at 900 C. Again, the full curve in Fig. 6 is in ex- 
cellent agreement with the NBS Tables—though these extend 
to only ~350 C. 

The correlation equation for carbon dioxide (see reference [1] ) 
was formulated from conductivity measurements in the range 
—50 to 300 C. Over this temperature range the heat capacity 
increases rapidly for, with carbon dioxide, the vibrational degrees 
of freedom become excited at quite low temperatures. Never- 
theless, the heat capacity rises much less rapidly beyond 300 C, 
and since no allowance is made for this in the correlation equa- 
tion, it is not surprising that the predicted high temperature con- 
ductivities are too large. With simpler gases the rate of increase 
in heat capacity is less marked at low temperatures; correlation 
equations should thus remain effective over a much greater 
temperature range—cf. the results given for argon, nitrogen, and 
air. 

Measurements on carbon dioxide/nitrogen mixtures at 760 C 
are shown in Fig. 7. This illustrates the complexities likely to 
be encountered in conductivity measurements at high tempera- 
tures for, as reported earlier by Keyes [1] and Rothman and 
Bromley [9], the conductivity versus composition curve exhibits 
a well-defined maximum. At considerably lower temperatures 


* These measurements of Keyes and Sandell for steam and nitrogen 
are known to be low at higher temperatures. If the results are ad- 
justed so that there is agreement with accepted conductivity values 
for nitrogen over the whole range, the corresponding steam meas- 
urements are then in accord with the values for steam reported 
here (see Fig. 5). 
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Fig. 7 Conducttvity measurements for carbon dioxide/nitrogen mixtures 
at 760C 


the corresponding curves are roughly linear, the slight negative 
deviations observed near room temperature being in good agree- 
ment with the values predicted from the empirical equation of 
Lindsay and Bromley [10]—cf. Bennett and Vines [11]. It has 
been shown [11] that the maxima often observed in the conduc- 
tivity versus composition curves for gaseous mixtures are due to 
differences in the polarity of the mixture constituents. On this 
basis, Fig. 7 suggests a definite lack of symmetry in the charge 
distribution within the carbon dioxide molecule at 760 C, in full 
accord with heat capacity measurements which indicate that, at 
this temperature, both the deformation and valence vibrations in 
carbon dioxide are activated. 


Conclusions 


This investigation shows that the thermal conductivities of 
simple gases may be calculated accurately from low temperature 
measurements by the use of correlation equations (see also ref- 
erence [16]). Even at very high temperatures up to 2000 C, val- 
ues so estimated are unlikely to exhibit gross errors, as is clear 
from a study of Figs. 2,3, and 4. This is of special importance 
to the engineer who can, with confidence, use extrapolated con- 
ductivity values in heat-transfer calculations for gases such as ar- 
gon, neon, hydrogen, nitrogen, oxygen, and so on. 

On the other hand, the method is not suitable for rather more 
complex gases possessing higher specific heats; here the use of 
correlation equations, based solely on low temperature measure- 
ments, could lead to enormous errors in the calculated high tem- 
perature values. Again in the case of water vapor, for which 
heat-transfer calculations are specially important, correlation 
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equations cannot be used; this is particularly unfortunate, as 
even experimental conductivity values are often unreliable owing 
to the difficulties of steam measurements. Existing knowledge 
of steam properties is insufficient for the extreme requirements 
of modern power stations and atomic plants, and The American 
Society of Mechanical Engineers has recently instituted an ex- 
tensive research program to extend the present data. The 
measurements reported here for the thermal conductivity of 
steam were actually part of the investigations carried out during 
the initial phases of this program. 


Acknowledgment 


The author is grateful to Prof. F. G. Keyes for the opportunity 
to work in his laboratory, and for his continual interest and help 
while the work was in progress. Thanks are due to the U. 8. 
Educational Foundation in Australia for the award of a Fulbright 
Travel Grant. 


Comments by Frederick G. Keyes 


The attempt to actively develop high temperature equipment 
for the measurement of the thermal conductivity of gases under 
Project Squid began in January, 1952, with Dr. William T. 
Lindsay, Jr., collaborating in the design and construction of the 
equipment. The thermometry depended for the greater part on 
the use of platinum and platinum rhodium thermocouples and 
much effort was expended in creating special facilities for thermo- 
couple calibrations and thermometry controls. 

The development of a furnace design for temperatures over 
the range 300 to 900 C included provision for temperature regu- 
lation. The necessary requirement of constancy of temperature 
in the tube furnace and mid-portion for a distance of several 
inches proved difficult of achievement. The cell employed was 
of platinum and its functioning suggested numerous design 
changes, in both cell and furnace, which were incorporated in the 
equipment during 1953. The cell was entirely rebuilt and im- 
provements made in temperature control apparatus. The 
thermocouple wire was increased in size from 0.1 mm to 0.2 mm 
because of inhomogeneities which developed in the finer wire. 
The alteration of the thermoelectric properties of wire caused by 
cold working was also reduced by the change to a larger size. 
A large variety of detail changes and improvements were in- 
corporated in the second cell which underwent tests in June, 1953. 
In spite, however, of the care and attention to detail in construc- 
tion of the second cell it was not altogether satisfying. A series 
of observations were made in vacuum and on air, argon, helium, 
and nitrogen. The highest temperature employed was 700 C. 
Different temperature distributions led to variations in results 
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for air, nitrogen, and argon; in each case amounting to as much 
as ten per cent. 

In April, 1954, Mr. Robert G. Vines on leave from the Univer- 
sity of Melbourne continued the work and began the construc- 
tion of the third ceil based on the experience accumulated by Dr. 
W. T. Lindsay. The record of the performance of the third cell 
is given in the present paper. 
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Laminar Transfer From Isothermal 
Spanwise Strips on a Flat Plate 


The Rubesin-Klein-Tribus method is used to calculate rates of heat transfer by 


forced convection from two isothermal spanwise strips in tandem. The general results 
are applicable to either laminar or turbulent flow. Predictions of laminar transfer are 
based on Eckert's step function. 


In experimentation, 
while the inert material between them simulates the adiabatic wall. 


naphthalene cast in trays is used to simulate the isothermal strips 
The velocities 


range from about 40 to 90 ft/sec, and the air temperature from 72 to 88 F. The boundary 
layer is shown to be laminar. 

The theoretical and experimental results are found to be in good agreement, confirming 
the calculation method and Eckert’s step function. 

The experimentation is extended to an array of several equally pene strips in tandem. 


1 Introduction 


= TO DEAL with surfaces of arbitrary tem- 
perature distribution have been developed by Chapman and 
Rubesin [1] ,? Lighthill [2], Rubesin [3], and others. An excellent 
review has been provided by Klein and Tribus [4] who have 
shown, also, how to calculate wall temperature distributions from 
prescribed distributions of heat flux. Of all the methods, the 
Rubesin-Klein-Tribus (RKT) probably is the most convenient in 
practice, and further simplification for engineering usage has 
recently been provided by Eckert, et al. [5, 6]. 

Eckert’s [7] step function, which the RKT method uses for the 
laminar boundary layer has been confirmed by comparison of 
several theoretical solutions [5]. However, until now it has been 
only indirectly confirmed by experimentation [8, 9]. 

In the present investigation, the RKT method is employed to 


1 Numbers in brackets designate references at end of paper. 
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Nomenclature 


Effects of roughness and spanwise diffusion are noted. 


predict the’ performance of isothermal spanwise strips on a flat 
plate, and the results are checked against experimentation based 
on the heat-mass analog. 

With reference to the flat plate sketched in Fig. 1(a), isothermal 
heating occurs from? z,’ to z, and from 2,’ to z, as shown by the 


temperature distribution sketched in Fig. 1(b); these regions 
are the isothermal strips. Elsewhere the surface is adiabatic, 
as shown in Fig. 1(c), which is a sketch of the local rate of heat 
transfer q"(z). Returning to Fig. 1(b), we see that the tempera- 
ture of the fluid at the surfaces of the adiabatic elements rapidly 
decays downstream from each strip, as is shown later in the cal- 
culations. The areas under the curves AB and CD in Fig. 1(c) 
are proportional to the total rates of heat transfer from Strips 1 
and 2, respectively. The rate of heat transfer at the leading 
edge of each strip is very high. One of the questions to be answered 
is whether or not the large rates of heat transfer at the leading 
edge of Strip 2 and of subsequent strips would possibly compen- 
sate for the loss of heated area between them. 


? The prime (’) used with a symbol related to a length refers to the 
leading edge of a strip. The same symbol without the prime refers 
to the trailing edge of the strip. The prime used with g or m (rates of 
heat or mass transfer) refers toa unit span. Similarly, a double prime 
(”) refers to a unit area. 


A = area Pr = Prandtl number, v/a 6 = boundary layer thickness 
B,, B: = constants p = pressure 6 = temperature excess 
b = coefficient of mass transfer q = rate of heat transfer A = defined by Equation (25) 
C = dimensionless coefficient R = gas constant vy = kinematic viscosity 
c, = specific heat at constant pres- Se = Schmidt number, v/D € = variable distance from leading 
sure s = z* edge when z. is considered 
D = coefficient of diffusion T = absolute temperature constant 
F = ratio of total rates of transfer t = temperature p = density; concentration 
h = coefficient of heat transfer by U = velocity at outer edge of bound- o = 
convection ary layer T = time 
I, = incomplete beta function u = velocity component in z-direc- xX = defined by Equation (20) 
K = dimensional coefficient,Ck(U/v)" tion Subscript 
k = thermal conductivity w = spanwise width of strip : : : 
j = index of strip 
1 = chordwise length of a strip z= distance from leading edge n = total number of strips 
m = mass transfer y = distance from surface v = vapor 
mh = rate of mass transfer @ = exponent; thermal diffusivity w = wall 
n = exponent 8 = exponent 0 = mainstream 
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Fig. 2 Simple isothermal step 


In the experimentation, trays of naphthalene are assembled on 
a flat plate, and rates of mass transfer are measured by weighing 
the trays before and after exposure to the air stream. Since the 
strips of naphthalene are virtually isothermal, their surface con- 
centrations are very nearly uniform, and they simulate isothermal 
heated strips [10]. The advantage of employing the analog in 
this case is that thermal isolation can be completely simulated 
by inert surfaces and that a large variety of surface distributions, 
or “‘mosaics,’’ can be tested without any modification of the in- 
strumentation. 

In Sections 2 through 4, which deal with the analysis and solu- 
tion of the problem, a heated plate is considered, and the term 
“strip’’ is used for the longer term “spanwise isothermal strip.” 
In Sections 5 through 8, which deal with the experimentation, 
the term strip refers to the analogical equivalent of the heated 
strip, namely, the naphthalene in the trays. 


2 Analytical Statement of the Problem 


It is assumed that the flow is steady, two-dimensional, and in- 
compressible and that the temperature differences are suf- 
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ficiently small to permit the use of constant values for the fluid 
properties. Since the flat plate is being studied, the pressure on 
the surface and throughout the boundary layer is uniform. Fur- 
ther, the velocity of the fluid is low enough so that effects of fric- 
tional dissipation may be neglected. 

Under these assumptions we obtain the well-known system of 
equations which express the conservation of mass, momentum, 
and energy in the boundary layer. The only one that need con- 
cern us here is the energy equation 


(1) 
“> 


For under the assumptions given, the equations are coupled uni- 
laterally, the functions u = u(z, y) and v = v(z, y) being solutions 
only of the momentum and continuity equations. These func- 
tions may be tabulated, for instance, after the general Blasius 
solution. 

The boundary conditions that ¢ = ¢(z, y) must satisfy have al- 
ready been discussed in conjunction with Fig. 1. In regarding the 
distribution of the flux, it should be observed that q°(z) = 
—k(dt/2y)y-e. Hence, in the first, third, and fifth intervals of z 
in Fig. 1, (d¢/dy) 0 = 0. 


3 Scheme of the Calculations 


The temperature distribution of a simple isothermal step is 
shown in Fig. 2. If g’(z) denotes the rate of heat transfer per 
unit area at distance z from the leading edge, then a local co- 
efficient of heat transfer A(£, x) may be defined by the equation 


q"(z) = h(E, z)-A(z) (2) 
where 
0 for 
Kz) = (3) 
6, for z>& 


and 6,, denotes the constant difference t, — t. 
The relation between the usual isothermal coefficient A(0, z) 
and the coefficient h(£, z) is taken to be of the form 


Typically, for either laminar or turbulent transfer on a flat plate, 
h(0, z) = (5) 
z 


and 


The coefficient C is dimensionless; it depends upon the Prandtl 
number and the nature of the flow in the boundary layer. Also, 
the exponents n, a, and 8 are constants depending upon the 
nature of the flow. 

Hence, for the simple step under particular flow conditions of a 
given fluid, 


0 (z< &) 
q"(z) = (7) 
— pap (z > &) 


The coefficient K has the units of kz~*. 

The significance of Rubesin’s work [3] is that for an arbitrary 
distribution of the excess temperature at the wall, 0(z), the local 
rate of heat transfer can be computed by superposing the solu- 
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tions of infinitesimal as well as finite temperature steps. 
formula has been concisely expressed as follows: 


1 
q'(z) =K — (8) 
the integration being performed in the Stieltjes sense. 
The three exponents satisfy the relation 
n=(1 — B)-a (9) 


Some sets of n, a, and B are tabulated here: 


Nature of the 


boundary layer n a 8 Reference 
laminar 0.5 3/4 1/3 Eckert [7 
turbulent 0.8 39/40 7/39 Rubesin 
turbulent 0.8 6/10 1/9 Seban [11] 


The values for the turbulent boundary layer, which are based on 
semiempirical analyses, are quite well substantiated by ex- 
perimentation. Some fairly recent comparisons have been made 
by Jakob in his communication regarding the work of Furber [12] 
and by Eckert, et al. [5]. 

As already mentioned, Klein and Tribus [4] provide a method 
to calculate the distribution of surface temperatures from a 
known distribution of the heat transfer q’(z). Without repeating 
the details of the derivation, which may be found in [4], we may 
write on the basis of the power relations given that 


asinBe (* 
The following minor transformation is introduced for con- 
3 venience: 
s=2z* and o = (11) 
Then (8) and (10) become 
q'(s) = K (12) 
and 
sin Bx (* 
= f (13) 


4 Application to the Isothermal Heated Strips 


Equations (12) and (13) are used alternately as follows: First, 
(12) is employed to calculate the heat transfer from Strip 1. 
Second, (13) is used to calculate the temperature distribution at 
the adiabatic wall between the two strips. And third, (12) is 
again employed to calculate the heat transfer from Strip 2. It 
is assumed that the boundary layer is either laminar or turbulent. 

Since the total instead of the local transfer is measured, an ex- 
pression for the total rate of heat transfer from Strip 7 per unit 
span is needed. It is obtained from the following relationships: 


(14) 


The results of this calculation will be compared with the ex- 
perimental results. 

4.1 Heat Transfer From Strip 1. For Strip 1, (12) reduces to(7). 
Hence, 


= 


3 In referring to [4] it may be useful to observe that I'(8)-I'(1 — 8) 
= x/sin Br. 
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After replacing s with ¢, sudstituting into (14), integrating, and 
eliminating @ through (9), 


K-6 
= (8, — (16) 


In the particular case that heating begins from the leading edge 


n n 


qo’ (17) 
Hence 


(18) 


4.2 Distribution of the Wall Temperature Downstream of Strip 1. 
Upon substitution into (13) for q"(o) in the interval co = Otog = 


sin Ba do 
Introducing 
x = x(s, 0) = (20) 
we obtain 
sin Br x(s, #1) dx 


or, observing that the integral is the incomplete beta function, 
Hs) = Oui a)(1 B, B) (8 s< 82") (22) 


For the laminar boundary layer 8 = '/;. The normalized in- 
complete beta function /,(*/s, '/;) has been calculated, and the 
results are presented in Table 1. It is evident that the wall tem- 
perature falls very abruptly downstream of Strip 1. Neverthe- 
less, the theory continues on the assumption that the diffusion in 
the z-direction is zero everywhere. 

4.3 Heat Transfer From Strip 2. Upon substitution into (12) for 
Hc) in the interval from o = Oto a = &, 


1 


(s) 6,Ks [s 


d 


‘a 
+f x 4 
[s — 


= — B, B) 
[s — 


The integral is readily evaluated with the result that 


+ < 8 <8) (23) 


& Tacs) 
(a—1)/a J__*A@) _ 
21) 1 


where 
x(82’, 81) — x(8, 81) 
1 — #1) 


and J, is used for the full expression 7,(1 — 8, 8). 


A(s) = 


(25) 
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— 


7 
3426 
4ous 
5477 5550 


7810 


8336 
8896 
9595 


= %3,62760 


In accordance with (14) the total rate of heat transfer from 
Strip 2 is 


Ke, 


n 


+ (1 — B) dot (26) 


The first term in the braces represents the heat transfer that would 
occur if the wall temperature from z = 0 to z = 2,’ were uni- 
formly 6,-[1 — I The second term represents the effect 
of the temperature gradient along the adiabatic portion. A 
closed expression for the integral has not been found. Calcula- 
tions have shown that the contribution of the second term is rela- 
tively very small unless the distance between the two strips is 
quite short. 


5 Description of the Experimental Apparatus 


The wind-tunnel has an open circuit with a fan at the intake 
end. The test section is 32 in. wide by 22 in. high. The contrac- 
tion ratio is approximately 7:1. The air speed is regulated by 
varying the pitch of the blades with the aid of a hydraulic mecha- 
nism that can be operated while the tunnel is running. 

The mainstream turbulence intensity was measured by means 
of a hot-wire anemometer. Expressed as the ratio of the root- 
mean-square of the longitudinal velocity fluctuation to the main- 
stream mean velocity, its value was found to vary from 0.9 to 
1.35 per cent, depending upon the station and the air speed. At 
71 ft/sec, the velocity for most of the mass-transfer tests, the 
turbulence intensity varied from about 1.1 to 1.3 per cent. 

The flat plate spanned the tunnel test section in a horizontal 
position. It was made of */,-in. thick aluminum; its chord and 
over-all span were 15 and 38*/, in., respectively. The leading-edge 
radius was in. 

The plate, m unted on heavy stanchions outside the tunnel, 
could be rolled into and out of the test section through the vertical 
tunnel walls on special tracks. Thus the test mosaics could be in- 
stalled and disassembled outside the tunnel, and the undesirable 
transient losses that occurred while the wind was being ac- 
celerated to or decelerated from test conditions were minimized. 
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Some details of the test plate are illustrated in Fig. 3. Test 
trays slide into the channels from the rear. The two rows of trays 
allow two simultaneous tests to be run at the same speed 
and temperature. Each row of trays is provided with two rows 
of diffusion guards, one row on each side of the channels. The 
configuration of the guards matches only one set of test mosaics. 

The mosaics were assembled from steel trays of the type shown 
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Fig. 3 Plan and profile views of the test plate 


Transactions of the ASME 


3°43 

3 0 1 2 3 4 5 6 a — 
0,0 0,0000 0,1780 2246 2573 2635 3057 3252 3732 

0.1 3869 3998 4119 423k 4833 
0.2 4921 5006 5090 «5250 «5327 S403 5621 
0.3 5691 5759 5827 5893 5959 6023 6087 6150 6212 6273 : E 
0.4 6334 6453 6512 6570 6627 6684 6742 6797 6852 
0.5 6908 6963 7017 7072 7126 7179 7233 7286 7339 7392 
0.6 7497 7550 7602 7655 7706 7758 MEM 7862 7915 
0.7 7967 8019 6072 612) 8177 8230 68263 8390 Bo 
0,8 8498 8553 8609 8665 8722-8779 «8837 8957-9028 
0.9 9081 «95508 9691 9805 

1,0 1.0000 + 
4 
: 

| 

- 15" 


U (ft/sec) 


4 
x (in) 


Fig.6 Pressure distributions on the test plate 
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Fig. 5 Orientation of the test plate in the wind tunnel 


in Fig 4. There were 14 trays 1.000 in. in length. Their effective 
spanwise widths were 4.754 in. +1.1 per cent; the individual 
widths were used in the calculations. Blanks 1.000 in. in length 
were employed to introduce any additional inert (analogically un- 
heated) lengths. Some longer trays are briefly described later. 

A plate assembled from many individually machined and cast 
pieces must, of course, be expected to have a sort of waviness 
superposed on whatever roughness the basic materials possess. 
In the present case this difficulty is augmented by the regression 
of the naphthalene as it sublimes. This regression is not uniform, 
and one must expect that it impresses its peculiar type of rough- 
ness in different ways as time progresses during any run. The 
step-type roughnesses of contiguous trays was as much as 0.003 
in., and in some cases the regression of the naphthalene appears 
to have been as much as 0.015 in. at the leading edges of some 
strips. But on the average the change of depth in a leading tray 
of '/2 in. subliming length was in the order of 0.001in. A system- 
atic effect of the surface roughness due to the protrusion of the 
naphthalene in the “‘as cast’’ (unsanded) condition was observed. 
The variation of the laminar transfer with respect to this rough- 
ness was small; and it was minimized by sanding the specimens 
flush. 

It will be seen from the boundary-layer velocity profiles and the 
mass-transfer results themselves that, within the scope of the 
experiments, most of the assembled surfaces on the plate were 
smooth in the hydraulic sense. 

The plate was oriented as shown schematically in Fig. 5. After 
establishing that the flow over the front guide vane was com- 
pletely attached and laminar, the baffle was lowered until no 
separation occurred at the leading edge of the test plate. 

Because the flow was split by the front guide vane, the main- 
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stream velocity U was measured above the test plate. In the 
region of the plate where the specimens were used (the leading 
7.5 inches) the average value of the dynamic head was found to 
be 31.4 per cent of the total head, po, scattering within +3 per cent 
over the entire range. This value was used to reduce the pressure 
distributions represented in Fig. 6. 

The pressure coefficient is defined by 


P — Pm (27) 
Pam 
where p is the local pressure measured at taps along the center 
of the plate, p,, is the mean static pressure over the main part of 
the plate, and p,,, is the corresponding mean dynamic pressure, 
0.314po. Consequently, 


C, = 3.18 2.18 (28) 
Po 

Fig. 6 shows typical pressure distributions observed at velocities 

of about 41, 71, and 92 ft/sec. The drop in C, over the leading 

12 in. is almost linear, and its value is about 0.045 ft~*. In cor- 

relating the mass-transfer data, a mean velocity above each tray 

was used. 

In order to ascertain that the boundary layer was laminar, 
velocity profiles were measured at several positions above the 
plate. The profiles were measured by means of a 0.035-in. thick 
total-pressure probe in conjunction with a Betz manometer. 
The trays were filled with clay and made smooth. In one test the 
clay in a leading tray was gouged to simulate a naphthalene speci- 
men that would have sublimed for an exceedingly long time; the 
downstream flow exhibited no effect within the sensitivity of the 
instrumentation. 
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The result of eight profiles located at distances of 3.78 and 7.5 
in. from the leading edge are presented in Figs. 7 and 8. The 
flags refer to the right and left-hand sides of the test plate.‘ 

By comparison with the Blasius profile in Fig. 7, it is concluded 
that the boundary layers were laminar. 

The boundary-layer thickness 6 in Fig. 8 is defined to be the 
distance from the surface to the point where the velocity reached 
99 per cent of ite local mainstream value. The experimental re- 
sults are well bracketed by two analytical profiles which have 
been employed in approximate boundary-layer calculations for 
laminar transfer on a flat plate and which lead to Eckert’s step 
function. 


6 Procedure and Scope of the Experimentation 


The procedure employed in this investigation is essentially the 
same as that described in [13 and 14]. Each tray was weighed 
in an analytical balance before and after it was exposed to the air 
stream. The time of exposure was measured and so were the air 
speed and temperature. To account for the transient losses that 
occur during handling, auxiliary tests were performed at the be- 
ginning and the end of each main test. Thus each tray was 
weighed at least 4 times. The corrections for the transients 
amounted to about 1 mg per '/; in. subliming length. 

Since the losses of weight were determined to +0.1 milligram 
and the minimum amount of mass transferred from any tray dur- 
ing the main tests was about 10 milligrams, the maximum relative 
error due to the repeated measurements was approximately 4 per 
cent. Usually, the net transfer was 20 to 30 milligrams, and so 
the maximum relative error usually was about 2 per cent. 


* The observer is oriented with his back to the stream. 
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Fig. 8 Boundary-layer velocity profiles compared with polynomial approximations 
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Most of the tests were performed at a mainstream velocity of 
about 71 ft/sec corresponding to a maximum Reynolds number of 
about 250,000. A few tests were run at 41 and 92 ft/sec. 

The mean stagnation temperature of the air was evaluated by 
arithmetically averaging the emf readings taken during each 
run. The temperatures of all the tests were between 72 and 88 F. 
They were measured with calibrated thermocouples in conjunc- 
tion with a potentiometer, accurate to within 0.2 F. The rate of 
temperature rise was about 0.15 Fahrenheit degree per minute, 
and the maximum variation during any run was about 3 F. The 
exposure times were varied from 6 to 20 minutes. In some cases, 
the specimens were exposed more than once. 

Data were obtained for single trays, two trays in tandem, and 
multiple (more than two) trays in tandem. 


7 Reduction of the Data 


The observed rate of mass transfer per unit span of the jth tray 
or strip is 


where Am, is the net mass transferred from the jth strip after a 
correction for the handling transient, Ar is the exposure time, and 
w, is the spanwise width of the strip. 

In order to compare the observed and theoretical rates of mass 
transfer, and, in particular, to check the simple step function (18) 
for the total transfer, the so-called “absolute” transfer ratio 
th ;'/tio;' was formed. The quantity ri,’ is the transfer rate per 
unit span that would occur if the entire surface area from z = 0 
to x = z; were subliming at the conditions of velocity and tem- 
perature in the test. It is calculated according to the following 
equation, which defines the average coefficient of mass transfer 
z,). 


= 0(0, "(Pow — (30) 


The factors in the right-hand member are evaluated as discussed 
in the next paragraphs. 

The mass-transfer coefficient is calculated by means of the 
equation 


4) an 4 
U Se ( ui) (31) 

v 
where U is the velocity observed at the outer edge of the boundary 
layer. From Pohlhausen’s theory, C = 0.664andn = 0.5. The 
distance z, is measured from the forward tip of the plate to the 

trailing edge of the jth strip. 

The vapor pressure in the mainstream is nil, and so p..o = 0. 
The vapor at the surface of the subliming wall is assumed to be in 


equilibrium with the solid and, also, to be thermally perfect. 
Thus 


Pow 
RT, (32) 
where 
B. 


As in [13], Sc = 2.5 and R, = 12.05 ft lb/lby F; B, = 11.884 
and B,; = 6713 for absolute temperature in degrees Rankine and 
pressure in lb/ft*. Further, the wall temperature is taken to be 
equal to the observed temperature; for the effects of cooling by 
sublimation and heating by friction are both small and, besides, 
virtually cancel themselves at the air temperatures and speeds 
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employed in these experiments. Hence (30) may be written 


0.5 
tha,’ = ‘Powe (34) 


8 Discussion of the Results 


The main tests are comprised of runs in which the specimens 
were sanded. In the next three sections the results on the single, 
double, and multiple strips are discussed. The influences of 
roughness and lateral (spanwise) diffusion are noted. 

8.1 Single Strips. The results of the leading strips are presented 
in Fig. 9. The solid line represents (18). Except for the results 
on the longest strips, the theoretical curve correlates the main 
daia within the expected scatter due to experimental error. 

The effect of the roughness of the unsanded specimens is to 
raise the rates of mass transfer above the theoretical values. It 
was particularly large for a specimen which was made to protrude 
0.005 in. above the general level of the plate. The deviation from 
the theoretical value in that case was about 7 per cent. The effect 
systematically diminished as the protuberance regressed during 
repeated exposures. 

There was a residual effect of velocity in the experimental setup. 
It is shown particularly by the abnormally high rate of mass 
transfer from a tray of 6.5 in. subliming length (Tray 10) at 41 
ft/sec. The effect evidently was a property of the tunnel and 
of mounting the plate independently of the guide vanes. At the 
higher speeds the residual effect disappeared. 

The deviation of the experimental points from the theoretical 
curve at z,’/z,; = 0.926 may be ascribed to the fact that the pres- 
sure gradient at 92 ft/sec was slightly adverse as shown in Fig. 6, 
thus reducing the shear and the rate of mass transfer. The scat- 
ter of the points at 2’/2, = 0.038 probably is a consequence of 
the difficulty of keeping such large surfaces free of waviness, 
which could have been avoided, it is believed, if the trays had 
been made deeper. 

The lateral diffusion increased the rate of mass transfer by 
thinning the diffusional boundary layer. Since its influence could 
be most easily found if it were acting along large chordwise dis- 
tances, Tray 10 was run with a reduced width, the edges being 
filled flush with an inert clay. The result is shown by the par- 
tially shaded point marked “w ~ 2.5 in.’’ which lies about 5 per 
cent higher than the theoretical value. Accordingly, the cor- 
responding deviation for the wider strip would be less than 5 per 
cent, but a definite conclusion cannot be drawn on account of the 
afore-mentioned scatter. 

With regard to the trays of '/; in. length, simultaneous tests 
on guarded and unguarded strips and on strips of 4.75 in. and 
2.5 in. span show no significant differences in the rates of mass 
transfer from corresponding specimens. In view of these results 
and of those on the long tray, it is concluded that the effect of the 
lateral diffusion from the 1/:-in. trays is considerably less than 5 
per cent. Hence, for the present purpose it has been neglected, 
and the '/;-in. specimens are regarded as though they are strips of 
infinite span. Further, two strips in tandem as well as the single 
strips can be viewed in this way. 

8.2 Two Strips in Tandem. As soon as two strips in tandem are 
used, the number of possible mosaics becomes very large. In order 
to keep the program within bounds, only two sets of such mosaics, 
comprised of !/;-in. strips, were tested. 

In one set of tests the distance between them was kept at '/sin., 
and they were moved together from one position to another in 
different runs. The results are shown in Fig. 10, where the 
abscissa is the distance to the center of the configuration expressed 
in units of the subliming length. 

A noticeable effect of the roughness in the interposed steel 
blanks was observed when the trays were at their greatest dis- 
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Fig. 10 Comparison of experimental and theoretical relative ratios (fixed 
distance between strips) 


tance downstream. Thus, when the 1-in. blanks were carefully 
selected so that the forward inert length was smooth to the sense 
of touch, the experimental result fell in good agreement with the 
theory, as shown by the lowest point at z,,/l = 6.5. 

In the second set of tests on two '/,-in. strips in tandem, Strip 1 
was held fixed at z,; = 0.75 in. and Strip 2 was moved to different 
positions downstream. The results are shown in Figs. 11 and 12 
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Comparison of experimental and theoretical results for single strips 


where the abscissas are distances between the strips measured in 
units of the length of one strip. 

The ordinate in the graph of Fig. 11 is the relative mass transfer 
ratio ti:'/m,’. Interestingly, the theory shows that the mass 
transfer from Strip 2 is increased 20 per cent when it is moved to 
the distance | downstream. It reaches a maximum at about 2.51 
and thereafter it hardly changes. The experimental results sub- 
stantially verify the results predicted by the theory. Again, the 
low-speed runs appear high. The scatter is somewhat greater 
than it was for the single strips, especially for the prelimi- 
nary results, undoubtedly because the roughness of both strips 
influences the transfer process. 

The absolute mass-transfer ratios, 2’/1ig’, for the same series 
of tests is presented in Fig. 12. Theoretically, the absolute ratio 
reaches a maximum when the distance between the trays is 1/2, 
and again the agreement between the experimental and theoretical 
results is good. 

8.3 Multiple Strips in Tandem. Several tests were performed on 
mosaics comprised of seven '/,-in. strips in tandem. In each 
test the inert distance to the leading edge of the first strip was 
x,’ = 0.25 in., and the total distance to the trailing edge of the last 
strip was z; = 6.75 inch. The inert distances between the strips 
were '/, in., and each strip was situated on an individual tray. 
The velocity was approximately 71 ft/sec in all the tests. 

In order to determine whether the nature of the flow in the 
boundary layer was grossly altered by the changing roughness, 
the pressure behind the last strips was continually observed by 
laying the impact tube flat on the test surface. A pressure 
variation occurred during the first few minutes of each initial ex- 
posure, but its magnitude was quite small and no particular sig- 
nificance could be attached to it. In one instance a complete 
boundary-layer velocity profile, measured about '/; hour after 
the sublimation test had ended, virtually coincided with the pro- 
files described previously. Hence, it was concluded that within 
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Fig. 11 Comparison of experimental and theoretical relative ratios (fixed leading strip and variable 
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the limitation of the instrumentation the profiles remained sub- 
stantially laminar during the main tests and that no major change 
in the nature of the boundary layer occurred during ablation. 
However, significant changes probably would have been found 
with finer instrumentation as indicated by the mass-transfer re- 
sults themselves. 


n 
In Fig. 13 the total mass-transfer ratio se m,'/ton’ is plotted 
j=l 
against the number of strips from n = 1 to 7. Curves are not 
drawn through the points because the functional relationship is 
discontinuous. The dotted curve represents the upper bound of 
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Comparison of experimental and theoretical absolute ratios (fixed leading strip and variable 


the experimental results of very rough specimens, details of which 
are omitted from the graph in order to avoid obscuring the per- 
tinent features. Almost all exposure times were 10 minutes. The 
suffixes A and B refer, respectively, to the first and second ex- 
posure. The points are staggered to indicate the side of the test 
plate on which the strips were situated. 

Clearly, there is a remarkable spread of the results, increasing 
as the number of the strips is increased. This spreading is ascribed 
to the effect of the changing roughness. 

The results of the second exposure are almost consistently 
higher than those of the first, the possible exception being those 
at Strip 1. Hence, as one expects intuitively, the effect of the 
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increasing roughness is to increase the rate of mass transfer. This 
is similarly confirmed by the results on unsanded specimens 
which were rough from the outset. 

Since the surfaces ablate during the entirety of each exposure, 
the results presented in Fig. 13 must be regarded as averages for 
the duration of each exposure. If it is assumed that the change 
is linear with respect to time, the initial values may be estimated 
by extrapolating back to half the difference between two ex- 
posures on the same mosaic. Accordingly, the initial values cor- 
responding to Strips 6 and 7 in the tests represented by the open 
symbols would be about 0.02 less than the values plotted for the 
first exposure. 

The lowest points on the graph were obtained with diffusion 
guards, so that the simulation to strips of infinite span was better 
than with the others. It appears that the total mass-transfer 
ratio reaches a maximum of about 0.74. Supposing that the 
effects of roughness are the same with or without the diffusion 
guards, we may conclude that for the given configuration in 
laminar flow the maximum value of the total mass-transfer ratio 
for smooth strips of infinite span is about 72 per cent. The rela- 
tively smooth, unguarded strips possibly would reach maximum 
values beyond n = 7, but the rough ones probably lead to early 
transition. 

What happens when the inert distance between successive 
strips is reduced? Tests were performed on the single aluminum 
tray with 11 strips 0.100 in. apart (Tray 11). Like the 7-strip 
mosaic, the leading edge of the first strip was at x,’ = 0.25 in. and 
the trailing edge of the last strip was at z,, = 6.75inch. The data 
were obtained without guards, and all the strips could not be 
finished as well as desired on account of their large number and 
of the softness of the aluminum. For these reasons the results on 
the 11-strip mosaic, which lay between 0.88 and 0.95, must be 
regarded as preliminary. As could be expected, the increase from 
a nominally 50 to 80 per cent active area increased the total mass 
transfer. The trend is similar as before, namely, the ratio in- 
creases with increasing roughness; but without additional tests 
using smoother specimens and diffusion guards, it would be im- 
possible to separate the effects of roughness and spanwise diffusion. 


The conclusions are based on the results from the initial ex- 
posures of specimens that had been made relatively smooth by 
sanding. The scatter introduced by the remaining roughnesses 
is well within the maximum experimental error, which is 5 per 
cent. 

Single and double strips of aspect ratio 9.5 show no significant 
effect of spanwise diffusion and can be considered of infinite span. 
The influence of spanwise diffusion is observed, however, on the 
long narrow subliming surfaces as well as on equally long mosaics 
of multiple strips, each of aspect ratio 9.5. It increases the total 
rate of mass transfer about 5 per cent of the values predicted 
for a single strip of the same length but of infinite span. No 
attempt has been made to account for the lateral diffusion by cor- 
relating against the aspect ratio in a general way. Here, a theory 
is needed that would, also, account for the fact that the thermal 
diffusivity of the air is considerably greater than its diffusivity 
with respect to naphthalene vapor. 

The results on 7 equally spaced strips of virtually infinite span 
indicate that an active area of nominally 50 per cent would yield 
about 72 per cent of the mass transfer from a surface of uniform 
concentration functioning otherwise under like conditions. Some 
preliminary results on 11 strips indicate, however, that when the 
active area along the same distance is increased to nominally 80 
per cent, the corresponding mass transfer is in the order of 85 per 
cent. Clearly, the surface discontinuities of concentration and 
flux due to the short inert surfaces between the strips neither 
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considerably deter nor appreciably improve the rate of mass 
transfer. 

By analogy the surface discontinuities of temperature and flux 
in similar heating arrangements, comprised of alternate isothermal 
and adiabatic surfaces, would not significantly affect the rate of 
heat transfer. 

The good agreement between the theoretical and experimental 
results on the strips of infinite span indicates that the RKT 
method using Eckert’s step function for the laminar boundary 
layer on a flat plate at zero incidence is satisfactory for engineering 
purposes. 
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Introduction 


| = EVALUATION of the heat transfer to and from 
clothing requires a knowledge of the absorptance, reflectance, 
transmittance, and emittance of the cloth. Such information 
may consist of total radiant property measurements which are 
applicable only for specific conditions, or it may consist of spectral 
data which can be utilized to predict the total radiant properties 
for sources of any spectral distribution. This investigation was 
concerned with the determination of the spectral reflectance of 
several types of natural, synthetic, and regenerated textiles used 
as clothing materials. The types of textiles are defined by the 
fibers used in their construction. Natural fibers are those which 
are derived directly from animals, plants, or insects, e.g., wool, 
cotton, silk. Synthetic fibers are produced in the laboratory from 
simple chemical units and have no relationship to natural fibers. 


Contributed by the Heat Transfer Division of Tae American 
Sociery or MecHANnicat ENGINEERS and presented at the ASME- 
AIChE Heat Transfer Conference, Storrs, Conn., August 9-12, 1959. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, May 4, 
1959. Paper No. 59—HT-9. 


Table 1 Description of cloth samples 


Spectral Characteristics of Fabrics 
From 1 to 23 Microns 


Nylon and Orlon are two examples of this type. Regenerated 
fibers are made from relatively complex chemical units, e.g., 
cellulose or protein, and are frequently made by treating natural 
fibers chemically to produce a new type of fiber. Rayon and 
acetate are examples of regenerated fibers. Table 1 lists the 
samples tested, groups them into their main categories, gives a 
short description, and indicates whether the samples were the 
same on both sides. 

The work described in this paper was performed at the Uni- 
versity of California at Berkeley as an outgrowth of a research 
project sponsored by the United States Army Quartermaster 
Corps Research and Development Center. The spectral re- 
flectance as reported here is actually the system reflectance, since, 
in addition to energy reflected from the sample, it also includes 
energy that passes through the sample and is reflected from the 
low reflectance background and transmitted back through the 
sample. A low reflectance background was used in these tests; 
98 a consequence, the effect of the background was negligible ex- 
cept in the region of 1.0 to 3.0 microns where the textile material 
is relatively transparent to infrared radiation. It was found ex- 
perimentally, for the materials reported in this paper, that three 
layers were sufficient to eliminate the effect of the background, so 
data were reported for both single and triple layers of cloth. 


No. of 
layers 
Fig. Material Class Color run Remarks 
3 Cotton Natural 1&3 Close weave avg. of sam- 
ples, dyed, and undyed 
4 Cotton Natural Medium gray 1&3 Sateen, treated w/Car- 
boxy Methyl cellulose 
and impregnated with 
carbon 
4 Cotton Natural Dark gray 1&3 
4 Cotton Natural Deep Black 1&3 
4 Milium Natural Aluminum 1&3 Aluminized cloth, close 
weave 
5 Linen Natural Blue 1&3 Close weave, thin 
6 Wool Natural Spotted pink 1&3 Close weave, thin 
7 Wool-nylon (85-15) Nat.-syn. ark green 1 Close weave, thick 
7 Wool-cotton (50-50) Natural White 1&3 Close weave, thick 
8 Silk Natural White 1&3 Close weave, thin 
9 Acetate Regenerated White 1&3 Close weave, thin 
10 Arnel Synthetic Rust 1&3 Close weave, thin 
ll Rayon Regenerated Green 1&3 Med. weave, thin 
12 Dacron Synthetic White 1&3 £4Med. weave, thin 
12 Dacron Synthetic Green 1&3 Med. weave, thin 
13 Orlon Synthetic White 1 Fuzzy ae, suit ma- 
teria! 
13 Orlon Synthetic Blue 1 Fuzzy matin, suit ma- 
teria 
14 Nylon Synthetic White 1&3 Med. weave, thin 
14 Nylon Synthetic Green 1&3 # £Med. weave, thin 
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Fig. 1(b) Block diagram of system 


Apparatus and Equipment 

A Gier-Dunkle black body reflectometer [1]! is used, which 
operates by exposing the sample to black body radiation within a 
“hohlraum”’ at a temperature of approximately 1400 F, Fig. 1(a). 
Since the sample is cooled to approximately detector tempera- 
ture, the reflectance can be evaluated from the ratio of the radios- 
ity of the sample to the radiosity of the hohlraum. The radiosity 
is defined as the flux density leaving a surface, and includes 
transmitted, reflected, and emitted radiation. For opaque ma- 
terials, the radiosity is equal to the sum of the emitted and re- 
flected radiation fluxes. A block diagram of the system is shown 
in Fig. 1(6). Radiation from the black body is sent through the 
optical system to the Perkin-Elmer monochromator where a 
prism breaks the infrared radiation into its spectrum. The radia- 
tion is detected by means of a sensitive vacuum thermocouple, the 
signal of which is transmitted through an amplification system 
and then to a recorder. 

A slight modification was made in the sample cooling method 
described in reference [1], Fig. 1(c). The center tube of the sam- 
ple holder which normally directs water on the rear of a metal 
sample, was connected to a tank of gas through a flow-metering 
system. This allowed gas to flow through the cloth samples, 
thus cooling them by convection. 

An inert cooling gas, which is transparent to thermal radiation, 
and is relatively inexpensive, is desirable for this purpose. The 
majority of texts were run using dry nitrogen as a coolant. 


Experimental Procedure 


The reflectance of the sample holder was taken without a 
sample in place to provide some basis for evaluating the back- 
ground effect. The reflectance is found to be reasonably constant 
ever the wave-length range from 3 to 23 microns (see Fig. 2). 
The variations shown in the figure are attributed to the existence 


1 Numbers in brackets designate References at end of paper. 
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Fig. 1(c) Sketch of sample holder 


of surface films on the central cooling tube in the sample holder 
The central tube was constructed of copper and had an enlarged 
end built up of solder to direct the water flow onto the sample for 
normal operation with metallic samples. This tube should have 
been blackened for use with the cloth samples to reduce reflec- 
tions, but this was not done. The nature of the corrosion prod- 
ucts in this tube was not known, but thev were present in small 
amounts. The effect of the variations in the constancy of the 
background reflectance is negligible. 

Samples used were cut from large swatches of material and 
trimmed to a circular shape of 7/s in. diameter. These were 
placed in the retaining cap and held in place by an O-ring. 

The gas flow was started and the sample holder inserted in the 
black body reflectometer. The coolant gas entered the sample 
holder at room temperature, which was controlled at 68 F. The 
sample temperature could not be reliably measured with the cool 
gas flowing through it while exposed to the high radiation flux on 
one side. The gas flow rate through the sample was found to be 
critical. If the flow rate was too large, the bottom of the cavity 
would be cooled, resulting in a low value for the reflectance. The 
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temperature distribution in the cavity was unknown due to ip- 
adequate instrumentation. Several thermocouples had been 
mounted in the walls of the cavity, but were inoperative at this 
time. Visual inspection of the cavity, which is very sensitive to 
small temperature differences in this temperature range, indi- 
cated only small variations in brightness over the cavity walls. 
If the gas flow rate was too small, the sample would heat up and, 
due to emission from the sample, a high value of the reflectance 
would be obtained, and/or the sample was affected by the high 
temperature. The adjustment of the flow rate was made by set- 
ting the monochromator in a low reflectance region of the spec- 
trum (usually at the longest wave length available with the prism 
used), and then slowly reducing the gas flow rate until a slight 
increase occurred in the recorder deflection. A slight increase in 
the flow rate above this point yielded the optimum flow rate, 
which was found to be approximately 0.3 cubic feet per minute. 
This flow rate was satisfactory for most materials, except those 
with fibrous or fuzzy surfaces, such as mohair (not included in this 
paper ). 

After adjusting the gas flow rate through the sample, three re- 
corder deflections were obtained for each wave length; namely, 
“viewing” the sample, viewing the cavity wall, and viewing 
a black surface at room temperature (a zero reading). 

The reflectance was determined at every 0.25 micron from 1.25 
microns to 23 microns. The lower wave-length limit was fixed by 
the monochromatic emissive power avsilable from the 1400 F 
hohlraum and the sensitivity of the system. The upper limit was 
chosen for several reasons: One, that 78 per cent of the energy lies 
below 23 microns for a temperature of 60 F; two, the energy 
from the black body became so small that the slit widths used be- 
came too large for good spectral resolution; and three, the uni- 
formity of the curves in the far infrared indicated the possibility 
of extrapolation to longer wave lengths with a small amount of 
error. This method of taking data at 0.25 micron intervals was 
determined by Gier, Dunkle, and Bevans [1] to be the most 
feasible since it allowed the amplifier gain and slit widths to be ad- 
justed for maximum deflection, thus affording less error than 
would be the case with a continuous scanning and comparison 
process. Because of a fairly high noise level the recorder deflec- 
tion was rounded off to the nearest integer on the scale of 100. 

In addition to the error in the recorder reading there are small 
errors due to emission from the cloth and cooling of the bottom of 
the hohlraum cavity. The error is not computed as a percentage, 
but only as units of reflectance (see Appendix). 


Discussion of Results 


The spectral reflectance curves for the various materials tabu- 
lated in Table 1 are plotted in Figs. 3 to 14. For most materials, 
two curves are plotted, one for a single layer of cloth with the low 
reflectance sample holder as a background (Fig. 2), and a second 
curve for a triple layer of cloth which was found to be sufficiently 
thick to eliminate the effect of the background for all materia!s 
tested. As all the cloths tested had low reflectances (below 50 
per cent) beyond 3 microns, the curves for single and triple layers 
should be nearly identical at wave lengths greater than 3 microns, 
and this correspondence acts as a check on the consistency of the 
results. 

The infrared spectrum is unique for each particular molecule, 
hence in the textile samples observed one would expect to find a 
unique spectrum for each particular type of textile. However, 
since we are dealing with a specialized section of organic com- 
pounds where certain basic molecules are inevitably encountered, 
one expects and finds a number of similarities in the infrared 
spectrum which should be readily apparent. All the textiles 
show absorption lines or bands at 3-3.5 microns, 5.75-6.00 mi- 
crons, 7.25 and 7.75 microns, 7.00-9.00 microns, and 14-17 
microns. Since this investigation is not concerned with infrared 
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spectroscopic identification we will not endeavor to identify the 
causative agents beyond the fact that certain alkanes, alkenes, 
alkynes, and esters, and the carbon-hydrogen bond are probably 
the causative molecular groups. Colthrup [3] gives a chart of the 
correlation between the molecular structure and infrared spectra. 

The reflectance of a cloth consists of the contributions from in- 
terreflections within the cloth itself as well as surface reflectance 
(6). This means that in regions where the cloth material is 
transparent a higher reflectance will result due to these interreflec- 
tions or due to background reflectance. In most materials the 
region below 2 microns corresponds to a “window” in the 
material so that a high reflectance is observed in this region. 
On the other hand, in regions where a bigh absorptance occurs 
(absorption bands in the material) there will be little or no con- 
tribution due to the interreflections, and a low reflectance will 
result, which will be little affected by the background. 

It should be noted that the use of dyes on textiles does not materially 
affect the infrared spectrum past 8 microns on any of the textiles ex- 
cept in the case of Dacron (Fig. 12). This is primarily due to the 
fact that the dye is present in relatively small quantities com- 
pared to the base material so that the amount of dye the fiber 
takes up is not sufficient to cause any perceptible effect in the in- 
frared region. The dye has a large effect in the visible region due 
to the resonant electronic structure of the dye molecule. From a 
study of the gross spectrum shown it would appear that the dyes 
present do not alter the molecular structure of the textile mole- 
cule. In the case of the Dacron, the green dyed specimen de- 
viates considerably from the white sample, and it is not known 
whether this is actually due to the effect of the dve or to a dif- 
ference in the chemical structure of the two types of Dacron. 

It should also be noted that several other processes on textiles 
do change the infrared spectrum considerably. Metallizing the 
cloth by coating the fiber of the fabric with some metal such as 
aluminum, gold, or silver causes the cloth to exhibit more of the 
properties of the metal than of the basic textile. This is evidenced 
by the almost constant, relatively high, reflectance of Milium 
(aluminized cloth, Fig. 4) from 1 to 23 microns, broken only in 
portions of the spectrum by the strong carbon hydrogen bands. 
Cloth samples impregnated with carbon black (Fig. 4) show 
effects through the complete infrared region as well as the visible. 
The effect of carbon black is to minimize by absorption the con- 
tribution of the interreflections within the cloth, so that the 
principle result is to reduce the reflectance in regions where the 
reflectance is high. 

It is also of interest to note that the three main classifications of 
textiles (i.e., natural, regenerated, and synthetic) show certain 
basic similarities and differences in their infrared spectra. 

For example, in the natural textiles two major maxima occur in 
the reflectance curve at approximately 5.25 and 11.5 microns. 
The regenerated fibers also have two maximums, although they 
are located at 4.25 and 12.75 microns. The synthetic textiles 
differ in their spectral characteristics from the regenerated and 
natural groups and also differ from each other. 

From the nature of these basic similarities and dissimilarities it 
is easy to classify an unknown textile sample into one of the three 
major groups. 

The spectral reflectance curves can be used to compute the 
total emissivity of the cloth [4, 5]. The curves can also be 
utilized to compute the reflectance for sources of any temperature 
or energy distribution [4]. If the spectral reflectance curves are 
extended through the visible region, it is found that the reflec- 
tance in this region is determined almost exclusively by the type 
of dye used on the cloth. These data can be used to calculate the 
reflectance of materials for solar radiation. A paper is now under 
preparation describing apparatus for direct measurements of solar 
transmittance and reflectance of textiles. Good agreement is ob- 
tained between the measured and calculated solar reflectance. 
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Fig. 10 Reflectance, Arnel 


FEBRUARY 1960 / 6] 


ny 
| 
70) 
| 
so 
ara 
Fig. 3 Reflectance, cotton, close weave 
| = 
3 / 
2 
oo 
4 
ie 
1 ( 
te? 
of 
| 
i 
Fig.6 Reflectance, wool 
| Journal of Heat Transfer 
? 


(GREEN) 


SNGLE Laven 
Laver 


wom 
Wo veleng!® — Mrerene 
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Fig. 12 Reflectance, Dacron, white and green 


In the spectral region reported herein, it should be noted that 

the atmosphere contains two absorbing-radiating substances, 
water vapor, and carbon dioxide (CO,). Ozone also has absorp- 
tion bands in the infrared region, but its influence in the lower 
atmosphere is negligible. The CO, and water vapor bands are 
shown schematically in Fig. 15 [2], but it must be remembered 
that this schematic is not a rigorous delineation of the absorption 
bands since the concentration of water vapor and CO, in the at- 
mosphere varies between wide limits. This wide variation can 
produce broadening or attenuation of the band limits as well as 
increasing or diminishing the percentage of absorption and emis- 
sion [2]. 

If we neglect radiation in the shorter wave lengths (night time), 
such as soler radiation, we see that the net gain or loss of heat 
from the textile is strongly influenced by the atmospheric infrared 
spectrum. For instance, if a textile emits radiation in one of the 
atmospheric absorption bands the radiation will be absorbed and 
some of it reradiated beck to the textile. This, of course, means 
a reduced net heat transfer as compared to the net heat transfer 
when the same smount of radiation is emitted in the window re- 
gions of the spectrum. Analogously, the same is seen to be true 
of radiation reflected by the cloth for the different spectrum re- 
gions of the atmosphere. 

It can be seen by inspection of the data that certain textiles 
have spectral distributions with relatively high reflectance in the 
absorption bands of the atmospheric infrared spectrum and pre- 
dominately low relative reflectance (high emittance) in the 
window regions. Such cloths as linen, Dacron, and Arnel fulfill 
these specifications. Such cloths as wool-cotton mixtures, wool- 
nylon combinations, and pure wool have the reverse spectral dis- 
tribution. Finally, some cloths have their high and low reflec- 
tance bands distributed randomly in and out of the atmospheric 
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Fig. 13 Reflectance, Orion 


(WHITE) 
DUPONT WOK 


oo — SINGLE LAYER 
~~DOUBLE LAYER 


NYLON (GREEN) 
DUPONT NO. 7430~-C 


- 


~"SINGLE LAYER 
TRIPLE LAYER 


Retiectonce — Percent 
& 


Wovelength — Microns 


Fig. 14 Reflectance, nylon 


THE FAR SPECTRUM OF 
THE ATMOSPHERE TING 
OZONE) 


¢ 


Mic one 


Fig. 15 Spectral distribution of atmospheric radiation 


Absorption 


absorption bands. This includes the majority of the textile 
samples tested. The textiles in this category are cotton (and 
high percentage cotton mixtures), silk, rayon, acetate, Orion, 
Nylon, and the metallized and carbon impregnated fabrics. 

The work reported herein should be very useful to individuals 
concerned with heat-trensfer problems. The data can be used to 
compute [4] the absorptivity or emissivity of cloth as long as the 
transmittance is small. As no data are available on true trans- 
mittance of the samples, the reader should be strongly cautioned 
about assuming values in the expression 


a+r+7r=0 


in order to compute the absorptance, a. For solar radiation, or 
radiation from other high temperature sources, further data must 
be obtained in the short wave region in absorptivity and trans- 
missivity. This region was not covered in this paper as the cloth 
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properties in this spectral range are highly dependent upon the 
nature of the dyes used and are readily measured with conven- 
tional methods. 


1 J. T. Gier, R. V. Dunkle, and J. T. Bevans, “Measurement of 
Spectral Reflectivity From 1.0 to 15 Microns,” Journal of the Optical 
Society, vol. 44, 1954, pp. 558-562. 

2 W. M. Elsasser, ‘Heat Transfer by Infrared Radiation in the 
Air,”’ Harvard Meteorological Studies, no. 6, 1942. 

3 N. B. Colthrup, “Spectra Structure Correlations in the Infra- 
red,”” Journal of the Optical Society, vol. 6, 1950, p. 397. 

4 R. V. Dunkle, “Thermal-Radiation Tables and Applications,” 
Trans. ASME, vol. 76, 1954, p. 549. 

5 J. T. Bevans, J. T. Gier, and R. V. Dunkle, ‘Comparison of 
Total Emittances With Values Computed From Spectral Measure- 
ments,” Trans. ASME, vol. 80, 1958, pp. 1405-1416. 

6 R. V. Dunkle, J. T. Gier, J. T. Bevans, and D. K. Edwards, 
“Study of the Various Modes of Heat Transfer From and Through 
Clothing Materials,” Report Series 80, Issue 1, Project 92-06-001, 
University of California, Berkeley, Calif. 


7 R. V. Dunkle and J. T. Bevans, “An Approximate Analysis of 
the Solar Reflectance and Transmittance of Snow,’ Journal of 
Meteorology, vol. 13, 1956, p. 121. 


APPENDIX 
Discussion of Errors 


Errors Due to Heating of the Sample and Cooling of the Cavity 

Two significant temperature-dependent errors can occur when 
the cloth sample is cooled by passing inert gas through it. One 
error arises as a result of increased thermal emission from the 
sample if its temperature rises. A second error is due to cooling of 
the bottom of the hohlraum by the inert gas. 

Two other radiation errors can occur within the black body. 
One error is due to net radiation exchange between the sample | 
and reference area 2. This error can be shown to be very small. 
The fourth error is due to the opening in the bottom of the 
hohlraum which is important for diffuse surfaces of high reflec- 
tance but is normally negligible. 

In order to analyze the radiation exchange within the hohlraum, 
it is necessary to break up the boundaries into a number of iso- 
thermal and constant radiosity zones. 

The sample is designated as 1, the reference area on the top of 
the cavity as 2, the view port in the bottom of the cavity as 3, and 
the remainder of the zones on the cavity walls from 4 to N. In 
taking reflectance measurements, the optical system first views 
the sample 1, and then the reference area 2. The ratio of the re 
corder reading when viewing the sample to the recorder reading 
when viewing the reference area is taken as the measured re- 
flectance r,,. (Zero readings are also taken and corrected.) The 
recorder reading is proportional, at a given slit width, wave length, 
and amplifier gain, to the difference in radiosity between the 
radiation source viewed and the chopper blade. Thus, the 
measured reflectance can be expressed: 


Ju — Jw 
(1) 


The spectral radiosities of surfaces 1 and 2 are: 


Tim 


N 
= +m D> 


n=1 


N 


n=1 


Although any cooling of the wall or bottom of the hohlraum 
due to the inert gas will be greater on the sample side of the 
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cavity, the form factors of the cooled areas relative to the sample 
and reference are not greatly different. The sample and reference 
area are located symmetrically in the cavity top, and have the 
same geometrical relationship with regard to the exit opening. 

If equations (2) and (3) are substituted into (1), the expression 
for the measured reflectance becomes: 


N 

n=1 


n=1 


Thm = 


Rearrangement of equation (4) with a little algebraic juggling, 
which is not included here for reasons of brevity, the assumption 
of symmetry for the sample and reference area so that Fi, = 
F 2, and with the aid of the relationship for opaque substances 
that a + r = 1.0, yields: 


Trm — Tu = + + (5) 
where 
6, = error due sample heating 
exn(Ex — Jr) 
— Jno 


6, = error due to nonuniform hohlraum temperatures 


N 


n=4 


(Siz — Jno) 
6; = error due to view port opening in bottom of hohlraum 


(7) 


— Jrs) 


(Jrz — 


6, = error due to radiation exchange between sample and 


reference area 


(Jr2 — Jno) 


5, = (9) 


Error 1. The effect of sample heating is to increase the measured 
reflectance above the actual reflectance. This error is proportional 
to the spectral emittance and is a maximum for materials of low 
reflectance. The value of 5,/e,: is plotted in Fig. 16 as a function 
of wave length for a 50 F rise in sample temperature and a 
hohlraum temperature of 1400 F. This error is negligible at wave 
lengths less than about 10 microns. If the zero reading surface 
is at the sample temperature and black, this error is zero. 


FOR MEAN HOLRAUM TEMPERATURE So F 
REFERENCE Ty 


FOR MEAN HOLRAUM TEMPERATURE 
St.ow REFERENCE 


ron sor mse samme TEMPERATURE 
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Fig. 16 Reflectance errors occurring in hohiraum 
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Error 2. The effect of cooling the bottom part of the hohlraum 
is to reduce the irradiation of both the sample and the reference 
area, which causes the measured reflectance to be low. It is seen 
that this error tends to oppose the first error. This error is pro- 
portional to the reflectance of the sample and to the emittance of 
the reference area. This indicates that a polished metal reference 
area would be advantageous, although the rest of the hohlraum 
should have a high emittance. The part of equation (7) within 
the summation represents the effect of cavity wall cooling. If 
the walls were all black at the same temperature as the reference 
area, this term would be zero. The magnitude of this error is 
shown in Fig. 16 for the case where the cooling of the walls is 
equivalent to a 50 F drop in temperature relative to the reference 
area of 1400 F. This error is seen to increase rapidly in the short 
waves. A curve is also given for a 10-deg F temperature differ- 
ence. 

Error 3. The effect of this error is also to reduce the measured 
value of the reflectance. As J); is nearly equal to J)», this error 
becomes independent of hohlraum temperature, 


6; = (10) 


The value of F., is about 0.022 for the cavity, so this repre- 
sents the limiting error for a diffuse reflector with r,,; = 1 and 
The error is only important for diffuse reflectors and is 
negligible for specular reflectors. The reason for this is that the 
diffuse surface views the opening with a firm factor of 0.022, 
whereas a specular surface only includes radiation reflected from 
the hohlraum interior. 

Error 4. The effect of direct radiation exchange between the 
sample and reference area is negligible for the present design. 


xn = 


Other Sources of Error 


Noise Level: The noise level of the detection amplification and 
recording system was of the order of 1 per cent of full scale de- 
flection. This resulted in a possible error in reflectance of 
+0.01. 

Scattered Energy: This error is of the greatest magnitude at the 
two ends of the spectral range covered when the spectra’ energy 
in the beam is small and hence the contribution due to scattered 
energy may be significant. This source of error was controlled in 
the short waves by use of a plexiglas filter and was shown to be 
smal! at 23 microns by use of a NaC] filter. 

Most of the scattered energy is of wave lengths less than 15 
microns; hence, if a sodium-chloride filter is placed in the beam 
while the monochromator is set at 23 microns, the energy reach- 
ing the detector should be zero as sodium chloride is opaque at 
this wave length. Any signal from the system would then repre- 
sent scattered short wave energy. 

Effect of Background of Sample: This error cannot be ac- 
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curately defined since no data are available on the true transmit- 
tance of the samples; however, from the Final Clothing Report 
[6] it is seen that for a single layer of cloth with a low reflectance 
background of magnitude r, (Fig. 2): 


= Tm — — 


whence it should be apparent, since r, < 0.15 and in practically 
no case is r,, > 0.40 and is generally much lower, that: 


and in the single layer case where as a worst case we assume T = 
0.60 


Tm — T, & 0.05 


and in the triple layer case where it would be unusual for rT = 
0.20 to occur, 


Tm — T, & 0.006 (negligible) 


Therefore as a first approximation we can assume that in the 
triple layer case the reflectance measured is very close to the 
actual reflectivity of the 3 layers of cloth and that due to inter- 
reflections from the low reflectance background the system re- 
flectance of thinner samples will closely approach the reflectance 
of the triple layer samples. This is seen to be true from the 
measured results, and indicates that the background is negligible 
for the triple layer of cloth. 


Summary of Error 


The actual data error is estimated to be on the order of +0.03 in 
the reported reflectance. 


Nomenclature 
€, = spectral emittance 
E, = black body spectral emissive power Btu/(hr ft*u) 
F = form or configuration factor 
J, = spectral radiosity (Btu/(hr 
rh = spectral reflectance 
T = transmittance of cloth sample 
Subscripts 
0 = chopper blade 
1 = sample 
2 = reference area 
3 = opening in base of hohlraum 
b = background of sample 
m = measured value of reflectance 
n = zone on boundaries of hohlraum 
S = sample corrected for background effect 
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The Effect of Vapor Drag 
on Rotating Condensation 


E. M. SPARROW: and J. L. GREGG: 
Introduction 


IN THE ANALYsis of film condensation on external surfaces, 
it is common to assume that the entire process is governed 
by force and energy balances within the condensed liquid, and 
that the vapor plays a very passive role. In particular, it is al- 
most standard to suppose that the vapor offers no shearing re- 
sistance to the condensate flow and, hence, the shear stress at the 
liquid-vapor interface is zero. To lift this assumption and attack 
the problem in a rigorous way would require simultaneous solution 
of the flow problem in both the liquid and vapor. This poses an 
exceedingly difficult task which has not heretofore been carried 
out.? 

Here, we consider the problem of laminar film condensation on 
a cooled rotating disk. The centrifugal force field associated with 
the rotation sweeps the condensed liquid outward over the disk 
surface—in much the same way as gravity causes the condensate 
to flow downward in a natural condensation process. Our goal is 
to formulate and solve this problem without imposing restricting 
assumptions on the interfacial shear—and this will lead to simul- 
taneous consideration of the flow problem for the liquid and 


r. 

A detailed treatment of the condensation problem for a rotating 
disk has already been presented in a previous paper (reference 
(2]);* but there, the assumption of zero shear at the interface was 
used. The present work thus constitutes an extension and, to 
achieve a concise treatment here, we will lean heavily on the ma- 
terial contained in reference [2]. The analysis is carried out for 
condensation of a pure, saturated vapor at temperature 7',., on 
an isothermal disk at temperature 7. The only motions in the 
vapor are those induced by the rotation of the disk. 


Analysis 

The condensation process is governed by the basic conservation 
laws of mass, momentum, and energy as given by equations (1) 
through (5) of reference [2]. These laws apply equally well in the 
vapor and in the liquid. 

The fluid motions which are induced in the vapor will, in turn, 
affect the flow within the liquid; and this mutual interaction re- 
quires simultaneous solution of the equations of motion in the 
liquid and vapor. But, since we are dealing with saturated 


1 NASA Lewis Research Center, Cleveland, Ohio. Assoc. Mem. 
ASME. E. M. Sparrow is now Professor of Mechanical Engineer- 
ing at the University of Minnesota. 

? An estimate of the effects of vapor drag in the vertical plate prob- 
lem has been made by Rohsenow, Webber, and Ling (ref. [1]); but 
their method continues to deal only with the liquid film to which there 
is applied arbitrary amounts of interfacial shear. 

* Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of Taz AmeRIcAN 
Society or MecHANIcaL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, 
November 13, 1959. 
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vapor, the vapor temperature is essentially constant,‘ and the 
energy equation need not be considered in the vapor. So, once 
the velocity solutions are known, we can attack the energy equa- 
tion for the liquid and obtain the heat-transfer results. 

The first step in the analysis is to obtain solutions for the flow. 
Noting that z represents the distance measured normal to the disk 
surface, the velocity problem may be formulated as follows: In 
the condensate layer, which extends from z = 0 to z = 6, the 
mass and momentum conservation equations (1) through (4) of 
reference [2] are used with fluid properties appropriate to the 
liquid. In the vapor region, which extends from z = dtoz = @, 
the same equations apply, but now the properties are those of the 
vapor. At the interface, z = 6, there are conditions of continuity 
which tie together the liquid and vapor velocities and their 
derivatives. Using subscripts / and v to denote liquid and vapor, 
these conditions are: 


Conservation of mass across interface: (pV,), = (pV.), 
Equality of tangential velocity: (Vg); = (V¢), 

Equality of radial velocity: (V,); = (V,), 

Continuity of tangential shear: (T2g); = (T2¢)» 
Continuity of radial shear: (7,,); = (Ts)» 


Finally, the statement of the problem is completed by giving the 

boundary conditions: 
V,=V,=0 V, = 0 

z=0 


(liquid ) Ve = 0) (vapor) 


where w is the angular velocity of the disk. 

The formidable mathematical problem just outlined is made 
tractable by use of von Karman’s transformation, which reduces 
the partial differential equations to ordinary differential equa- 
tions. The transformation consists of a new independent variable 
» and new dependent variables F, G, and H which have the fol- 
lowing definitions 


n = (w/v)'/*2, F(n) = V,/rw, G(n) = Ve/rw, H(n) = V,/(wv)'/* 


The transformation applies to both the liquid and vapor provided 
that appropriate fluid properties are used and additional sub- 
scripts / and v are affixed. So, the transformed conservation equa- 
tions (7) and (8) of reference [2] apply to the liquid when subscript 
lis added and to the vapor when subscript v is added. The con- 
ditions at the liquid-vapor interface 7 = 3 (z = 5) now take the 
form® 


Vs 


where the primes denote differentiation with respect to 7. It is 
seen that the jump in velocity functions depends on the ratio of 
the pu products. From physical reasoning, it would be ex- 
pected that when the vapor density and viscosity are very small, 


Ve = TW 


4 The change in saturation pressure in the vapor due to the rota- 
tion is very small. 

5 We also use the fact that F = —H’/2 from Equation (1a), ref- 
erence [2]. 
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the vapor drag would play a minor role. Conversely, the effects 


of vapor drag should increase as (py),/(pu), decreases. The 
boundary conditions may also be transformed into the new 
variables as follows: 


H = H’ =0) H’'=0 
| 


G=1 (liquid) G=0 (vapor) 


So, with this, the simultaneous velocity problem for the liquid 
and vapor is reduced to a form where solutions can be obtained 
on an electronic computer. Before the numerical integration of 
the governing equations can be undertaken, it is necessary to 
specify values for two parameters: (1) The dimensionless film 
thickness, 9s; (2) the ratio (pu),/(pu),, which relates the liquid 
and vapor velocities at the interface. 

With the velocity problem now formulated, we turn to the 
energy equation (5a) for the liquid. The velocity solutions are 
needed as input data to this equation and, in addition, the 
Prandtl number Pr must be prescribed. So, there are thus a 
total of three parameters which characterize each solution of our 
condensation problem. In reference [2], it is shown that yg can be 
replaced by the more convenient group: ¢,;(7'sat — T\,)/hig where 
he is the latent heat. With this, the governing parameters 
may be stated as: (pu),/(pu),, Pr, ¢,,AT/hig. 


Results 


Since there are three parameters to be prescribed for each case, 
care must be exercised in choosing cases to be studied in order to 
keep the computations within reasonable bounds. Table 1 gives* 
the cases which have been considered, as well as the results for 
the heat-transfer coefficient h and for the torque M required to 
rotate the disk. As a point of departure in making these choices, 
cognizance was taken of the experiments on rotating condensation 
by Nandapurkar and Beatty [3], who worked with methanol, 
ethanol, and freon —113 at atmospheric pressure. The py ratio 
of 150 which appears in Table 1 is typical of these experiments, 
but it also applies to other fluids; e.g., water near atmospheric 
pressure. The choice of Pr = 10 is made to typify ordinary 
liquids (i.e., nonmetals), and it is not expected that the precise 
choice of Pr (for Pr > 1) will significantly influence the com- 
parison between results with and without drag. The range of 
c, AT /hy covers current practice. To test the effect of decreasing 
pu ratio (e.g., for water, this ratio decreases as saturation tem- 
perature increases), a value of 50 was used, while maintaining 
c, AT /hy at a practical value. Finally, a typical operating con- 
dition for liquid metals was studied. 

To provide perspective for the results to be presented, it is 
useful to note the following characteristics of the vapor drag: 
First, the tangential component of the vapor drag always op- 
poses the rotational motion of the condensate. So, taking ac- 
count of the vapor drag will give rise to a larger prediction for 
the torque requirement M, and this is demonstrated numerically 
in the last two columns of Table 1. Since the rotational motions 


* Unsubseripted properties are those of the liquid. 


do not transport energy, the effect of the tangential drag upon the 
heat transfer is felt only indirectly through modification of the 
radial and axial velocities. Secondly, depending on the operating 
conditions and on the film thickness in particular, the radial com- 
ponent of the vapor drag may either retard the radial flow of the 
condensate or else may augment it. It is expected that the re- 
tardation would occur for thicker films and the augmentation for 
thinner films. Since the radial flow in the condensate is directly 
involved in the transfer of heat to the surface, it is seen that 
accounting for vapor drag may either increase or decrease the 
heat transfer, depending on operating conditions. This is il- 
lustrated in Table 1. The first entry is for very thin films and 
shows an increased heat transfer. The other entries are for thicker 
films and show a decreased heat transfer. To understand the in- 
teresting behavior of the radial drag, we recall that the radial 
velocity is zero both at the disk surface (no slip) and in the am- 
bient fluid. So, a maximum value occurs somewhere between, 
and it is at this point that the direction of the radial shear force 
changes sign (r,, ~ 0V,/dz). For very thin condensate layers, the 
velocity maximum occurs in the vapor; and under this condition, 
the vapor exerts a positive radial shear force on the condensate 
layer. For thicker condensate layers, the velocity maximum 
occurs within the condensate. In this event, the vapor exerts a 
negative (retarding) shear force. 

Now, looking in greater detail at the heat-transfer listing of 
Table 1, we first turn our attention to the results for Pr = 10 
(nonmetallic liquids). For all the cases considered, the effects of 
vapor drag on heat transfer are seen to be quite small. At the 
fixed pu ratio of 150, larger effects were found for either thin films 
(c, AT /ht, — 0) or very thick films (c,AT/hy = 0.855). The 
change of the py ratio from 150 to 50 (at a fixed value of c,AT/ 
hy, and Pr) had almost no effect on the case considered, although 
larger effects might conceivably have occurred at the extremes 
of thick and thin films. For the single liquid metal case con- 
sidered (Pr = 0.008), the effect of vapor drag was larger than 
that for the nonmetallic liquids. This finding is associated with 
the fact that the liquid metal case had the largest (dimension- 
less) film thickness of all the entries of the table. Apparently, 
this was sufficient to override the expected tendency of the high 
pu ratio (650) to diminish the coupling between liquid and vapor. 

The experiment of reference [3] corresponded approximately 
to the second entry in the table. Since the effect of vapor drag 
on the analytical prediction of heat transfer as given in reference 
{2|is not significant, the 30 per cent deviation between theory 
and experiment must be laid to other causes. 
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Table | Effects of vapor drag 


cp AT Pr No drag? Drag No drag? Drag 
150 All 0.903 0.9147 
150 0.170 10 0.924 0.9176 0.041 0.08005 
150 0.855 10 0.995 0.9654 0.126 0.2174 
50 0.172 10 0.924* 0.9208 0.041* 0.08283 
650 0.00244 0 0.31 0.4627 


* Thin film theory, Appendix, reference [2]. 
+ Interpolated values from reference [2]. 
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Radiation Fin Effectiveness 


J. G. BARTAS' and W. H. SELLERS? 
Introduction 


POWER PLANTs operating in space for long periods of time have 
one heat sink, space, and one means of transferring heat to the 
heat sink, a radiator. This radiator is a “‘true’’ radiator in that 
all heat leaving its surfaces does so by radiation. There is, ap- 
parently, no reference in the available literature which provides 
a rapid method for estimating the radiator effectiveness and for 
optimizing fin designs. 

This analysis considers the radiator configuration shown in Fig. 
1. In this design, parallel pipes are joined by webs, which act 
as radiation fins. Heat flows by conduction from the pipes down 
the fin and radiates from both surfaces. The fin surface may ‘‘see’’ 
other portions of the radiator if the radiator is conical, cylindrical, 
spherical, and so forth. The radiator may see a sun or other 
body, as well as deep space. The absorbed incident energy 
from these sources can be determined, and a mean “‘space’’ tem- 
perature 7, can be computed. Thus the problem is reduced to 
consideration of a fin which radiates to a constant temperature 
environment 

Results for radiator fin effectiveness have been computed for a 
limiting class of radiation fins—the case where the radiator is per- 
fectly flat, i.e., there is no spacial interference of radiation by the 
pipes. This situation also is mathematically equivalent to a 
condition where the pipes are perfect specular reflectors. 


Assumptions 
The following assumptions are made: 


1 Material properties k, € are independent of the tempera- 
ture. 

2 Emissivity € is independent of the angle. 

3 The fin thickness 2) is constant. 

4 The fin base temperature 7’, is constant across the base of 
the fin. 

! Senior Heat Transfer Engineer, General Electric Company, West 
Lynn, Mass. Assoc. Mem. ASME. 

* Supervisor, Calculations and Special Problems Unit, General 
Electric Company, West Lynn, Mass. 

Contributed by the Heat Transfer Division of Tae American 
Society or MecHanicat ENGINEERS and presented at the ASME- 
AIChE Heat Transfer Conference Storrs, Conn., August 9-12, 1959. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, 
November 26, 1959. Condensed from paper No. 59—-HT-17. 
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5 The fin base temperature 7’, is constant along the pipe axis. 
Conduction heat flow is one-dimensional. 
6 The surface heat transfer is entirely by radiation. 


The foregoing assumptions are those made in conventional fin 
analyses, except for the added restrictions due to the radiation 
effects and the lack of convection. 


Energy Balance Equation 
Consider the geometry shown in Fig. 2. The energy balance 
equation for the element dz is 


Q. + = 0. + + Qrovt (1) 


Conduction heat flow Q, is 


where the cross-sectional flow area A, is equal to 2b per unit fin 
depth. The net radiant heat flux is given by 


Qrovr — Qain = — 
where é is the emissivity function, and dA, is the incremental sur- 
face area 2dr. The parameter 7'gz is defined here as 
T sr‘ Ts(1 F) + T,'F 


where 7’, is the space temperature, 7’, is the fin root temperature, 
and F is the radiant interchange configuration factor of the ele- 
ment surface dz with respect to the two nearest pipes. The net 
radiant heat flux equation is, therefore, 


Qrovt — = oedA,[T* — — F) — (3) 


Combining Equations (1), (2), (3), and canceling common 
terms, we have 


aT 
kb — — F) — (4) 


zr? 


Equation (4) is subject to the boundary conditions 
at 
dT 
=0 — 
dr 
at 


z= L, T = Tr 
Equation (4) is nondimensionalized to 
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Fig. 3 Radiator fin effectiveness versus parameter \ 
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d*6 
dp? = — 0,4 — F(1 — 
where 
6 = T/T, 
p = 2/L 
X = o&L*T,*/kb 


The boundary conditions are 
at 


p =0, 


at 
p=1, 
Effectiveness 
The fin effectiveness 7 is defined as 


(8) 


where Q,, and Q, are the actual and ideal heat flows, respectively. 
The actual heat flux is evaluated at z = L 


1 


The ideal heat flux is defined by assuming the entire fin is at the 
fin root temperature and radiates without interference, 


(9) 


Combining Equations (8), (9), (10) 


(3) 
dp p=l 


(11) 


Results 

A limiting case of radiator design has been analyzed in detail- 
It was assumed that F = 0. That is, there is no special inter- 
ference of the radiation by the pipes. 

Equations (6) and (11) were solved numerically on an IBM 704 


digital computer. The results for the fin effectiveness are 
plotted in Fig. 3. 


Optimum Fin Geometry 
The premium for weight saving in space hardware makes it im- 
perative to optimize radiator designs. The following is limited 
to finding the optimum design for a rectangular fin such as that 
shown in Fig. 2. We assume that fin weight remains constant, 
while heat flow Q, is maximized. 
With F equal to zero, combine Equation (8), (10), and solve for 
Q, = — 8,4) 
Since fin weight is constant, we have 
bL = c’ (13) 
The definition of \, combined with Equations (12), (13), leads to 


(14) 


(12) 


Q, = constant 
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where the constant is equal to 
2T 1 — 


We can plot Q, versus A and find the maximum Q, directly, 9 
being a function of \ as shown in Fig. 3. 

A graph of the effectiveness 7 and parameter A that correspond 
to a minimum weight fin at any @, is plotted in Fig. 4. 


Concluding Comment 
Note the relationship between the present parameter \ and the 
parameter ml in conventional convection fin analyses. 


kb 


L? 
(ml? = 


It is apparent that 
h ~ 


This result is expected, since the equivalent radiation heat-trans- 
fer coefficient is commonly defined as 


_h, = 


if T's — 0. 

The radiation fin performance curve follows the classical pat- 
tern but with considerable variation due to the mode of heat 
transfer. It is unfortunate that the profuse data available on 


convective fins are not convertible to the evaluation of radiation 
fins. 


A Tabulation of the Function 
exp erfe X 


J. E. SUNDERLAND' and R. J. GROSH? 


THe FUNCTION exp X?* erfe X is tabulated for values of the argu- 
ment running from 0.10 to 7.00 in increments of 0.01. The 
values reported were obtained by hand and machine computation 
and are given to four, five, or six significant figures with an error 
of plus or minus one in the last figure. For values of the argu- 
ment less than 0.10 or greater than 7.00, convenient expressions 
are given for the function. The tabulation of the function has 
been placed on file at the American Documentation Institute.* 

The function exp X? erfe X often appears in the solution of 
transient heat conduction and diffusion problems. Carslaw and 
Jaeger [1]* give a table for values of the argument ranging from 
0.00 to 3.00 in increments of 0.05 or 0.10. The values are reported 

1 Assistant Professor of Mechanical Engineering, Northwestern 
University, Evanston, Illinois. 

2 Professor of Mechanical Engineering, Purdue University, Lafay- 
ette, Ind. 

* The tabulation of the function has been deposited as Document 
number 6157 with the ADI Auxiliary Publications Project, Photo- 
duplication Service, Library of Congress, Washington 25, D.C. A 
copy may be secured by citing the Document number and by re- 
mitting $1.25 for photoprints, for $1.25 for 35 mm microfilm. Advance 
payment is required. Make checks or money orders payable to: 
Chief, Photoduplication Service, Library of Congress. 

* Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of Tae AMERICAN 
Society or MEecHANICAL ENGINEERS. 

Norte: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, Decem- 
ber 9, 1959. 
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with four significant figures. Radau [2] tabulated 1/2/2 exp X? 
erfe X to twelve significant figures for integral values of X* from 
3 to 16. The function 1/#/2 exp X*/2 erfe X/+/2 has been 
tabulated by Sheppard [3] to twenty-four places for values of X 
ranging from 0.00 to 10.0 in increments of 0.01, and by Mills [4] 
to five places for values of X ranging from 0.00 to 4.00 in incre- 
ments of 0.01, 4.00 to 5.00 in increments of 0.05, and 5.00 to 10.0 
in increments of 0.10. All these functions are, of course, related. 
The purpose of the tabulation reported here is to provide a con- 
venient and accurate table of the function exp X* erfe X. Before 
the results are given, the functions will be defined, the method of 
calculation described, and certain useful, but not original ap- 
proximations indicated. 

The symbol exp X* is defined as the value of the number 2.71828 
(e) raised to the X* power. 

The symbol erfe X is defined as the complement of the error 
function with argument X, erf X; that is 


erfe X = 1 — erf X 


where 


2 z 
erf X = we f exp (—u?)du 


It follows that 


erfe X = +f exp (—1*)du 


Finally, it can be remarked that these integrals are closely related 
to the probability or Gauss error integral [5, 6]. 

Rather than make use of the previous tabulations, which differ 
from the results reported here by constants in the argument 
and/or preceding the function, it was found more convenient to 
use an ElectroData digital computer for obtaining values of the 
function when X was in the range 0.10 to 3.00. The programming 
routine employed polynomial expansions for the exponential and 
error functions. 

For values of X in the range 3.00 to 7.00 the function was 
evaluated by means of the twelve to eighteen place tables given in 
references [7 and 8]. 

The consistency of the table was checked by a simple differenc- 
ing technique which involved taking the first and second dif- 
ferences between each successive number on the table. Where the 
second differences showed no abrupt changes, the table was con- 
sidered to be consistent. 

For values of X greater than one the function can be deter- 
mined from the expansion 


16 / 1960 


1 13 


1.3.5... (2n — 


and the error is less than the last term used. As X becomes larger, 
the series converges faster and for values of X greater than seven, 
the error is less than 1 per cent when only the first term of the ex- 
pansion is used. Using two terms of the expansion, the error is 
less than 0.04 per cent. 

For all finite values of X it is noted that 


(2n + 1)n! 


f(z) = exp (2?) E 


For values of X less than 0.1, using only the first term of the 
series expansion, the function becomes 


f(z) = exp (z*) 


and gives an error of less than 0.05 per cent. 
The tabulation of the function is on file at the American 
Documentation Institute. 
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STRUCTURAL 
DAMPING 


In this book engineers and scien- 
tists will find a wealth of information 
useful in engineering design, such 
as: 


the discussion concerning the 
principal mechanisms by which 
structure dissipates damping 
energy; 


the identification and classifica- 
tion of the various types of ma- 
terial and system mechanisms; 


the reviews of current research 
in energy dissipation through in- 
terfacial slip at contact surfaces 
of structure and machine joints; 


a method for analyzing systems 
composed of viscoelastic and 
elastic materials applied to an 
elastic plate, along with specific 
results for a number of damping 
treatments such as, automotive 
undercoat, damping tape, and 
sandwich plates with dissipative 
cores; 


the results of a preliminary ex- 
perimental study of vibrational 
energy dissipation at the support 
junctions of built-in beams and 
plates; 


details of structural damping 
measurement methods and the 
limitations that are imposed by 
the assumptions associated with 
each method; 


a fatigue strength damping cri- 
teria for selecting the proper 
structural materials for the de- 
sign of members which are sub- 
jected to fluctuating resonant 
stresses; 


a hundred references to articles 
on structural damping, arranged 
alphabetically according to 
authors and chronologically for a 
given author. 


This information was prepared by 
authorities in the field and presented 
at a colloquium held at the 1959 
ASME Annual Meeting. 
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SPRAY 
LITERATURE 
ABSTRACTS 


Here, for the first time, is a com- 
prehensive source-book on the 
multitudinous applications § of 
sprays and their scientific, tech- 
nical and industrial aspects. The 
compilation comprises over 1300 
items collected from the technical 
literature of 24 countries, published 
from 1880 to 1958 inclusive. The 
abstracts are detailed, contain full 
bibliographic data, as well as in- 
formation on the affiliation of each 
author wherever available. Entries 
include: 


fuel sprays for furnaces, internal com- 
bustion engines, gas turbines, and 
rockets; 

sprays in industrial processes for evap- 
oration, drying, humidification, cool- 
ing, air conditioning, and chemical 
reactions; 

atmospheric sprays, as rain, fog, haze, 
sleet and snow, their ‘formation, 
evaporation and freezing; 

agricultural sprays for sprinkling, pest, 
fungus, and weed control; 

fire- ‘Ss, as in sprinkler 
systems and int ire-smothering fogs; 

aerosols for medical, disinfecting, ther- 
apeutic, and air conditioning pur- 
poses; 

sprays and fogs for ~ “gee purposesin 
defense and attack 


Other entries direct attention to 
published information having a 
bearing on the theoretical and 
experimental study of sprays; in 
particular: 


mathematics, from arithmetics 
to statistics and probability theory; 

« - ics and chemistry, dealing with 

e pesgerses of liquids and gases, 

especially the science of colloids and 
other disperse systems; 

hydro-, aero-, and thermodynamics, 
and rheology, dealing with the 
motion of liquids in spray-generating 
devices, the breakup of liquid stream 
into drops, and their motion in air, 
their evaporation and combustion; 

optical and electronic devices for the 
experimental study of sprays; 

branches of mechanical and produc- 
tion engineering concerned with the 
design, production, and testing of 
pumps, nozzles, and other elements 
of spray generating equipment; 

borderline and associated disciplines, 
such as the science and technology 
of powders and dusts. 
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TECHNICAL 
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This is a comprehensive treat- 
ment that can be used profitably 
by engineers, scientists, and tech- 
nicians who wish to improve their 
skill and effectiveness as well as 
by those who have never written 
before. It shows how to look for 
ideas, how to evaluate them, how 
to build an idea into an outline, 
how to write-up the idea, and how 
to work with editors and publishers 
in getting the idea into print. In 
addition to sections on technical 
articles, reports, and books, includ- 
ing text books and handbooks, 
this book also gives useful material 
on military and industrial training 
manuals and industrial advertising. 
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McGraw-Hill Book Company, ASME 


members may Guschese cop es from 
the Society at a 20% discount. 
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With the help of this Ameri- 
can Standard clear, legible, and 
effective illustrations can be pre- 
pared with a minimum of effort. 
Here will be found descriptions of 
preferred illustrations and layout 
practices, a discussion of the 
factors influencing legibility of 
reproductions, suggestions on the 
preparation of illustrations suitable 
for technical publications, recom- 
mendations on preparing copy for 
either still projections or publica- 
tions, and general guidance infor- 
mation on drafting practices and on 
material useful in illustration work. 
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